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FILLMORE SPRINGER CNOPS CONSTRUCTION 
IMPLEMENTED IN GilMaC 

VLADIMIR V. KISIL 


Abstract. This is an implementation of the Fill more Springer—Cnops construction (FSCc) based on the Clifford 
algebra capacities [11] of the GiNaC computer algebra system. FSCc linearises the linear-fraction action of the Mobius 
group. This turns to be very useful in several theoretical and applied fields including engineering. 

The core of this realisation of FSCc is done for an arbitrary dimension, while a subclass for two dimensional cycles 
add some 2D-specific routines including a visualisation to PostScript files through the MetaPost or Asymptote software. 

This library is a backbone of many result published in [10], which serve as illustrations of its usage. It can be ported 
(with various level of required changes) to other CAS with Clifford algebras capabilities. 
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1. Introduction 

The usage of computer algebra system (CAS) in Clifford Algebra research has an established history with the 
famous “Green book” [3] already accompanied by a floppy disk with a REDUCE package. This tradition is very much 
alive, see for example the recent book [6] accompanied by a software CD. Numerous new packages are developed by 
various research teams across the world to work with Clifford algebras generally or address specific tasks, see on-line 
proceedings of the recent IKM-2006 conference [5]. 

Along this lines the present paper presents an implementation of the Fillmore-Springer-Cnops construction (FSCc) 
along with illustrations of its usage. FSCc [2,13] linearises the linear-fraction action of the Mdbius group in R". This 
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has clear advantages in several theoretical and applied fields including engineering. Our implementation is based on 
the Clifford algebra capacities of the GiNaC computer algebra system]!], which were described in [ ]. 

The core of this realisation of FSCc is done for an arbitrary dimension of R” with a metric given by an arbitrary 
bilinear form. We also present a subclass for two dimensional cycles (i.e. circles, parabolas and hyperbolas), which 
add some 2D specific routines including a visualisation to PostScript hies through the MetaPost [8] or Asymptote [7] 
packages. This software is the backbone of many results published in [10] and we use its application to [10] for the 
demonstration purpose. 

There is a Python wrapper [9] for this library. It is based on BoostPython and pyGiNaC packages. The wrapper 
allows to use all functions and methods from the library in Python scripts or Python interactive shell. The drawing 
of object from cycle2D may be instantly seen in the interactive mode through the Asymptote. 

The present package can be ported (with various level of required changes) to other CAS with Clifford algebras 
capabilities similar to GiNaC. 


2. User interface to classes cycle and cycle2D 
The cycle class describes loci of points x £ R" defined by a quadratic equation 
(2.1) fcx 2 — 2 (1, x) + m = 0, where k, m £ R, 1 £ R". 

The class cycle correspondingly has member variables k , l , m to describe the equation (2.1) and the Clifford algebra 
unit to describe the metric of surrounding space. The plenty of methods are supplied for various tasks within FSCc. 
We also define a subclass cycle2D which has more methods specific to two dimensional environment. 

2.1. Constructors of cycle. Here is various constructors for the cycles. The first one takes values of k, 1, m as well 
as metric supplied directly. Note that 1 is admitted either in form of a 1st, matrix or indexed objects from GiNaC. 
Similarly metric can be given by an object from either tensor, indexed, matrix or Clifford classes exactly in the 
same way as metric is provided for a clifford-unit() constructors [11]. 

2a (cycle class constructors 2a) = (38b) 2b o 

public: 

cycle(const ex & kl, const ex & 11, const ex & ml, 

const ex & metr = -(new tensdelta)^>setflag(status-flags::dynallocated)); 

Defines: 

cycle, used in chunks 2-9, 11-17, 20e, 24, 25, 33a, 38-47, 49-56, and 58. 

Uses ex 4a lib 39d 48a 48b 66a 66b 66c 67a. 

Constructor for a cycle (2.1) with k = 1 and given l defined by the condition that square of its “radius” (which is 
detC, see [10, Defn. 5.1]) is rsquared. Note that for the default value of the metric the value of l coincides with the 
centre of this cycle. 

2b (cycle class constructors 2a) += (38b) <2a 2c > 

cycle(const 1st & /, 

const ex & metr = -(new tensdelta)-^setflag(status-flags::dynallocated), 
const ex & rsquared = 0, const ex & e = 0, 

const ex & sign = (new tensdelta)^setflag(status-flags::dynallocated)); 


Defines: 

cycle, used in chunks 2-9, 11-17, 20e, 24, 25, 33a, 38-47, 49-56, and 58. 

Uses ex 4a lib 39d 48a 48b 66a 66b 66c 67a and 1st lie. 

If we want to have a cycle identical to to a given one C up to a space metric which should be replaced by a new one 
metr , we can use the next constructor. 

2c (cycle class constructors 2a) += (38b) «2b 2do 

cycle(const cycle & C, const ex & metr); 


Uses cycle 2a 2b 3e 3e 3e 3e 3e 46a 46b 47b 47b 48a 48a 48a 48a 48a and ex 4a lib 39d 48a 48b 66a 66b 66c 67a. 

To any cycle FSCc associates a matrix, which is of the form (2.2) [10, (3.2)]. The following constructor make a cycle 
from its matrix representation, i.e. it is the realisation of the inverse of the map Q [10, (3.2)]. 

2d (cycle class constructors 2a) += (38b) «2c 

cycle(const matrix & M, const ex & metr, const ex & e = 0, const ex & sign = 0); 


Uses cycle 2a 2b 3e 3e 3e 3e 3e 46a 46b 47b 47b 48a 48a 48a 48a 48a, ex 4a lib 39d 48a 48b 66a 66b 66c 67a, and matrix 8c lOe Ilf 12a 25a. 


2.2. Accessing parameters of a cycle. The following set of methods get-*{) provide a reading access to the various 
data in the class. 

3a (accessing the data of a cycle 3a) = (38b) 3b o 

public: 

virtual inline ex get-dim{ ) const { return ex-to<idx>{unit.op{l)).get_dim {); } 
virtual inline ex get-metric() const { return ex-to<clifford>(unit).get-metric (); } 
virtual inline ex gef_metnc(const ex SziO , const ex fail) const 
{ return exAo<c\idord> {unit), getjmetric{iO, il); } 

virtual inline ex get-k{) const { return k; } 


Uses ex 4a lib 39d 48a 48b 66a 66b 66c 67a. 

The member l can be obtained as the whole by the call getJQ, or its individual component is read, for example, by 
getJ(l). 

3b (accessing the data of a cycle 3a) += (38b) «3a 3c o 

inline ex getJQ const { return l; } 
inline ex getJ{ const ex & i) const 
{ return l.subs{l.op{ 1) = i, subs-options::no-pattern ); } 
inline ex get-m{) const {return m;} 
inline ex get-unit{) const {return unit ;} 


Uses ex 4a lib 39d 48a 48b 66a 66b 66c 67a. 

Methods nops {), op{), let-op {), is-equal{), subs{) are standard for expression in GiNaC and described in the GiNaC 
tutorial. The first three methods are rarely called by a user. In many cases the method subs{) may replaced by more 
suitable subject-to{) 2.4. 

3c (accessing the data of a cycle 3a) += (38b) «3b 

sizeJ, nops{) const {return 4;} 
ex op{size-t i) const; 
ex & let-op{size-t i); 

bool is-equal( const basic & other ) const; 
bool is-zero{) const; 

cycle sw6s(const ex & e, unsigned options = 0) const; 
inline cycle normalQ const 

{ return cycle(k.normal{), LnormalQ, m.normal{) 7 unit.normal{));y 

Uses bool lie, cycle 2a 2b 3e 3e 3e 3e 3e 46a 46b 47b 47b 48a 48a 48a 48a 48a, and ex 4a lib 39d 48a 48b 66a 66b 66c 67a. 

2.3. Linear Operations on Cycles. Cycles are represented by a points in a projective vector space, thus we wish 
to have a full set of linear operation on them. The metric is inherited from the first cycle object. First we define it 
as an methods of the cycle class. 

3d (Linear operation as cycle methods 3d)= (38b) 

cycle add{ const cycle & rh) const; 
cycle sw&(const cycle & rh) const; 
cycle eOTW^const ex & rh) const; 
cycle div{ const ex & rh) const; 


Uses cycle 2a 2b 3e 3e 3e 3e 3e 46a 46b 47b 47b 48a 48a 48a 48a 48a and ex 4a lib 39d 48a 48b 66a 66b 66c 67a. 
After that we overload standard binary operations for cycle. 

3e (Linear operation on cycles 3e)= (38b) 4a c> 

const cycle operator+(const cycle & Ih, const cycle & rh); 
const cycle operator-(const cycle & Ih, const cycle & rh); 
const cycle operator*(const cycle & Ih, const ex & rh); 
const cycle operator*(const ex & Ih, const cycle & rh); 
const cycle operator-^(const cycle & Ih, const ex & rh); 


Defines: 

cycle, used in chunks 2-9, 11-17, 20e, 24, 25, 33a, 38-47, 49-56, and 58. 
Uses ex 4a lib 39d 48a 48b 66a 66b 66c 67a. 


We also define a product of two cycles through their matrix representation (2.2). 

4a (Linear operation on cycles 3e)+= (38b) «3e 

const ex operator* (const cycle & Ih, const cycle & r/i); 

Defines: 

ex, used in chunks 2-8, 11a, 13, 14, 17, 19-25, 33b, 38-40, 42, 44-47, 49-60, 63, 65a, 67, and 68. 

Uses cycle 2a 2b 3e 3e 3e 3e 3e 46a 46b 47b 47b 48a 48a 48a 48a 48a. 

2.4. Geometric methods in cycle. We start from some general methods which deal with cycle. The next method 
is needed to get rid of the homogeneous ambiguity in the projective space of cycles. If k_new= 0 the cycle is normalised 
such that its det becomes 1. Otherwise the first non-zero coefficient among k, m, Iq, l±, ... is set to k_new. 

4b (specific methods of the class cycle 4b) = (38b) 4c > 

public: 

cycle normalize( const ex & k-new = numeric(l), const ex & e = 0) const; 
cycle normalize-det(const ex & e = 0) const; 

Uses cycle 2a 2b 3e 3e 3e 3e 3e 46a 46b 47b 47b 48a 48a 48a 48a 48a, ex 4a lib 39d 48a 48b 66a 66b 66c 67a, and numeric lOd. 

The method center () returns a list of components of the cycle centre or the corresponding vector (Hxl matrix) if the 
dimension is not symbolic. The metric, if not supplied is taken from the cycle. 

4c (specific methods of the class cycle 4b) += (38b) a 4b 4d > 

virtual ex center( const ex & metr = 0, bool returnjmatrix = false) const; 

Uses bool lie and ex 4a lib 39d 48a 48b 66a 66b 66c 67a. 

The next method returns the value of the expression —fcx 2 — 2 (1, x) + m for the given cycle and point x. Obviously 
it should be 0 if x belongs to the cycle. 

4d (specific methods of the class cycle 4b)+= (38b) a4c 4e> 

virtual ex val(const ex & y) const; 

Uses ex 4a lib 39d 48a 48b 66a 66b 66c 67a. 

Then method passing () returns a relational defined by the identity fcx 2 — 2 (1, x) + m = 0, i.e this relational describes 
incidence of point to a cycle. 

4e (specific methods of the class cycle 4b)+= (38b) «4d 4f> 

inline relational passing^ const ex & y) const {return val(y).numer().normal() = 0;> 

Uses ex 4a lib 39d 48a 48b 66a 66b 66c 67a and relational lie. 

We oftenly need to consider a cycle which satisfies some additional conditions, this can be done by the following 
method subject-to. Its typical application looks like: 

Cl = C.subject-to(lst(C.passing(P ), C.is-orthogonal(Cl ))); 

The second parameters vars specifies which components of the cycle are considered as unknown. Its default value 
represents all of them which are symbols. 

4f (specific methods of the class cycle 4b) += (38b) <4e 4go 

cycle subject-to( const ex & condition , const ex & vars = 0) const; 

Uses cycle 2a 2b 3e 3e 3e 3e 3e 46a 46b 47b 47b 48a 48a 48a 48a 48a and ex 4a lib 39d 48a 48b 66a 66b 66c 67a. 

2.5. Methods representing FSCc. There is a set of specific methods which represent mathematical side of FSCc. 
The next method is the main gateway to the FSCc, it generates the 2x2 matrix 

(2.2) ^ l< fc^ ba^e ij fr° m the cycle fcx 2 — 2 (1, x) + m = 0. 

Note, that the Clifford unit e has an arbitrary metric unrelated to the initial metric stored in the unit member variable. 

4g (specific methods of the class cycle 4b) += (38b) <4f 5a o 

virtual matrix to-matrix( const ex & e = 0, 

const ex & sign = (new tensdelta)^>setflag(status-flags::dynallocated )) const; 


Uses ex 4a lib 39d 48a 48b 66a 66b 66c 67a and matrix 8c lOe Ilf 12a 25a. 


The next method returns the value of determinant of the matrix (2.2) corresponding to the cycle. It has explicit 
geometric meaning, see [10, § 5.1]. Before calculation the cycle is normalised by the condition k^k-norm, if k_norm is 
zero then no normalisation is done. 

(specific methods of the class cycle 4b) += (38b) <4g 5b o 

virtual ex def(const ex & e = 0, 

const ex & sign = (new tensdelta)^setflag(status-flags::dynallocated ), 
const ex & k-norm = numeric(l)) const; 


Uses ex 4a lib 39d 48a 48b 66a 66b 66c 67a and numeric lOd. 

The matrix (2.2) corresponding to a cycle may be multiplied by another matrix, which in turn may be either generated 
by another cycle or be of a different origin. The next methods multiplies a cycle by another cycle of matrix supplied 
in C. 

(specific methods of the class cycle 4b)+= (38b) <5a 5c o 

virtual matrix mul{ const ex & C, const ex & e = 0, 

const ex & sign = (new tensdelta)-^setflag(status-flags::dynallocated ), 
const ex & signl =0) const; 

Uses ex 4a lib 39d 48a 48b 66a 66b 66c 67a and matrix 8c lOe Ilf 12a 25a. 

Having a matrix C which represents a cycle and another matrix M we can consider a similar matrix M~ X CM. The 
later matrix will correspond to a cycle as well, which may be obtained by the following three methods. In the case 
then M belongs to the SL 2 {R) group the next two methods make a proper conversion of M into Clifford-valued form. 

(specific methods of the class cycle 4b) += (38b) «5b 5do 

cycle sl2similarity(const ex & a, const ex & b, const ex & c, const ex & d, 
const ex & e = 0, 

const ex & sign = (new tensdelta)^setflag(status-flags::dynallocated ), 
bool not-inverse= true) const; 
cycle sl2similarity( const ex & M, const ex & e = 0, 

const ex & sign = (new tensdelta)^>setflag(status-flags::dynallocated), 
bool not-inverse= true) const; 

Uses bool lie, cycle 2a 2b 3e 3e 3e 3e 3e 46a 46b 47b 47b 48a 48a 48a 48a 48a, and ex 4a lib 39d 48a 48b 66a 66b 66c 67a. 

If M is a generic 2 x 2-matrix of another sort then it is used in the similarity in the unchanged form by the next 
method. 

(specific methods of the class cycle 4b)+= (38b) <5c 5ei> 

cycle matrixsimilarity( const ex & M, const ex & e = 0, 

const ex & sign = (new tensdelta)—>setflag(status-flags::dynallocated), 
bool not-inverse= true) const; 


Uses bool lie, cycle 2a 2b 3e 3e 3e 3e 3e 46a 46b 47b 47b 48a 48a 48a 48a 48a, and ex 4a lib 39d 48a 48b 66a 66b 66c 67a. 


The 2 x 2-matrix M = 


\ 

r I can be also defined by the collection of its elements. 


(specific methods of the class cycle 4b) += (38b) <5d 5f> 

cycle matrixsimilarity( const ex & a, const ex & b, const ex & c, const ex & d, 
const ex & e = 0, 

const ex & sign = (new tensdelta)^setflag(status-flags::dynallocated), 
bool not_inverse= true) const; 


Uses bool lie, cycle 2a 2b 3e 3e 3e 3e 3e 46a 46b 47b 47b 48a 48a 48a 48a 48a, and ex 4a lib 39d 48a 48b 66a 66b 66c 67a. 

Finally, we have a method for reflection of a cycle in another cycle (7, which is given by the similarity of the representing 
matrices: CC\C, see [10, § 4.2]. 

(specific methods of the class cycle 4b) += (38b) <5e 6a o 

cycle cycle-similarity (const cycle & C, const ex & e = 0, 

const ex & sign = (new tensdelta)^>setflag(status-flags::dynallocated), 
const ex & signl = 0) const; 


Uses cycle 2a 2b 3e 3e 3e 3e 3e 46a 46b 47b 47b 48a 48a 48a 48a 48a and ex 4a lib 39d 48a 48b 66a 66b 66c 67a. 


A cycle in the matrix form (2.2) naturally defines a Mobius transformations of the points: 


(2.3) 


m \ 

V k -^)&J 


The following methods realised this transformations. 


ljcrje^x + TO 
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(specific methods of the class cycle 4b) += (38b) <5f 6b o 

virtual inline ex moebius-map( const ex & P, const ex & e = 0, 

const ex & sign = (new tensdelta)^>setflag(status-flags::dynallocated )) const 
{return clifford-moebius-map(to-matrix(e, sign), P, (e.is-zero()7get-metric():e))\} 


Uses ex 4a lib 39d 48a 48b 66a 66b 66c 67a. 

For two matrices C\ and C 2 obtained from cycles the expression 

(2.4) {C 1 ,C 2 )=mr{C l C 2 ) 

naturally defines an inner product in the space of cycles. The follwong methods realised it. 

(specific methods of the class cycle 4b)+= (38b) <6a 6c o 

virtual inline ex inner_product( const cycle & C, const ex & e = 0, 

const ex & sign = (new tensdelta)^setflag(status-flags::dynallocated)) const 
{return scalar_part(mul(C, e, sign) .trace())-^ numeric(2);} 


Uses cycle 2a 2b 3e 3e 3e 3e 3e 46a 46b 47b 47b 48a 48a 48a 48a 48a, ex 4a lib 39d 48a 48b 66a 66b 66c 67a, and numeric lOd. 

The inner product (2.4) defines an orthogonality relation (C\, C 2 ) = 0 in the space of cycles which returned by the 
method is-orthogonal(). 

(specific methods of the class cycle 4b) += (38b) <1 6b 6d > 

virtual inline relational is-orthogonal( const cycle & G, const ex & e = 0, 
const ex & sign = (new tensdelta)^>setflag(status-flags::dynallocated)) const 
{return ( inner-product(C, e, sign) = 0);} 


Uses cycle 2a 2b 3e 3e 3e 3e 3e 46a 46b 47b 47b 48a 48a 48a 48a 48a, ex 4a lib 39d 48a 48b 66a 66b 66c 67a, and relational lie. 

In many cases we need a higher order orthogonal relation between cycles— so called s-orthogonality, see [10, § 4.3], 
which is given by the relation: 

mr{ClCiClRl)=Q. 

(specific methods of the class cycle 4b)+= (38b) <6c 6eo 

relational iss-orthogonal( const cycle & C, const ex & e = 0, 

const ex & sign = (new tensdelta)^>setflag(status-flags::dynallocated)) const; 

Uses cycle 2a 2b 3e 3e 3e 3e 3e 46a 46b 47b 47b 48a 48a 48a 48a 48a, ex 4a lib 39d 48a 48b 66a 66b 66c 67a, and relational lie. 

The remaining to methods check if a cycle is a liner object and if it is normalised to k = 1. 

(specific methods of the class cycle 4b) += (38b) <6d 

inline relational isJinearQ const {return (k = 0);} 
inline relational is-normalized() const {return (k = 1);} 


Uses relational lie. 

2.6. Two dimensional cycles. Two dimensional cycle cycle2D is a derived class of cycle. We need to add only 
very few specific methods for two dimensions, notably for the visualisation. 

This a specialisation of the constructors from cycle class to cycle2D. Here is the main constructor. 

(constructors of the class cycle2D 6f)= (39b) 6go 

public: 

cycle2D(const ex & kl , const ex & 11, const ex & ml, 

const ex & metr = -unit-matrix{2))\ 

Uses cycle2D 6g 11c 11c 34b 35b 39d 39d 39d 39d 39d 48b 48b 48b 48b 48b 56b 56b 57d and ex 4a lib 39d 48a 48b 66a 66b 66c 67a. 
Constructor for the cycle from / and square of its radius. 

(constructors of the class cycle2D 6f)+= (39b) o6f 7ao 

cycle2D(const 1st & l, const ex & metr = -unitjmatrix(2), const ex & msquared =0, 
const ex & e =0, const ex & sign = unitjmatrix(2))\ 

Defines: 

cycle2D, used in chunks 6f, 7a, 12c, 14-18, 20a, 21g, 24a, 25b, 30a, 32, 33, 35-37, 39-41, 55-57, and 62-64. 



Make a two dimensional cycle out of a general one, if the dimensionality of the space permits. The metric of point 
space can be replaced as well if a valid metr is supplied. 

(constructors of the class cycle2D 6f}+= (39b) «6g 

cycle2D( const cycle & C, const ex & metr = 0); 

Uses cycle 2a 2b 3e 3e 3e 3e 3e 46a 46b 47b 47b 48a 48a 48a 48a 48a, cycle2D 6g 11c 11c 34b 35b 39d 39d 39d 39d 39d 48b 48b 48b 48b 
48b 56b 56b 57d, and ex 4a lib 39d 48a 48b 66a 66b 66c 67a. 

The method focusfi) returns list of the focus coordinates and the focal length is provided by focaLlengthfi). This turns 
to be meaningful not only for parabolas, see [10]. 

(methods specific for class cycle2D 7b) = (39b) 7c o 

public: 

ex focus( const ex & e = diag-matrix(\st{-l, 1)), bool return-matrix = false) const; 
inline ex focaLlengthQ const {return (gefJ(l)-P2-bfc);} // focal length of the cycle 


Uses bool lie, cycle 2a 2b 3e 3e 3e 3e 3e 46a 46b 47b 47b 48a 48a 48a 48a 48a, ex 4a lib 39d 48a 48b 66a 66b 66c 67a, and 1st lie. 
The methods rootsQ returns values of u (if first true) such that k(u 2 + ey 2 ) — 2l\u — 2l2y + m = 0, i.e. solves a 
quadratic equations. If first = false then values of v satisfying to k{y 2 + ev 2 ) — 2l\y — 2l 2 v + m = 0 are returned, 
(methods specific for class cycle2D 7b) += (39b) <7b 7do 

1st roofs(const ex & y = 0, bool first = true) const; 

Uses bool lie, ex 4a lib 39d 48a 48b 66a 66b 66c 67a, and 1st lie. 

The next methods is a generalisation of the previous one: it returns intersection points with the line ax + b. 

(methods specific for class cycle2D 7b) += (39b) <7c 7eo 

1st line-intersect (const ex & a, const ex & b) const; 

Uses ex 4a lib 39d 48a 48b 66a 66b 66c 67a and 1st lie. 

The method metapost-drawfi) outputs to the stream ost MetaPost comands to draw parts of two the cycle2D within 
the rectangle with the lower left vertex ( xmin, ymin) and upper right ( xmax, ymax). The colour of drawing is specified 
by color (the default is black) and any additional MetaPostoptions can be provided in the string more-options. By 
default each set of the drawing commands is preceded a comment line giving description of the cycle, this can be 
suppressed by setting with-header = false. The default number of points per arc is reasonable in most cases, however 
user can override this with supplying a value to points-per-arc. The last parameter is for internal use. 

(methods specific for class cycle2D 7b) += (39b) «7d 7f> 

void metapost-draw(ostream & ost, const ex & xmin = -5, const ex & xmax = 5, 
const ex & ymin = -5, const ex & ymax = 5, const 1st & color = lst(), 
const char * more-options = "", 

bool with-header = true, int points-per-arc = 0, bool asymptote = false, 
char * picture = bool only-path= false, bool is-Continuation=i&lse) const; 

Uses bool lie, ex 4a lib 39d 48a 48b 66a 66b 66c 67a, and 1st lie. 

The similar method provides a drawing output for Asymptote [7] with the same meaning of parameters. However 
format of more-options should be adjusted correspondingly. Currently asy-drawfi) is realised as a wrapper around 
metapost-drawfi) but this may be changed. 

(methods specific for class cycle2D 7b) += (39b) <7e 8a > 

inline void asy-draw(ostream & ost, char * picture, 
const ex & xmin = -5, const ex & xmax = 5, 

const ex & ymin = -5, const ex & ymax = 5, const 1st & color = lst(), 

const char * more-options = "", 

bool with-header = true, int points-per-arc = 0) const 
{ metapost-draw(ost, xmin, xmax, ymin, ymax, color, more-options, with-header, 
points-per-arc, true, picture)', } 

inline void asy-draw(ostream & ost = stdwcout, 
const ex & xmin = -5, const ex & xmax = 5, 

const ex & ymin = -5, const ex & ymax = 5, const 1st & color = lst(), 

const char * more-options = "", 

bool with-header = true, int points-per-arc = 0) const 
{ metapost-draw{ost, xmin, xmax, ymin, ymax, color, more-options, with-header, 
points-per-arc, true); } 


Uses bool lie, ex 4a lib 39d 48a 48b 66a 66b 66c 67a, and 1st lie. 


Finally, we have a similar method which does not issue drawing command, instead it writes a definition for a (array 
of) path, which may be manipulated later. 

8a (methods specific for class cycle2D 7b) += (39b) <7f 

inline void asy-path(ostream & ost = stdr.cout, 

const ex & xmin = -5, const ex & xmax = 5, 
const ex & ymin = -5, const ex & ymax = 5, 
int points-per_arc = 0, bool is-continuation = false) const 
{ metapost-draw(ost : xmin , xmax , ymin , ymax , lst(), false, 
points-per_arc , true, true, is-continuation ); } 


Uses bool lie, ex 4a lib 39d 48a 48b 66a 66b 66c 67a, and 1st lie. 

2.7. An Example: Mobius Invariance of cycles. A quick illustration of the library usage is the symbolic calcu¬ 
lation which proves the Lem. 3.2 from [12]: We check that a Mobius transformation g £ SL 2 (&) acts on cycles by 
similarity g : C —> gCg~ x . We use the following predefined objects: 

cycle2D C(k,lst(l,n),m,e); 
ex W=lst(u,v); 

Firstly we define a cycle2D C2 by the condition between k, l and m in the generic cycle2D C that C passes 
through some point W. 

8b (Moebius transformation of cycles 8b) = (8e) 8d > 

C2 = C.subject-to(\st(C.passing(W ))); 


Uses 1st lie. 

The point gW is defined to be the Mobius transform of Why an arbitrary g. 

8c (Moebius transforms of W 8c) = (12a) 

const matrix gW=ex-to<vnsLtYvx>{clifford-moebius-map{sl2-clifford{a 1 b , c, d, e), W, e).subs(sl2-relationl, 
subs-options::algebraic \ subs-options::no-pattern).normal())\ 

Defines: 

matrix, used in chunks 2d, 4g, 5b, 16d, 18b, 25d, 38a, 40-44, 50-56, and 68. 

Finally we verify that the new cycle gCg~ x passes through P. This proves Lem. 3.2 from [10]. 

8d (Moebius transformation of cycles 8b) += (8e) <8b 

cout <C "Conjugation of a cycle comes through Moebius transformation: " 

<C C2.sl2similarity(a , b, c, d , es, S2).val(gW).subs(sl2-relationl , 

subs-options\:algebraic \ subs-options:'.no-pattern).normal().is-zero() 

<C endl <IC endl; 

Uses cycle 2a 2b 3e 3e 3e 3e 3e 46a 46b 47b 47b 48a 48a 48a 48a 48a. 

3. Demonstration through example 

We illustrate the library usage by the complete program which was used for computer-assisted proofs in the pa¬ 
per [10]. The numerous cross-references between these two papers are active hyperlinks. It is recommended to obtain 
PDF files for both of them from http: //arXiv. org and put into the same local directory. In this case clicking on a 
reference in a PDF reader will automatically transfer to the appropriate place (even in the other paper). 

3.1. Outline of the main(). The mainQ procedure does several things: 

(i) Makes symbolic calculations related to Mobius invariance; 

8e (List of symbolic calculations 8e)= (10b) 9a > 

(Moebius transformation of cycles 8b) 

(K-orbit invariance 12c) 

(Check Moebius transformations of zero cycles I2e) 

(Check transformations of zero cycles by conjugation 13c) 
cout <C end/; 


(ii) Calculates properties of orthogonality conditions and corresponding inversion in cycles; 

9a (List of symbolic calculations 8e)+= (10b) <8e 9b c> 

(Orthogonality conditions I3e) 

(Two points and orthogonality 14c) 

(One point and orthogonality I4e) 

(Orthogonal line 15b) 

(Inversion in cycle I5e) 

(Reflection in cycle 16b) 

(Yaglom inversion I6d) 
cout -C endl ; 

(iii) Calculates properties of s-orthogonality conditions and second type of inversion; 

9b (List of symbolic calculations 8e)+= (10b) <9a 9c o 

(Second orthogonality conditions 17a) 

(One point and s-orthogonality I7g) 

(s-Orthogonal line 18a) 

(s-Inversion in cycle 18c) 
cout -C endl ; 

(iv) Calculates various length formulae; 

9c (List of symbolic calculations 8e)+= (10b) <9b 9dc> 

(Distances from cycles I8g) 

(Lengths from centre 22b) 

(Lengths from focus 22d) 

(Infinitesimal cycle 23e) 
cout <C endl ; 

(v) Generates Asymptote output of the for illustrations. 

9d (List of symbolic calculations 8e)+= (10b) «9c 

Since we aiming into two targets simultaneously—validate our software and use it for mathematical proofs—there 
are many double checks and superfluous calculations. The positive aspect of this—a better illustration of the library 
usage. 

3.1.1. Program’s outline. Here is the main entry into the program and its outline. We start from some inclusions, 
note that GiNaC is included through <cycle./i>. 

9e (* 9e)= 9fo 

^include <cycle.h> 

^include <fstream> 

^define par_matr diag_matrix(lst (-1, 0)) 

^define hyp_matr diag_matrix(lst (-1, 1)) 

Defines: 

hypjnatr, used in chunk 35b. 
par_matr, used in chunks 33-35. 

Uses cycle 2a 2b 3e 3e 3e 3e 3e 46a 46b 47b 47b 48a 48a 48a 48a 48a and 1st lie. 

We try to make the output more readable both in simple text and modes. 

9f (* 9e)+= < 9e 10a > 

#define math_string << (output_latex?"$" : " ") 

#define wspaces (output J.atex?"\\quad ") 

Defines: 

math_string, used in chunks 12-15, 17-22, 24, and 25. 
wspaces, used in chunks 12-14 and 16-25. 


<9f 10b > 


The structure of the program is transparent. We declare all variables. 

10a (* 9e)+= 

(Declaration of variables lOd) 

(Subroutines definitions 12b) 
int main(){ 
cout <C boolalpha ; 

if {output-latex) cout <C latex ; 

Defines: 

main, never used. 

Then we make all symbolic calculations listed above. The exception catcher helps to identify the possible problems. 
10b (* 9e)+= < 10a 10c > 

try { 

(List of symbolic calculations 8e) 

} catch {exception &p) { 

cerr <C "***** Got probleml: " -C p.what{) -C endl ; 

> 


Defines: 

catch, used in chunks 42d and 50a. 

We end up with drawing illustration to our paper [10]. 

10c (* 9e)+= < 10b 

(Draw Asymptote pictures 26a) 

> 


3.1.2. Declaration of variables. First we declare all variables from the standard GiNaC classes here. 
lOd (Declaration of variables I0d)= (10a) I0ei> 

const char* eph-names=" eph." ■ 
const numeric half{ 1,2); 

const realsymbol a("a"), 6("b"), c("c"), d("d"), :r("x"), y{" y"), z(" z"), t{" t"), 
fc("k"), Z("L","1"), n("n"), // Cycles parameters 

fci("kl","\\tilde{k}"), D("ll","\\tilde{l>"), ml("ml" J "\\tilde{m>"), nl{" nl","\\tilde{n}"), 

m("u"), u("v"), ul{''u’ "), i >l{"v ’"), // Coordinates of points in R 2 
epsilon{ "eps", "Wepsilon"); // The "infinitesimal” number 

const varidx mt(symbol("nu", "\\nu"), 2), mtt(symbol("mu", "\\mu"), 2); 


Defines: 

eph_names, used in chunks 30a, 32a, 36d, and 37b. 

numeric, used in chunks 4-6, 11c, 18g, 20a, 21a, 30-32, 34a, 35c, 38a, 39c, 41-43, 46c, 47c, 49-51, and 56-68. 
realsymbol, never used. 

varidx, used in chunks 40-42, 44a, 50, 51b, 53a, and 56-58. 

Uses points 65b, u 62b, and v 62b. 

We need a plenty of symbols which will hold various parameters like ef, ef, s for the FSCc. 

I0e (Declaration of variables l0d)+= (10a) <l0d llai> 

const realsymbol sign{" s", "Wsigma"), signl{" si", "\\breve{\\sigma>"), //Signs of ef of ef 

sign2{" s2", "\\sigma_2"), sign3{" s3", "\\sigma_3"), signJ^{" s4", "Wmathring-fWsigma}"); 
int si , sil ; // Values of ef and ef for substitutions 

const matrix S2{ 2, 2, lst(l, 0, 0, jump-fnct{sign2 ))), 

55(2, 2, 1st(1, 0, 0, jump-fnct{sign3 ))), 

5/(2, 2, 1st(1, 0, 0, jump-fnct{sign4 ))); //Signs of l in the matrix representations of cycles 


Defines: 

matrix, used in chunks 2d, 4g, 5b, 16d, 18b, 25d, 38a, 40-44, 50-56, and 68. 
realsymbol, never used. 

si, used in chunks 16b, 20, 21, 26a, 30-32, 36, and 37. 
sil, used in chunks 20, 21, 26a, 30-32, and 36. 

Uses jump_fnct 38a and 1st lie. 


Here are several expressions which will keep results of calculations. 

(Declaration of variables lOd)+= (10a) «lOe lib > 

ex u2, v2, // Coordinates of the Moebius transform of {u, v) 
u3, v3, u4, v4, u5, v5, 

P, PI, // points on the plain 
K, LO, LI, /I Parameters of cycles 
Leri-C, // Expressions of Lengths 
P\ 


Uses ex 4a lib 39d 48a 48b 66a 66b 66c 67a, points 65b, u 62b, and v 62b. 

Next we define metrics (through Clifford units) for the space of points {M, e) and space of spheres {Ml, es). 
(Declaration of variables lOd)+= (10a) <lla llc> 

const ex M = diag-matrix{ lst(-l, sign)), // Metrics of point spaces 

e = clifford_unit(mu, M, 0), // Clifford algebra generators in the point space 

Ml = diag-matrix{ 1st(-1, signl)), // Metrics of cycles spaces 

es = clifford-unit{nu, Ml, 1); // Clifford algebra generators in the sphere space 


Defines: 

ex, used in chunks 2-8, 11a, 13, 14, 17, 19-25, 33b, 38-40, 42, 44-47, 49-60, 63, 65a, 67, and 68. 

Uses 1st lie. 

Now we define instances of cycle2D class. Some of them (like real-line or generic cycles C and Cl) are constants. 
(Declaration of variables i0d)+= ( 10 a) <iib lid > 

cycle2D C2, C3, C4, C5, C6, C7, C8, C9, CIO, Cll\ 

const cycle2D real-line{Q, lst(0, numeric(l)), 0, e), // the real line 

C{k, 1st {l, n), m, e), Cl{kl, 1st{11, nl), ml, e); // two generic cycles 
const cycle2D Zinf{ 0, lst(0, 0), 1, e), // the zero-radius cycle at infinity 

Z(lst(w, v), e), Zl{lst{u, v), e, 0, es), // two generic cycles of zero-radius 
Z2{ lst(u, v), e, 0, es, S2)\ 


Defines: 

cycle2D, used in chunks 6f, 7a, 12c, 14-18, 20a, 21g, 24a, 25b, 30a, 32, 33, 35-37, 39-41, 55-57, and 62-64. 

Uses cycle 2a 2b 3e 3e 3e 3e 3e 46a 46b 47b 47b 48a 48a 48a 48a 48a, 1st lie, numeric lOd, u 62b, and v 62b. 

For solution of various systems of linear equations we need the followings lsts. 

(Declaration of variables l0d)+= (10a) <a lie llec> 

1st eqns, eqnsl, 

vars=\st{kl, 11, ml, nl), 

solns, solnsl, // Solutions of linear systems 

signjval ; 


Uses 1st lie. 

Here are relational and lists of relational which will be used for automatic simplifications in calculations. They 
are based on properties of SL 2 {M) and values of the parameters. 

(Declaration of variables I0d)+s (10a) <lld llf > 

const relational sl2-relation = {c*b = a*d- 1), sl2-relationl = {a = (1 +b*c)-^d); // since ad — bc= 1 
const 1st signs-cube = 1st {pow{sign, 3) = sign, pow{signl, 3) = signl)', // s| = sJ, since Sj = —1, 0, 1 

const int debug = 0; 

const bool output-latex = true; 

Defines: 

bool, used in chunks 3-5, 7, 8a, 12d, 21f, 38a, 39a, 46c, 47a, 50a, 54-58, 62b, and 68. 
debug, used in chunks 14d, 15d, 17, 18, 20a, and 21f. 

1st, used in chunks 2b, 6-12, 14-26, 30, 32-37, 40c, 42a, 43a, 45, 46a, 50—58, 60-64, and 68b. 
relational, used in chunks 4e, 6, 45b, and 55c. 

Two generic points on the plain are defined as constant vectors (2 x 1 matrices). 

(Declaration of variables l0d)+= (10a) cite I2ac> 

const matrix W{2,1, lst(u, v)), Wl{2,1, 1st {ul, vl))', 


Defines: 

matrix, used in chunks 2d, 4g, 5b, 16d, 18b, 25d, 38a, 40-44, 50-56, and 68. 


We will also frequently use their Mobius transforms. 

12a (Declaration of variables lOd)+= (10a) <llf 22ac> 

const matrix gWl=ex-to<ma.trix.>(clifford-moebius-map(sl2-clifford(a, b , c, d, e), Wl, e) .subs(sl2-relationl, 
subs-optionsv.algebraic \ subs-options::no-pattern).normal ()); 

(Moebius transforms of W 8c) 


Defines: 

matrix, used in chunks 2d, 4g, 5b, 16d, 18b, 25d, 38a, 40—44, 50—56, and 68. 

We repeat some calculations several times for carious values of parameters, such calculations are gathered here as 
subroutines. 

12b (Subroutines definitions I2b)= (10a) 

(Parabolic Cayley transform of cycles 25b) 

(Check conformal property 20c) 

(Print perpendicular 2lg) 

(Focal length checks 22g) 

(Infinitesimal cycle calculations 23g) 


3.2. Mobius Transformation and Conjugation of Cycles. 

3.2.1. Transformations of K-orbits. As a simple check we verify that cycles given by the equation ( u 2 — av 2 ) — 
2v- —+ 1 = 0, see [10, Lem. 2.4] are IP-invariant, i.e. are AT-orbits. To this end we make a similarity of a cycle C2 
of this from with a matrix from K and check that the result coincides with C2. 

12c (K-orbit invariance 12c) = (8e) I2d > 

C2 = cycle2D(l, lst(0, (pow(t,-T)-sign*t)^r2), 1, e); 

cout <C "A K-orbit is preserved: " -C C2.sl2similarity(cos(x ), sin(x), -sin(x), cos(x), e).is-equal(C2) 

Uses cycle2D 6g 11c 11c 34b 35b 39d 39d 39d 39d 39d 48b 48b 48b 48b 48b 56b 56b 57d and 1st lie. 

We also check that C2 passing the point (0, t). 

I2d (K-orbit invariance 12c) += (8e) «l2c 

<C ", and passing (0, t) : " <C (bool) C2.passing(\st(0, f)) <C endl; 

Uses bool lie and 1st lie. 

3.2.2. Transformation of Zero-Radius Cycles. Firstly, we check some basic information about the zero-radius cycles. 
This mainly done to verify our library. 

I2e (Check Moebius transformations of zero cycles I2e)= (8e) I2f > 

cout < wspaces <C "Determinant of zero-radius Z1 cycle in metric e is : " 

mathstring <C canonicalize-clifford(Zl.det(e , S2)) mathstring <C endl 

-C wspaces <C "Focus of zero-radius cycle is: " mathstring <C Zl.focus(e) mathstring <C endl 
< wspaces <C "Centre of zero-radius cycle is: " math_string <C Zl.center(e) math_string <C endl 
-C wspaces <C "Focal length of zero-radius cycle is: " math_string <C Zl.focaLlengthf) math_string <C endl] 

Uses cycle 2a 2b 3e 3e 3e 3e 3e 46a 46b 47b 47b 48a 48a 48a 48a 48a, math_string 9f, and wspaces 9f. 

This chunk checks that Mobius transformation of a zero-radius cycle is a zero-radius cycle with centre obtained from 
the first one by the same Mobius transformation. 

I2f (Check Moebius transformations of zero cycles l2e)+= (8e) «l2e I2g > 

C2 = Zl.sl2similarity(a , b, c, d, e, S2)\ 

cout <C "Image of the zero-radius cycle under Moebius transform has zero radius: " 

■C canonicalize-clijford(C2.det(es , S2)).subs(sl2jrelationl 1 

subs-options::algebraic \ subs-options::no_pattern).is-zero() <!C endl; 

Uses cycle 2a 2b 3e 3e 3e 3e 3e 46a 46b 47b 47b 48a 48a 48a 48a 48a. 

Here we find parameters of the transformed zero-radius cycle C 2 = gZg ~ l . 

I2g (Check Moebius transformations of zero cycles l2e)+= (8e) <l2f 13a c> 

C2 = Z.sl2similarity(a 1 b, c, d , e, S2); 

K = C2.get_k(); 

L0 = C2. getJ(0); 

LI = C2.getJ(l ); 



Now we calculate the Mobius transformation of the centre of Z 

13a (Check Moebius transformations of zero cycles l2e)+= (8e) «l2g 13b c> 

u2 = gW.op( 0); 
v2 = gW.op( 1); 

And we finally check that gW coincides with the centre of the transformed cycle C2. This proves [10, Lem. 3.12]. 

13b (Check Moebius transformations of zero cycles I2e) += (8e) < 13a 

cout <C"The centre of the Moebius transformed zero-radius cycle is: " 

•C equality((u2*K-L0).subs(sl2-relation , subs-options\:algebraic \ subs-options::no-pattern )) <C ", " 

<C equality((v2*K-Ll).subs(sl2-relation, subs-options::algebraic \ subs-options::no-pattern )) 

<C endl ; 

Uses cycle 2a 2b 3e 3e 3e 3e 3e 46a 46b 47b 47b 48a 48a 48a 48a 48a and equality 38a 67d. 

3.2.3. Cycles conjugation. This chunk checks that transformation of a zero-radius cycle by conjugation with a cycle 
is a zero-radius cycle with centre obtained from the first one by the same transformation. 

Firstly we calculate parameters of C 2 = CZC . 

13c (Check transformations of zero cycles by conjugation 13c) = (8e) 13d c> 

C2 = Z.cyclesimilarity{C, e, S2, S3 ); 

cout <C "Image of the zero-radius cycle under cycle similarity has zero radius: " 

<C canonicalize-difford(C2.det(e, S2)).subs(sl2-relationl, 

subs-options::algebraic | subs-options::no-pattern).normal^).is-zeroij) <C endl, 

K = C2.get_k()', 

L0 = C2.getJ(0); 

LI = C2.getJ{l)-, 

Uses cycle 2a 2b 3e 3e 3e 3e 3e 46a 46b 47b 47b 48a 48a 48a 48a 48a. 

Then we check that it coincides with transformation point P which is calculated in agreement with above used matrices 
S2 and S3. This proves the result [10, Lem. 4.11] 

13d (Check transformations of zero cycles by conjugation 13c) += (8e) «l3c 

P= C.moebius-map{W , e, S2.mul(S3)); 
u2 = P.op( 0); 
v2 = P.op( 1); 

cout <C"The centre of the conjugated zero-radius cycle coinsides with Moebius tr: " 

<C equality(u2*K-L0) <C ", " <C equality(v2*K-Ll ) <C endl ; 

Uses cycle 2a 2b 3e 3e 3e 3e 3e 46a 46b 47b 47b 48a 48a 48a 48a 48a and equality 38a 67d. 

3.3. Orthogonality of Cycles. 

3.3.1. Various orthogonality conditions. We calculate orthogonality condition between two cycle2Ds by the identity 
mtr{C\C- 2 ) = 0. The expression are stored in variables, which will be used later in our calculations. 

Here is the orthogonality of two generic cycle2Ds... 

I3e (Orthogonality conditions I3e)= (9a) I3f > 

cout <C wspaces < C "The orthogonality is: " math_string 
<!C ( ex)C.is-orthogonal(Cl , es, S2) math_string <C endl 

Uses ex 4a lib 39d 48a 48b 66a 66b 66c 67a, math.string 9f, and wspaces 9f. 

... and then its reduction to orthogonality of two straight lines. 

I3f (Orthogonality conditions l3e)+= (9a) <l3e 14a i> 

<C wspaces <C "The orthogonality of two lines is: " math_string 
<C (ex)C.subs(k = 0).is-orthogonal(Cl.subs(kl = 0), es, S2) math_string <C endl, 


Uses ex 4a lib 39d 48a 48b 66a 66b 66c 67a, math_string 9f, and wspaces 9f. 


Here is the orthogonality of a generic cycle2D to a zero-radius cycle2D. This reduces to concurrence of the centre 
the zero-radius and generic cycle. 

14a (Orthogonality conditions l3e)+= (9a) <l3f 14b > 

cout <C wspaces <C "The orthogonality to z-r-cycle is: " math_string 
-C (ex)C.is-orthogonal(Z, es) math_string <C endl ; 

Uses cycle 2a 2b 3e 3e 3e 3e 3e 46a 46b 47b 47b 48a 48a 48a 48a 48a, ex 4a lib 39d 48a 48b 66a 66b 66c 67a, mat lust ring 9f, 
and wspaces 9f. 

Here is the orthogonality of two zero-radius cycle2Ds. 

14b (Orthogonality conditions l3e)+= (9a) «l4a 

C2 = cycle2D(lst(wl, vl), e, 0, S2 ); 

cout < wspaces <C "The orthogonality of two z-r-cycle is: " math_string 
-C ( ex)C2.is-orthogonal(Zl, es) math_string -C endl ; 

Uses cycle 2a 2b 3e 3e 3e 3e 3e 46a 46b 47b 47b 48a 48a 48a 48a 48a, cycle2D 6g 11c 11c 34b 35b 39d 39d 39d 39d 39d 48b 48b 48b 48b 
48b 56b 56b 57d, ex 4a lib 39d 48a 48b 66a 66b 66c 67a, 1st lie, math_string 9f, and wspaces 9f. 

This chunk hnds the parameters of a cycle C2 passing through two points (u,v), (ui,v\) and orthogonal to the given 
cycle C. This gives three linear equations with four variables which are consistent in a generic position. 

14c (Two points and orthogonality 14c) = (9a) I4dc> 

C2 = Cl.subject-to(lst(Cl.passing(W ), 

Cl.passing( Wl), 

Cl.is-orthogonal(C, es)), vars); 

Uses 1st lie. 

To find the singularity condition of the above solution we analyse the denominator of k , which calculated to be: 

— 2{u'(a\n + vk) — vl + {—kv' — <j\n)u + lv')n\ 

— u' 2 l + u' 2 uk + alv' 2 — u'u 2 k + u'v 2 ak + u'm — uakv' 2 + u 2 l — v 2 al — um 
I4d (Two points and orthogonality 14c) += (9a) <l4c 

if ( debug > 0) 

cout <C "Cycle through two point is possible and unique if denominator is not zero: " <C endl 
math_string <C C2.get_k() math_string <C endl <C endl] 

Uses debug lie and math.string 9f. 

3.3.2. Orthogonality and Inversion. Now we check that any orthogonal cycle comes through the inverse of any its 
point. To this end we calculate a generic cycle C2 passing through a point (u,v) and orthogonal to a cycle C. 

I4e (One point and orthogonality I4e)= (9a) I4gc> 

C2 = Cl.subject_to(lst(Cl.passing( W), 

Cl.is-orthogonal(C, es))); 

Uses 1st lie. 

Then we calculate another cycle C3 with an additional condition that it passing through the Mobius transform P 
of (u, v). 

I4g (One point and orthogonality l4e)+= (9a) <l4e I4hc> 

P = C.moebius-map(W., e, -Ml)] 

C3 = Cl.subject_to(lst(Cl.passing(P), 

Cl.passing(W), 

Cl.is-orthogonal(C, es)))] 

Uses 1st lie. 

Then we check twice in different ways the same mathematical statement: 

(i) that both cycles C2 and C3 are identical, i.e. the addition of inverse point does not put more restrictions; 
I4h (One point and orthogonality l4e)+= (9a) «l4g I5ai> 

cowt<C "Both orthogonal cycles (through one point and through its inverse) are the same: 1 

<SC C2.is-equal(C3) <C endl 



(ii) that cycle C2 passes through the inversion P as well. 

15a (One point and orthogonality l4e)+= (9a) <l4h 

-C "Orthogonal cycle passes through the transformed point: " 

<C C2.val)P).normal)).is-zero)) <C endl <C endl; 

Uses cycle 2a 2b 3e 3e 3e 3e 3e 46a 46b 47b 47b 48a 48a 48a 48a 48a. 

3.3.3. Orthogonal Lines. This chunk checks that the straight line C4 passing through a point (u,v) and its inverse P 
in the cycle C is orthogonal to the initial cycle C. 

(Orthogonal line I5b}= (9a) I5ci> 

C4 = Cl.subject-to)lst)Cl.passing)W), 

Cl.passing)P ), 

Cl.isJinear {))); 

ccmf<?C " Line through point and its inverse is orthogonal: " -C C4-inner-product)C, es).is-Zero)) -C endl] 

Uses 1st lie. 

We also calculate that all such lines intersect in a single point ( 113 , V3), which is independent from ( u,v ). This point 
will be understood as centre of the cycle C5 in § 3.3.4. 

(Orthogonal line 15b) += (9a) <l5b I5d > 

u3 = C.center)).op) 0); 
v3 = C4-roots(u3, false). op)0) .normal)); 

cout<€. "All lines come through the point " math_string <C" (" <C u3 <C " , " -C v3 <C")" mathstring <C endl] 
Uses math_string 9f. 

The double check is done next: we calculate the inverse PI of a vector ( u3+u , v3+v) and check that Pl-)u3 , v3) is 
collinear to ( u , v). 

(Orthogonal line 15b) += (9a) <l5c 

PI = C.moebius-map(lst(u3+u , v3+ v), e, -Ml); 
cout <C "Conjugated vector is parallel to (u,v) : " 

< (( Pl.op(0)-u3)*v-(Pl.op(l)-v3)*u).normal().is-zero () -C endl; 
if (debug >1) 

cout <C "Conjugated vector to (u, v) is: " math_string <C "(" <C ( Pl.op(0)-u3).normal () <C ", " 

< (Pl.op{\)-v3).normal)) <S ")" math_string <C endl; 

Uses debug lie, 1st lie, math_string 9f, u 62b, and v 62b. 

3.3.4. The Ghost Cycle. We build now the cycle C5 which defines inversion. We build it from two conditions: 

(i) C5 has its centre in the point ( u3 , v3 ) which is the intersection of all orthogonal lines (see § 3.3.3). 

(ii) The determinant of C5 with delta-sign is equal to determinant of C with signs defined by Ml. 

(Inversion in cycle I5e)= (9a) I5f > 

C5 = cycle2D(lst(it3, - v3*jump-fnct(sign )), e, C.det{e , Ml)).subs{signs-cube , 
subS-options::algebraic \ subs-options\:no-pattern); 

Uses cycle2D 6g 11c 11c 34b 35b 39d 39d 39d 39d 39d 48b 48b 48b 48b 48b 56b 56b 57d, jump_fnct 38a, and 1st lie. 

As a consequence we find out that C5 has the same roots as C. 

(Inversion in cycle l5e)+= (9a) «l5e I5g > 

cout <C "C5 has common roots with C : " <C ( C5.val(lst)C.roots().op(0 ), 0 )).normal)).is-zero)) 

A C5.val)lst)C.roots)).op)l), 0)). normal)).is-zero))) -C endl 

<C "$\\chi(\\sigma)$-centre of C5 is equal to $\\breve{\\sigma}$-centre of C: " 

<C )C5.center)diag-matrix)lst)-l,jump-fnct)sign))), true)- C.center)es, true)) .normal)).is-zero)) <!C endl; 

Uses jump_fnct 38a and 1st lie. 

Finally we calculate point PI which is the inverse of (^ 3 ,^ 3 ) in C5. 

(Inversion in cycle l5e)+= (9a) <jl5f 16a c> 

PI = C5.moebius-map)W, e, diag-matrix)lst)l, -jump-fnct)sign)))); 


Uses jump_fnct 38a and 1 st lie. 


The final check: P (inversion in C5 in terms of sign) coincides with P —the inversion in C in terms of signl , see 
chunk 14f. 

16a (Inversion in cycle l5e)+= (9a) <i5g 

cout <C "Inversion in (C5, sign) coincides with inversion in (C, signl): " 

-C ( Pl-P).subs(signs-Cube , subs-options::algebraic \ subs-options::no-pattern) .normalf) .is-zerof) 

<C endl; 

3.3.5. The real line and reflection in cycles. We check that conjugation CiRCi maps the real-line to the cycle C and 
wise verse for the properly chosen Cl, see [10, Lem. 4.16]. 

16b (Reflection in cycle 16b) = (9a) 16c c> 

for (s*=-l; si<2; si+= 2) { 

C9 = cycle2D(fc* signl, 1st {hsignl, msignl 

+si*sqrt(- C. det( es, (new tens delta) ~^setflag(status-flags:: dynallocated), k) * signl)), m * signl, es ); 
cout <C "Inversion to the real line (with " <C (si=-l? : "+") <C " sign): " <C endl 

< wspaces <C "Conjugation of the real line is the cycle C: " 

<C real-line.cyclesimilarity(C9, es).is-equal(C) <C endl 

< wspaces <C "Conjugation of the cycle C is the real line: " 

<C C.cyclesimilarity(C9, es) .is-equal(reaLline) <C endl 

Uses cycle 2a 2b 3e 3e 3e 3e 3e 46a 46b 47b 47b 48a 48a 48a 48a 48a, cycle2D 6g 11c 11c 34b 35b 39d 39d 39d 39d 39d 48b 48b 48b 48b 
48b 56b 56b 57d, 1st lie, si lOe, and wspaces 9f. 

We also check two additional properties which caracterises the inversion cycle C9 in term of common roots of C 
[10, Lem. 4.16.(ii)] and C passing through C9 centre [10, Lem. 4.16.(iii)] . 

16c (Reflection in cycle 16b) += (9a) «l6b 

<C wspaces <C "Inversion cycle has common roots with C: " 

<C ( C9.val(lst(C.roots().op(0), 0)) .numerf).normalf).is-zerof) 

A C9.val(lst(C.roots().op(l), O)).numer().normal().is-Zero()) -C endl 

< wspaces <C "C passing the centre of inversion cycle: " 

cycle2D(C, es).val(C9.centerQ).numerf).subs(signl= sign, subs-options::no-pattern).normalf) 

.subs(pow(sign, 2)=1, subs-options::algebraic \ subs-options::no-pattern).is-zero() <C endl; 

> 

Uses cycle 2a 2b 3e 3e 3e 3e 3e 46a 46b 47b 47b 48a 48a 48a 48a 48a, cycle2D 6g 11c 11c 34b 35b 39d 39d 39d 39d 39d 48b 48b 48b 48b 
48b 56b 56b 57d, 1st lie, and wspaces 9f. 

3.3.6. Yaglom inversion of the second kind. In the book [15, § 10] the inversion of second kind related to a parabola 
v = k(u — l) 2 + m is defined by the map: 

(u, v) i—> (u, 2 (k(u — l) 2 + m) — v). 

We shows here that this is a composition of three inversions in two parabolas and the real line, see [12, Prop. 4.17]. 
I6d (Yaglom inversion I6d)= (9a) 

cout <C "Yaglom inversion of the second kind is three reflections in the cycles: " 

<C ( reaLline.moebius-map{cyc\e2T)(\st(l, 0), e, -m-^-k).moebius-map(cyc\e2D(\st(l, 2 *m), e, -m-h-k) 

,moebius-map( W))).subs(sign=0) 

-matrix(2,l,lst(if, 2*(k*pow(u-l,2)+m)-v))).normal().is-zero() <C endl; 


Uses cycle2D 6g 11c 11c 34b 35b 39d 39d 39d 39d 39d 48b 48b 48b 48b 48b 56b 56b 57d, 1st lie, matrix 8c lOe 1 if 12a 25a, u 62b, 
and v 62b. 


3.4. Orthogonality of the Second Kind. We study now the orthogonality condition of the second kind (s- 
orthogonality), [10, § 4.3]. 

3.4.1. Expressions for s-orthogonality. One more simple consistency check: the real-line is invariant under all Mdbius 
transformations. 

17a (Second orthogonality conditions 17a) = (9b) 17b > 

cout <C "The real line is Moebius invariant: " -C reaLline.is-equal(reaLline.sl2similarity(a, b, c, d, es)) <C endl 

-C "Reflection in the real line: " 

mathstring <C Z. cycle-similarity (real-line, es) .normalize)) mathstring <C endl 
<C "Reflection of the real line in cycle C: " <C endl 
mathstring <C real-line.cycle-similarity(C, es, S2, S3) mathstring <C endl ; 

Uses cycle 2a 2b 3e 3e 3e 3e 3e 46a 46b 47b 47b 48a 48a 48a 48a 48a and math.string 9f. 

The orthogonality condition of the second kind between two different cycles is calculated by the identity [10, § 4.3] 

mr (CiC 2 Ci,R) =0. 

Here is s-orthogonality of two generic cycle2Ds... 

17b (Second orthogonality conditions 17a) += (9b) <i7a 17c > 

cout < wspaces <C "The s-orthogonality is : " mathstring 
<C ( ex)C.iss-orthogonal(Cl, es, S2) mathstring <C endl 

Uses ex 4a lib 39d 48a 48b 66a 66b 66c 67a, math.string 9f, and wspaces 9f. 

... and its reduction to the straight lines case. 

17c (Second orthogonality conditions 17a) += (9b) <l7b I7dc> 

< wspaces -C "The s-orthogonality of two lines is: " mathstring 
-C (ex)C.subs(k = 0).iss-orthogonal(Cl.subs(kl=0), es, S2) mathstring <C endl, 

Uses ex 4a lib 39d 48a 48b 66a 66b 66c 67a, math.string 9f, and wspaces 9f. 

Here is s-orthogonality of a generic cycle2D to a zero-radius cycle2D. 

I7d (Second orthogonality conditions 17a) += (9b) <l7c I7ec> 

cout <C wspaces <C "The s-orthogonality to z-r-cycle is first way: " -C endl 
mathstring <C (ex)C.iss-orthogonal(Zl, es, S2) mathstring <C endl, 

Uses cycle 2a 2b 3e 3e 3e 3e 3e 46a 46b 47b 47b 48a 48a 48a 48a 48a, ex 4a lib 39d 48a 48b 66a 66b 66c 67a, math_string 9f, 
and wspaces 9f. 

Since s-orthogonality is not symmetric [10, § 4.3], we calculate separately s-orthogonality of a zero-radius cycle2D to 
a generic cycle2D. 

I7e (Second orthogonality conditions 17a) += (9b) <i7d l7fo 

cout <C wspaces <C "The s-orthogonality to z-r-cycle is second way: " <C endl 
mathstring <C (ex)Zl.iss-orthogonal(C, es, S2) mathstring <C endl, 

Uses cycle 2a 2b 3e 3e 3e 3e 3e 46a 46b 47b 47b 48a 48a 48a 48a 48a, ex 4a lib 39d 48a 48b 66a 66b 66c 67a, math_string 9f, 
and wspaces 9f. 

Here is s-orthogonality of two zero-radius cycle2Ds. 

I7f (Second orthogonality conditions 17a) += (9b) <i7e 

C9 = cycle2D(lst(Ml, vl), e); 

cout <C wspaces <C "The s-orthogonality of two z-r-cycle is: " <C endl 
mathstring <C (ex)Zl.iss-orthogonal(C9, es, S2) mathstring <C endl, 

Uses cycle 2a 2b 3e 3e 3e 3e 3e 46a 46b 47b 47b 48a 48a 48a 48a 48a, cycle2D 6g 11c 11c 34b 35b 39d 39d 39d 39d 39d 48b 48b 48b 48b 
48b 56b 56b 57d, ex 4a lib 39d 48a 48b 66a 66b 66c 67a, 1st lie, math_string 9f, and wspaces 9f. 

3.4.2. Properies of s-orthogonality. Find the parameters of cycle passing through a point and s-orthogonal to the given 
one 

I7g (One point and s-orthogonality I7g)= (9b) 

C6 = Cl.subject-to(lst(Cl.passing(W), C.iss-orthogonal(Cl, es, S2)))- 
if ( debug > 1) 

cout <C "Cycle s-orthogonal to (k, (1, n), m) is: " <C endl 
mathstring <C C6 mathstring <C endl; 


Uses debug lie, 1st lie, and math_string 9f. 


18a 


(9b) 18b > 


Check the orthogonality of the line through a point to the cycle. 

(s-Orthogonal line I8a)= 

Cl = C6.subject-to(lst(C6.isJinear())); 
u4 = C.center().op( 0); 
v4 = C7.roots(u4, false). op(0) .normalf)', 

Uses 1st lie. 

All orthogonal lines come through the same point, which the focus of the cycle Cwith respect to metric (-1, -signl). 
18b (s-Orthogonal line 18a) += (9b) <l8a 

cout <C "All lines come through the focus related $\\breve{e}$ : " 

<C {C.focus{diag-matrix{ 1st (-1, -signl)), true)-matrix(2, 1, lst(u4, v4)))-normal{) .is-zerof) <C endl; 

Uses 1st lie and matrix 8c lOe Ilf 12a 25a. 

3.4.3. Inversion from the s-orthogonality. We express s-orthogonality to a cycle C through the usual orthogonality to 
another cycle C8. This cycle is the reflection of the real line in C, see 3.3.5. 

18c (s-Inversion in cycle 18c) = (9b) I8dc> 

C8= real-line. cyclesimilarity(C, es, diag-matrix(lst(l, signl)), diag-matrix(lst(l, jump-fnct(sign)))).normalize(n*k); 
if ( debug > 1) 

cout -C "C8 is : " mathstring -C C8 mathstring <C endl ; 

Uses debug lie, jump_fnct 38a, 1st lie, and math_string 9f. 

We check that C8 has common roots with C. 

I8d (s-Inversion in cycle 18c) += (9b) «l8c I8ec> 

cout <C "C8 has common roots with C : " <C ( C8.val(lst(C.roots().op(0 ), O)).numer().normal().is-zero() 

A C8.val(lst(C.roots().op(l), 0)) .numer().normalf).is-zero()) <C endl, 

Uses 1st lie. 

This chunk checks that centre of C8 coincides with focus of C. 

I8e (s-Inversion in cycle 18c) += (9b) <l8d I8f > 

cout -C "$\\chi(Wsigma)$-center of C8 coincides with $\\breve{\\sigma}$-focus of C : " 

< (C8.center(diag-matrix(lst(-l,jump-fnct(sign))), true) 

- C.focus(diag-matrix( 1st(-1 , -signl)), true)) . evalmQ.normal (). is-zero-matrixlf) 

<C endl, 

Uses jump_fnct 38a and 1st lie. 

Finally we check that s-inversion in C defined through s-orthogonality coincides with inversion in C8. 

I8f (s-Inversion in cycle 18c) += (9b) «l8e 

PI = C8.moebius-map(W, e, diag-matrix{ lst(l, -jump-fnct(sign)))).subs(signs-cube, subs-options::algebraic | 
subs-options: :no-pattern).normal (); 

cout <C "s-Inversion in C coincides with inversion in C8 : " 

<C C6.val(Pl).normal().subs(signs-Cube, subs-options::algebraic | subs-options::no-pattern).normal().is-zero() 

<C endl ; 

Uses jump_fnct 38a and 1st lie. 

3.5. Distances and Lengths. 

3.5.1. Distances between points. We calculate several distances from the cycles. 

The distance is given by the extremal value of diameters for all possible cycles passing through the both points [12, 
Defn. 5.2]. Thus we first construct a generic cycle2d CIO passing through two points (u,v) and (u',v r ). 

I8g (Distances from cycles I8g)= (9c) 19a > 

CIO = cycle2D(numeric(l), 1st (Z, n), m, e); 

CIO = C10.subject-to(lst(C10.passing(W), ClO.passing(Wl)), lst(m, n, l )); 
if (debug > 0) cout <C wspaces <C "CIO is: " <C CIO <§C endl, 


Uses cycle2D 6g 11c 11c 34b 35b 39d 39d 39d 39d 39d 48b 48b 48b 48b 48b 56b 56b 57d, debug lie, 1st lie, numeric lOd, and wspaces 9f. 


Then we calculate the square of its radius as the value of the determinant D. The point l of extremum Len-C is 
calculated from the condition D[ = 0. 

19a (Distances from cycles l8g)+= (9c) <i8g 19b > 

ex D = 4*C10.det(es ); 

Len-C = D.subs(lsolve(\st(D.diff(l) = 0), 1st (l))). normalQ; 


Uses ex 4a lib 39d 48a 48b 66a 66b 66c 67a and 1st lie. 
Now we check that Len-C is equal to [10, Lem. 5.5] 


d 2 {y,y') = 


cr((u — u') 2 — a(v — v') 2 ) + 4(1 — cra)i 


-((u — u') 2 — a( v — v 1 ) 2 ), 


(u — u') 2 d — (v — v') 2 

19b (Distances from cycles l8g)+= (9c) <l9a 19c c> 

cout <C "Distance between (u,v) and (u^v’) in elliptic and hyperbolic spaces is " <C endl ; 


if (outputJatex) { 

ex dist = ( signl*(pow(u-ul,2)~sign*pow(v-vl : 2))+4:*(l-sign*signl)*v*vl)*(pow(u-ul,2 ) 

- sign* pow(v-v 1,2)) ~(pow(u-u 1,2)* sign l-pow(v-vl,2))\ 
cout <C "\\(Wdisplaystyle " <C dist <C "\\) : " <C (Leri-C-dist).normalQ.is-zero() <C endl; 
} else 
cout < endl 

< " s1*((u — u ; )~2-s*(v-v ; )~2)+4*(l-s*sl)*v*v’)*((u-u’)~2-s*(v-v ; )~2)" 

<C endl 

« " - : " 

<C (Len-C-(signl*(pow(u-ul,2)-sign*pow(v-vl,2))+4:*(l-sign*signl)*v*vl)*(pow(u-ul : 2) 
-sign*pow(v-vl,2))-^-(pow(u-ul,2)*signl-pow(v-vl,2))).normal().is-zero() <C endl 
(u-u’) ~2*sl-(v-v J ) ~2" endl <C endl; 


Uses ex 4a lib 39d 48a 48b 66a 66b 66c 67a, u 62b, and v 62b. 

Conformity is verified in the same chunk (see § 3.5.2) for this and all subsequent distances and lengths. Value si = -1 
initiates conformality checks only in elliptic and hyperbolic point spaces. 

19c (Distances from cycles l8g)+= (9c) <l9b I9d > 

check-conformality(Len-C, -1); 

Cll = C10.subs(lsolve(lst(D.diff(l) = 0), lst(/))); 
print-perpendicular( Cll); 

Uses check.conf ormality 20c, 1st lie, and print_perpendicular 21g. 

In parabolic space the extremal value is attained in the point i(it + itl), since it separates upward-branched parabolas 
from down-branched. 

I9d (Distances from cycles l8g)+= (9c) <l9c l9eo 

Len-C = D.sub s(lst (sign =0, l = (u+ul)*half)).normalQ; 

cout <C "Value at the middle point (parabolic point space):" <C endl <C wspaces 
mathstring <C Len-C mathstring <C endl; 

Uses 1st lie, math.string 9f, u 62b, and wspaces 9f. 

Value si = 0 initiates conformality checks only in the parabolic point space. 

I9e (Distances from cycles l8g)+= (9c) <l9d 20a > 

check-Conformality(Len-C , 0); 

Cll = C10.subs(\st(sign =0, 1= (u+ul)*half)); 
print-perpendicular( Cll); 


Uses check.conformality 20c, 1st lie, print_perpendicular 21g, and u 62b. 




We need to check the case v = v' separately, since it is not covered by the above chunk. This is done almost identically 
to the previous case, with replacement of l by n, since the value of l is now fixed. 

20a (Distances from cycles i8g)+= (9c) «i9e 20bi> 

CIO = cycle2D(numeric(l), 1st (Z, n), m , e); 

CIO = C10.subject-to(lst(C10.passing(W), 

C10.passing(lst(ul, v)))); 
if (debug >1) 

cout <€. wspaces -C CIO -C endl; 

Uses cycle2D 6g 11c 11c 34b 35b 39d 39d 39d 39d 39d 48b 48b 48b 48b 48b 56b 56b 57d, debug lie, 1st lie, numeric lOd, v 62b, 
and wspaces 9f. 

This time the extremal point n is found from the condition D' n = 0. 

20 b (Distances from cycles I8g) += (9c) «20 a 

D = 4 *C10.det(es); 

Leri-C = D.subs(lsolve(\st(D.diff(n ) = 0), lst(n))).normaZ(); 
cout <C "Distance between (u,v) and (u’,v): " -C endl 

<C wspaces -C "Value at critical point: " <C endl <C wspaces math_string <C Len-C math_string 
-C endl <C endl; 

Uses 1st lie, math_string 9f, u 62b, v 62b, and wspaces 9f. 

3.5.2. Check of the conformal property. We check conformal property of all distances and lengths. This is most time- 
consuming portion of the program and it took few minutes on my computer. The rest is calculated within twenty 
seconds. 

20 c (Check conformal property 20 c)= (12b) 20dc> 

void check-conformality(const ex & Ten_c, int si = 3) { 

(Evaluate the fraction 2ld) 

Defines: 

check_conf ormality, used in chunks 19, 22c, and 23b. 

Uses ex 4a lib 39d 48a 48b 66a 66b 66c 67a and si lOe. 

Several times we fork for two cases: the first one if the check is done for all signs combinations simultaneously. 

20 d (Check conformal property 20 c) += (12b) <20c 20ec> 

if (si > 2) 

cout <C wspaces <C "This distance/length is conformal:" ; 

Uses si lOe and wspaces 9f. 

The second case is we output coresponding results for different metric signs. 

20 e (Check conformal property 20 c) += (12b) «20d 20fc> 

else 

cout -C wspaces <C "Conformity in a cycle space with metric: E P H " <C endl; 

Uses cycle 2a 2b 3e 3e 3e 3e 3e 46a 46b 47b 47b 48a 48a 48a 48a 48a and wspaces 9f. 

However we make the substitution of all possible combinations of sign and signl (an initial value of si should be set 
before in order to separate parabolic case from others). The first loop is for point space metric sign. 

20 f (Check conformal property 20 c) += (12b) <20e 20gc> 

do { 
if (si > 1) 
sil = 2; 
else { 

cout <C wspaces <C "Point space is " <C eph_case(si ) <C ": 
sil = -1; 

> 

Uses eph_case 38a 68a, si lOe, sil lOe, and wspaces 9f. 

The second loop is for cycle space metric sign. 

20 g (Check conformal property 20 c) += (12b) <20f 2lac> 

do { 

if (si < 2) 


Uses si lOe. 


21a 


However the substition of signs is not done for dummy loops. 

(Check conformal property 20 c) += (12b) <20g 21b > 

Len-cD = Len-fD.sub s(lst (sign = numeric(sz), signl = numerical)), 
subs-options:: algebraic \ subs-options::no-pattern).normal()', 

Uses 1st lie, numeric lOd, si lOe, and sil lOe. 

But even for dummy loops we make a check the conformity. 

2lb (Check conformal property 20 c) += (12b) <2la 2lcc> 

(Find the limit 2 le) 

(Check independence 2lf) 

and then finalise all loops. 

2lc (Check conformal property 20 c) += (12b) <2lb 

sil++- 

} while (sil < 2); 
cout -C endl ; 
si+=2; 

} while (si < 2); 

> 

Uses si lOe and sil lOe. 

To this end we consider the ratio of distances between (u,v) and ( u + tx,v + ty ) and between their images gW and 
gWl under the generic Mobius transform. 

2ld (Evaluate the fraction 2 ld)= (20c) 

ex Len-cD= ((Len-C.subs(\st(u = gW.op(0), v^gW.op(l), ul = gWl.op( 0), 
vl=gWl.op( 1)), subs-options:: algebraic \ subs-options::no-pattern) 

^-Len-c).subs(lst(ul=u+t*x, vl=v+t*y ), subs-options::algebraic \ subs-options::no-pattern )); 
ex Len_fD = Len_cD; 

Uses ex 4a lib 39d 48a 48b 66a 66b 66c 67a, 1st lie, u 62b, and v 62b. 

If Leri-cD has the variable t , we take the limit t —► 0 using the power series expansions. 

2le (Find the limit 2 le)= (2lb) 

if ( Len_cD.has(t )) 

Len-cD = Len-cD. series(t^0,1) .op(0) .normalQ; 

The limit of this ratio for t, —> 0 should be independent from (x,y) (see [10, Defn. 5.9]). 

2lf (Check independence 2lf)= (2lb) 

bool is-Conformal = -> (Leri-cD.is-zero () V Leri-cD.has(t) 

V Leri-cD.has(x) V Len-cD.has(y))', 
cout <C " " <C is-conformal ; 
if ( debug > 0 V (-> is-conformal A (si > 2))) 

cout <C The factor is: " <C endl <C wspaces math_string <C Len-cD.normal () math_string ; 

Uses bool lie, debug lie, math_string 9f, si lOe, and wspaces 9f. 

3.5.3. Calculation of Perpendiculars. Lengths define corresponding perpendicular conditions in terms of shortest 
routes, see [10, Dcfn. 5.13]. 

2lg (Print perpendicular 2lg)= (12b) 

void print_perpendicular(const cycle2D & C) { 
cout <C wspaces <C "Perpendicular to ((u,v); (u^v 1 )) is: " 
math_string <C (C.getJ(l)+sign*C.get_k()*vl).normal() mathstring <C " 
mathstring <C (C.get-l(0)-C.get-k()*ul).normal() mathstring <C endl <C endl ; 

> 


Defines: 

print_perpendicular, used in chunks 19 and 22c. 

Uses cycle2D 6g 11c 11c 34b 35b 39d 39d 39d 39d 39d 48b 48b 48b 48b 48b 56b 56b 57d, math.string 9f, u 62b, v 62b, and wspaces 9f. 


3.5.4. Length of intervals from centre. We calculate the lengths derived from the cycle with a centre at one point and 
passing through the second, see [ 0, Defn. 5.3]. 

Firstly we need some more imaginary units, to accommodate different types of centres (foci). 

22a (Declaration of variables l0d}+= (10a) «l2a 23dc> 

ex sign5=signf\ 

ex e4 = clifford-unit(nu , diag-matrix{ lst(-l, signf )), 2); 

Uses ex 4a lib 39d 48a 48b 66a 66b 66c 67a and 1st lie. 

Then we build a cycle2D Cll which passes through (u',v') and has its centre at (u,v). 

22b (Lengths from centre 22 b) = (9c) 22c > 

Cll = C.subject-to(lst{C.passing(Wl), C.is-normalized {))); 

Cll = Cll.subject-to(lst(Cll.center().op(0) = u, Cll.center(ef).op( 1)= v))\ 

Uses 1st lie, u 62b, and v 62b. 

Then the distance is radius the Cll , see [10, Lem. 5.7. (i)]. We check conformity and calculate the perpendicular at 
the end. 

22c (Lengths from centre 22 b) += (9c) <22b 

Len-C = Cl 1. det(es) .normalf); 

cout <C "Length from *center* between (u,v) and (u’,v’):" <C endl <C wspaces 
mathstring <C Len-C mathstring <C endl ; 
check-conformality{Len-c ); 
print-perpendicular( Cll)', 

Uses check.conf ormality 20c, math.string 9f, print .perpendicular 21g, u 62b, v 62b, and wspaces 9f. 

3.5.5. Length of intervals from focus. We calculate the length derived from the cycle with a focus at one point. To 
use the linear solver in GiNaC we need to replace the condition C10.focus().op(l) = v by hand-made value for the 
parameter n. 

There are two suitable values of n which correspond upward and downward parabolas, which are expressed by plus 
or minus before the square root. After the value of length was found we master a simpler expression for it which 
utilises the focal length p of the parabola. 

22d (Lengths from focus 22 d)= (9c) 22ec> 

focaLlength_check(sign5*(-(vl-v)+sqrt(sign5*pow((ul-u ), 2 )+pow((vl-v), 2) -sign5*sign*pow(vl, 2)))); 

Uses f ocal.length.check 22g, u 62b, and v 62b. 

This chunk is similar to an above one but checks the second parabola (the minus sign before the square root). 

22e (Lengths from focus 22 d)+= (9c) <22d 22f > 

focaLlength_check(sign5*(-(vl-v)-sqrt(sign5*pow((ul-u), 2 )+pow((vl-v), 2) -sign5*sign*pow(vl, 2)))); 

Uses f ocal.length.check 22g, u 62b, and v 62b. 

We need to verify separately the case of sign5=0 , in this case p has a rational value. 

22f (Lengths from focus 22 d)+= (9c) <22e 

sign5= 0; 

focaLlength_check((pow(ul-u,2)-sign*pow(vl,2))-^-(vl-v)-i-2)', 

Uses f ocal.length.check 22g, u 62b, and v 62b. 

Again to avoid non-linearity of equation, we first construct a desired cycle. 

22g (Focal length checks 22 g)= (12b) 23a > 

void focalJength_check(const ex & p) { 

cout <C "Length from *focus* check for " math_string <C "p = " p mathstring <C endl ; 

Cll = C.subject-to(\st(C.passing(Wl ), &= 1, 1= u, n= p)); 

Defines: 

f ocal.length.check, used in chunk 22. 

Uses ex 4a lib 39d 48a 48b 66a 66b 66c 67a, 1st lie, math_string 9f, and u 62b. 


And now we verify that the length is equal to (1 — cri)p 2 — 2 vp, see [10, Lem. 5.7.(ii)] - 
23a (Focal length checks 22 g)+= (12b) <22g 23b > 

Len-C = Cll.det(es).subs(pow(sign4,2)=l,subs-options::algebraic \ subs-options::no-patterri).normal (); 
cout <C wspaces <C "Length between (u,v) and (u 1 , v’) is equal to " 

<C (outputJatexl "\\( (Wmathring(Wsigma)-\\breve{\\sigma})p~2-2vp\\) : "(s4-sl)*p~2-2vp: ") 

<C ( Len-C - (( sign5-signl)*pow(p , 2) - 2 *v*p)) .subs(signs-cube, subs-options::algebraic \ subs-options::no-pattern).expand() 
.subs(pow(sign4,2)=l,subs-options::algebraic \ subs-options::no-pattern).normal().is-zero () <C endl\ 

Uses u 62b, v 62b, and wspaces 9f. 

and we check all requested properties for Cll: it passes ( ul , vl) and has focus at ( u , v). 

23b (Focal length checks 22 g)+S (12b) «23a 23c i> 

cout <C wspaces <C "checks: Cll passes through (u’ , v’): " <C Cll.val(Wl).normal().is-zero() 

<C Cll focus is at (u, v) : " <C ( Cll.focus().op(0).normal()-u).is-zero( ) <C endl, 
check-Conformality(Leri-c ); 


Uses check.conf ormality 20c, u 62b, v 62b, and wspaces 9f. 

We finally verify that focal perpendiculars are multiples of the vector (av'+p, u—u'), see [10, E-it:focal-perpendicularity]. 
23c (Focal length checks 22g)+= (12b) <23b 

cout <C wspaces <C "Perpendicular to ((u,v); (u’.v 1 )) is " 

<C (outputJatex 7 "\\((\\sigma v’+p, u-u’)\\): " : "(s*v ; +p, u-u 1 ): ") 

■C (( Cll.getJ(l)+sign*Cll.get_k()*vl-(sign*vl+p)).normal().is-zero () 

A ( Cll.getJ(0)-Cll.get-k()*ul-(u-ul)).normal().is-zero ()) 

■C endl endl ; 

} 


Uses u 62b, v 62b, and wspaces 9f. 


3.6. Infinitesimal Cycles. The final bit of our calculation is related with the infinitesimal radius cycles, see [10, § 6.1]. 
Some additional parameters. 

23d (Declaration of variables l0d}+= (10a) <22a 25ac> 

possymbol vp("vp","v_p"); //the positive instance of symbol v 
ex displ ; //displacement of the focus 


Uses ex 4a lib 39d 48a 48b 66a 66b 66c 67a and v 62b. 

3.6.1. Basic properties of infinitesimal cycles. @We define an infinitesimal cycle CIO such that its squared radius (det) 
is an infinitesimal number e 2 and focus is at ( u , v ). This defined by the cycle (1, uo, n , u^ + 2nvQ — an 2 ) where n satisfies 
to the equation 

(3.1) (<7 — d)n 2 — 2 von + e 2 = 0. 

Only one root of the quadratic case produces a cycle with an infinitesimal focal length, and we consider it here: 

23e (Infinitesimal cycle 23e)= (9c) 23fc> 

inftnitesimaLcalculations(n=(vp-sqrt(pow(vp,2)-pow(epsilon,2)*(sign4-signl)))-^-(sign4-signl))-, 

Uses infinitesimal.calculations 23g. 

2 

The second expression for an infinitesimal cycle for the case a = a is given by the substitution n = — which the 
root of (3.1) in this case. 

23f (Infinitesimal cycle 23e)+= (9c) c23e 

infinitesimaLcalculations(\st ( n=pow( epsilon ,2)-P2 ^rvp, sign4= signl )); 

Uses infinitesimal.calculations 23g and 1st lie. 

We organise the infinitesimal cycles check as a separate subroutine and start it from several local variables definition. 
23g (Infinitesimal cycle calculations 23g)= (12b) 24a i> 

void infinitesimaLcalculations( const ex & nval) { 
exmap smap ; 
smap[v]=vp; 

Defines: 

inf initesimal.calculations, used in chunk 23. 

Uses ex 4a lib 39d 48a 48b 66a 66b 66c 67a and v 62b. 


24a (Infinitesimal cycle calculations 23g)+= (12b) «23g 24b > 

CIO = cycle2D(l, lst(it, n), pow(u,2)+2*n*vp-pow(n,2)*sign4-, e).subs(nval); 
cout <C "Inf cycle is: " mathstring <C CIO mathstring <C endl; 
cout <C "Square of radius of the infinitesimal cycle is: " 
mathstring <C C10.det(es). sub s(signs-cube, subs-options::algebraic 
| subs-options::no-pattern).normal() mathstring endl, 

Uses cycle 2a 2b 3e 3e 3e 3e 3e 46a 46b 47b 47b 48a 48a 48a 48a 48a, cycle2D 6g 11c 11c 34b 35b 39d 39d 39d 39d 39d 48b 48b 48b 48b 
48b 56b 56b 57d, 1st lie, math_string 9f, and u 62b. 

Then we verify that in parabolic space it focus is in the point (w, v) and the focal length is an infinitesimal. 

24b (Infinitesimal cycle calculations 23g)+= (12b) «24a 24c c> 

cout -C "Focus of infinitesimal cycle is: " mathstring <C C10.focus{eJf).subs{nval) mathstring <C endl 
•C "Focal length is: " mathstring -C C10.focaLlength().series(epsilon=0,3).normalQ mathstring <C endl; 

cout <C "Infinitesimal cycle passing points" mathstring -C "(u+" <C epsilomx <C" , vp+" 

<C lsolve(C10.subs(sign=0).passing(lst(u+epsilon*x,vp+y)),y).series(epsilon^0,3).normal() 

•C "), " mathstring <C endl, 

Uses cycle 2a 2b 3e 3e 3e 3e 3e 46a 46b 47b 47b 48a 48a 48a 48a 48a, 1st lie, math_string 9f, points 65b, and u 62b. 

3.6.2. Mobius transformations of infinitesimal cycles. Now we check that transformation of an infinitesimal cycle is 
an infinitesimal cycle again... 

24c (Infinitesimal cycle calculations 23g)+= (12b) <24b 24dc> 

Cll=C10.sl2similarity(a, b, c, d, es); 

cout <C "Image under SL2(R) of infinitesimal cycle has radius squared: " -C endl 
mathstring <C Cll.det(es).subs(sl2-relationl, 
subs-options::algebraic | subs-options::no-pattern).subs(signs-Cube, 
subs-options :: algebraic \ subs-options :: no-pattern) .series{ epsilon=0,3). normalf) 
mathstring <C endl 

-C "Image under cycle similarity of infinitesimal cycle has radius squared: " -C endl 
mathstring <C C10.cyclesimilarity(C, es).det(es).subs(signs-cube, subs-options::algebraic 
| subs-options::no-pattern).series(epsilon=0,3).normal() mathstring <C endl, 

Uses cycle 2a 2b 3e 3e 3e 3e 3e 46a 46b 47b 47b 48a 48a 48a 48a 48a and math_string 9f. 

... and focus of the transformed cycle is (up to infinitesimals) obtained from the focus of initial cycle by the same 
transformation. 

24d (Infinitesimal cycle calculations 23g)+= (12b) <24c 24ec> 

displ = ( Cll.focus(e4, true). subs{nval) - gW.subs{smap, subs-options::no-pattern)).evalm()', 
cout <C "Focus of the transormed cycle is from transformation of focus by: " 
mathstring <C displ.subs(sl2-relation, subs-options::algebraic 

| subs-options::no-patterri).subs(lst(sign=0,a=(l+b*c)+d)).series(epsilon^0,2).normal() 
mathstring <C endl, 

Uses cycle 2a 2b 3e 3e 3e 3e 3e 46a 46b 47b 47b 48a 48a 48a 48a 48a, 1st lie, and math_string 9f. 

3.6.3. Orthogonality with infinitesimal cycles. We also find expressions for the orthogonality (see § 3.3) with the 
infinitesimal radius cycle. 

24e (Infinitesimal cycle calculations 23g)+= (12b) <24d 24fc> 

cout <C "Orthogonality (leading term) to infinitesimal cycle is:" <C endl <C wspaces 
mathstring <C ex(C.is-orthogonal(C 10, es)).series(epsilon=0,l).normal () mathstring <C endl; 

Uses cycle 2a 2b 3e 3e 3e 3e 3e 46a 46b 47b 47b 48a 48a 48a 48a 48a, ex 4a lib 39d 48a 48b 66a 66b 66c 67a, math_string 9f, 
and wspaces 9f. 

And the both expressions for the s-orthogonality (see § 3.4) conditions with the infinitesimal radius cycle. The second 
relation verifies the Lem. 6.5 from [10]. 

24f (Infinitesimal cycle calculations 23g)+= (12b) <24e 25c i> 

cout <C "s-Orthogonality of other cycle to infinitesimal :" <C endl <C wspaces 
mathstring <C C.iss-orthogonal(C10, es).series(epsilon=0,l).normal() mathstring <C endl 
"s-Orthogonality of infinitesimal cycle to other:" -C endl <C wspaces 
mathstring <C C10.iss-orthogonal(C, es).series(epsilon=0,3).normal() mathstring <C endl; 


Uses cycle 2a 2b 3e 3e 3e 3e 3e 46a 46b 47b 47b 48a 48a 48a 48a 48a, math_string 9f, and wspaces 9f. 


3.6.4. Cayley transform of infinitesimal cycles. Here is two matrices which defines the Cayley transform and its 
inverses: 

25a (Declaration of variables lOd}+= (10a) «23d 

const matrix TC(2,2, 1st (dirac-ONEf), -e.subs(mu= 1), signl*e.subs{mu= 1), dirac-ONEQ)); 

II the inverse is TCI(2,2, lst(dirac_ONE(), e.subs(mu==l), - sign 1 *e. subs (rrm==1), dirac_ONE())); 


Defines: 

matrix, used in chunks 2d, 4g, 5b, 16d, 18b, 25d, 38a, 40-44, 50-56, and 68. 

Uses 1st lie. 

We conclude with calculations of the parabolic Cayley transform [10, § 8.3] on infinitesimal radius cycles. The 
parabolic Cayley transform on cycles is defined by the following transformation. 

25b (Parabolic Cayley transform of cycles 25b) = (12b) 

cycle2D cayley-parab (const cycle2D & C, const ex & sign = -1) 

{ 

return cycle2TD(C.get-k()~2*sign* C.getJ(l), C.getJQ, C.get-m()-2*C. getJ(l), C.get-uniti))-, 

} 


Uses cycle2D 6g 11c 11c 34b 35b 39d 39d 39d 39d 39d 48b 48b 48b 48b 48b 56b 56b 57d and ex 4a lib 39d 48a 48b 66a 66b 66c 67a. 
The image of an infinitesimal cycle is another infinitesimal radius cycle... 

25c (Infinitesimal cycle calculations 23g)+= (12b) <24f 25dc> 

Cll = cayley-parab(C 10, signl ); 

cout <C "Det of Cayley-transformed infinitesimal cycle: " 

math_string <C Cll.detf). sub s(\st (sign = 0), 

subs-options::algebraic \ subs-options\:no-patterri).series(epsilon=0,3).normal() 
math_string <C endl; 

Uses cycle 2a 2b 3e 3e 3e 3e 3e 46a 46b 47b 47b 48a 48a 48a 48a 48a, 1st lie, and math_string 9f. 

... with its focus mapped by the Cayley transform. 

25d (Infinitesimal cycle calculations 23g)+= (12b) <25c 25ec> 

displ = ( Cll.focus(e4, true). subs(nval) 

- clifford-moebius-map(TC, matrix(2,1,1st (u,vp)), e)).evalm().normal (); 
cout <C "Focus of the Cayley-transformed infinitesimal cycle displaced by: " math_string ; 
try{ 

cout -C displ.subs(lst(sign = 0), 

subs-options::algebraic \ subs_options::no-pattern).series(epsilon=0, 2).normal()', 

} catch (exception &p) { 
cout "(" -C displ.op{Q).subs{\st[sign = 0), 

subS-options::algebraic \ subs-options::no-patterri).series(epsilon=0 , 2 ).normal() 

< ", " displ.op(l).subs(lst(sign = 0), 

subs-options::algebraic \ subs-options::no-patterri).series(epsilon=0, 2).normal() 

> 


Defines: 

catch, used in chunks 42d and 50a. 

Uses cycle 2a 2b 3e 3e 3e 3e 3e 46a 46b 47b 47b 48a 48a 48a 48a 48a, 1st lie, math_string 9f, matrix 8c lOe Ilf 12a 25a, and u 62b. 
s-Orthogonality of 

25e (Infinitesimal cycle calculations 23g)+= (12b) < 25d 

cout math-string <!C endl 

-C "s-Orthogonality of Cayley transforms of infinitesimal cycle to other:" <C endl <C wspaces 
math_string <C Cll.iss-orthogonal(cayley-parab(C,signl), es).series(epsilon=0,3).normal() 
math_string <C endl <C endl, 

> 


Uses cycle 2a 2b 3e 3e 3e 3e 3e 46a 46b 47b 47b 48a 48a 48a 48a 48a, math_string 9f, and wspaces 9f. 


3.7. Drawing the Asymptote output. Although we use every possibility above to make double and cross checks one 
may still wish to see “by his own eyes” that the all calculations are correct. This may be done as follows. 

We draw some Asymptote pictures which are included in [10], see also Fig. 1. We start from illustration of the both 
orthogonality relations, see § 3.3 and 3.4. They are done for nine (=3x3) possible combinations of metrics (elliptic, 
parabolic and hyperbolic) for the space of points and space of cycles. 

26a (Draw Asymptote pictures 26a) = (10c) 26b > 

ofstream as?/mpfoie("parab-orthol .asy"); 
asymptote <C setprecision( 2); 
for (si = -1; si < 2; si++) { 
for (sil = -1; sil < 2; sil++ ) { 
sigri-val = lst(sj^n = si, signl = sil)', 

Uses 1st lie, si lOe, and sil lOe. 

For each of those combinations we produce pictures from the set of data which is almost identical. This help to see 
the influence of sign and signl parameters with constant other ones. All those graphics are mainly application of 
asy-draw() method (see § 2.6 mixed with some Asymptote drawing instructions. Since this is rather technical issue we 
put it separately in Appendix C. 

26b (Draw Asymptote pictures 26a)+= (10c) «26a 26c c> 

try { 

{(Drawing first orthogonality 30a)} 

{(Drawing second orthogonality 32a)} 

} catch ( exception &p) { 

cerr <C "***** Got problem2: " -C p.what() -C endl; 

} 

} 

} 


Defines: 

catch, used in chunks 42d and 50a. 

We finish the code with generation of some additional pictures for the paper [10]. 

26c (Draw Asymptote pictures 26a)+= (10c) «26b 

try { 

(Extra pictures from Asymptote 32b) 

} catch (exception Szp) { 

cerr <C "***** Got problem3: " -C p.what() -C endl; 

} 

asymptote. close(); 


Defines: 

catch, used in chunks 42d and 50a. 

4. How to Get the Code 

(i) Get the DT}]X source of this paper [12] from the arXiv.org. 

(ii) Run the source through DT}]X. Five new files will be created in the current directory. 

(a) noweb [14] sources. 

(b) Header file cycle.h of the library. 

(c) C++ source cycle. cpp of the library. 

(d) C++ source parab-orthol. cpp of the example. 

(e) Asymptote source of the graphics. 

(iii) Use it on your own risk under the GNU General Public License [4]. 

To obtain the Python wrapper [9] for this library proceed to the link http: //maths. leeds . ac .uk/~kisilv/pycycle .html 
and follow the instructions given there. 


Appendix A. Textual output of the program 


Conjugation of a cycle comes through Moebius transformation: true 
A K-orbit is preserved: true, and passing (0, t): true 

Determinant of zero-radius Z1 cycle in metric e is: —v 2 a + av 2 
Focus of zero-radius cycle is: it, \va — \av 
Centre of zero-radius cycle is: it, —va 
Focal length of zero-radius cycle is: 

Image of the zero-radius cycle under Moebius transform has zero radius: true 
The centre of the Moebius transformed zero-radius cycle is: -equal-, -equal- 
image of the zero-radius cycle under cycle similarity has zero radius: true 
The centre of the conjugated zero-radius cycle coinsides with Moebius tr: -equal-, -equal- 
The orthogonality is: i km + nan + ^mk — ll == 0 


The orthogonality of two lines is: nan — ll == 0 

The orthogonality to z-r-cycle is: — ul + ^ u 2 k — ^kv 2 a + anv == 0 

The orthogonality of two z-r-cycle is: v'av — ^v' 2 x(& 2 ) — jdn 2 + \u 2 — uu' — \u' 2 == 0 
Both orthogonal cycles (through one point and through its inverse) are the same: true 
Orthogonal cycle passes through the transformed point: true 
Line through point and its inverse is orthogonal: true 
All lines come through the point (I , — s r L ) 

Conjugated vector is parallel to (u,v): true 

C5 has common roots with C : true 

x(er)-centre of C5 is equal to d-centre of C: true 

Inversion in (C5, sign) coincides with inversion in (C, signl): true 

Inversion to the real line (with - sign): 

Conjugation of the real line is the cycle C: true 
Conjugation of the cycle C is the real line: true 
Inversion cycle has common roots with C: true 
C passing the centre of inversion cycle: true 
Inversion to the real line (with + sign): 

Conjugation of the real line is the cycle C: true 
Conjugation of the cycle C is the real line: true 
Inversion cycle has common roots with C: true 
C passing the centre of inversion cycle: true 
Yaglom inversion of the second kind is three reflections in the cycles: true 
The real line is Moebius invariant: true 

Reflection in the real line: (1, ( u -v ) symboi2696 , ~v 2 a + u 2 ) 

Reflection of the real line in cycle C: 

(2x(a 3 )akx(a 2 )n, ( 2 X (a 3 )lax(a 2 )n (dx(c r 2 )n 2 + Z 2 x(cr 2 ) - mkx(a 2 ))x(a 2 ) ) symbom5S , 2mx(a 3 )ax(a 2 )n) 


The s-orthogonality is: ( kx{a 2 )nrh — 2llx(a 2 )n — mnkx(a 2 ) + mkx(a 2 )n + l 2 nx( l J 2 ) + nax(a 2 )n 2 )x(o' 2 ) == 0 
The s-orthogonality of two lines is: x( a 2)(—2llx(& 2 )n + l 2 nx{a 2 ) + nax(a 2 )n 2 ) == 0 
The s-orthogonality to z-r-cycle is first way: 

X ( C7 2 ){—mkx{cr 2 )v + l 2 x(. a 2 )v — 2ulx(a 2 )n + u 2 kx(a 2 )n + mx{a 2 )n + <Tx(u 2 )n 2 D — akx(a 2 )nv 2 ) == 0 
The s-orthogonality to z-r-cycle is second way: 
x( tJ 2 )(— &kx(a 2 )v 3 + mx(a 2 )v + u 2 kx{a 2 )v — 2ulx(a 2 )v + 2ax(a 2 )nv 2 ) == 0 
The s-orthogonality of two z-r-cycle is: 

(2v'ax(& 2 )v 2 + u' 2 x(& 2 )v - 2uu'x{a 2 )v + u 2 x{& 2 )v - ax{cr 2 )v’ i - v' 2 x{^ 2 )va)x{^ 2 ) == 0 
All lines come through the focus related e: true 
C8 has common roots with C : true 
x(o')-center of C8 coincides with d-focus of C : true 
s-Inversion in C coincides with inversion in C8 : true 
Distance between (u,v) and (u’,v’) in elliptic and hyperbolic spaces is 
{—a(—v' + v ) 2 + (it — ii') 2 )(i/(4 — 4aa)v + (— a(—v' + v) 2 + (u — u') 2 )a) 


true 


—(—v 1 + v) 2 + (ii — n ') 2 d 
Conformity in a cycle space with metric: E P H 
Point space is Elliptic case (sign = -1): true false false 
Point space is Hyperbolic case (sign = 1): false false true 
Perpendicular to ((u,v); (u\v’)) is: -« 3 g-2«V»*+»V+« ,2 »+i/^ 7 ^^ 

— 2v' 2 u' — u' a v 2 cr-\-u av 2 a—2uv'v — u 3 a-\-2v' u' v—uv' 2 a <j — 3uu' 2 a+v' 2 u' <rcr-\-3u 2 u' a-\-2uv' 2 -\-u' 3 a 

— 4v'v-\-4uu'<j— 2u 2 <7+2 v’ 2 -\-2v r2 — 2u' 2 a 

Value at the middle point (parabolic point space): 
u 2 — 2 uu' + u' 2 

Conformity in a cycle space with metric: E P H 





Point space is Parabolic case (sign = 0): true true true 
Perpendicular to ((u,v); (u’,v’)) is: v'a; — \u! 
Distance between (u,v) and (u’,v): 

Value at critical point: 

— 2uu' a+u 2 a+Av 2 +u 12 ft—Aav 2 a 


Length from * center* between (u,v) and (u’,v’): 

2»2 o / o 2 /2 =2 i ' 2«2 

u a —Zuu a —v a a—av -\-Zv av-\-u a 


This distance/length is conformal: true 
Perpendicular to ((u,v); (u’,v’)) is: 


Length from *focus* check for p = b(^J—v l2 ba + ( v' — v) 2 + (—u + u') 2 1 1 — v' + v) 

Length between (u,v) and (u’, v’) is equal to ((a) — a)p 2 — 2 vp: true 
checks: Cll passes through (u’, v’): true; Cll focus is at (u, v): true 
This distance/length is conformal: true 
Perpendicular to ((u,v); (u’,v’)) is ( av' + p,u — u')\ true 

Length from *focus* check for p = b(—\J—v' 2 bo + (v' — v) 2 + {—u + u') 2 b — v' + v) 

Length between (u,v) and (u’, v’) is equal to ((cr) — b)p 2 — 2 vp: true 
checks: Cll passes through (u’, v’): true; Cll focus is at (u, v): true 
This distance/length is conformal: true 
Perpendicular to ((u,v); (u’,v’)) is (av' + P,u — u'): true 
Length from *focus* check for p = 5 ~ v a ~l!~Iy +u ^ 

Length between (u,v) and (u’, v’) is equal to ((cr) — b)p 2 — 2 vp: true 
checks: Cll passes through (u’, v’): true; Cll focus is at (u, v): true 
This distance/length is conformal: false. The factor is: 


(c 2 yv 2 a—2uc 2 vx-\-u 2 c 2 y-\-2udcy-\-d 2 y—2dcvx) 2 

Perpendicular to ((u,v); (u’,v’)) is ( av' +p,u — u'): true 
Inf cycle is: (1, ( u yv mb ° 13712 yA-y/-(-*+W+vl+v,) _ + ^ 

Square of radius of the infinitesimal cycle is: — e 2 
Focus of infinitesimal cycle is: u, v p 


Focal length is: (±d-)e 2 + 0(e 3 ) 


Infinitesimal cycle passing points(u + ex, vp + ( v p x 2 ) + (1 o+ax ^2 + 


/ 1 V 4 

Image under SL2(R) of infinitesimal cycle has radius squared: 


3^4 r -2 - 2-2 


/_ 4<7<7+4<t<t — a —ber a — a _ \ 2 1 /n/_3\ 

' (u 2 c 2 a 2 -\-u 2 c 2 a 2 — 2u 2 c 2 aa-\-d 2 a 2 -\-2udca 2 -2d 2 aa—4udcaa-\-2udca 2 -\-d 2 a 2 ) 2 ' ' ' 

Image under cycle similarity of infinitesimal cycle has radius squared: 

( 4m 2 ak 2 a 3 + 121 2 aa 2 n 2 —8ml 2 aka+Al 4 aa+2l 2 an 2 —2maka 4 n 2 —8l 2 a 2 an 2 +4aa 3 n 4 +2ml 2 ka A +4rn 2 ak 2 a+8rna 2 ka 3 n 2 — l 4 a 2 —a 2 n 4 —m 2 k 2 a 4 —2makn 2 —m 2 c 


0(e 3 ) 


( u 2 a 2 k 2 —an 2 — 2l 2 a a—2ula 2 k-\-l 2 a 2 -\-4ulaka — da 


Focus of the transormed cycle is from transformation of focus by: ( 


0 


) + ( 


)e + 0(e 2 ) 


Orthogonality (leading term) to infinitesimal cycle is: 

(^m — ul + \u 2 k == 0 ) + 0(e ) 
s-Orthogonality of other cycle to infinitesimal: 

(' u 2 kn — 2 uln + mn == 0) + 0(e) 
s-Orthogonality of infinitesimal cycle to other: 

(0 == 0 ) + (0 == 0 )e + (^( == 0 ))e 2 + 0(e 3 ) 

2 - _ 

Dct of Cayley-transformed infinitesimal cycle: ( v ^ ,p )e 2 + 0(e 3 ) 

Focus of the Cayley-transformed infinitesimal cycle displaced by: (0(e 2 ), 0(e 2 )) 


s-Orthogonality of Cayley transforms of infinitesimal cycle to other: 
(0 == 0 ) + (0 == 0 )e + a ^-2ul+n 2 k-2v p n == Q) y + ^3) 
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U - i~r+u‘) 

Square of radius of the infinitesimal cycle is: —e 2 
Focus of infinitesimal cycle is: u, v p 
Focal length is: (|^)e 2 

Infinitesimal cycle passing points(w + ex, vp + ( v p x 2 ) + (-j-fj-)e 2 ), 
Image under SL2(R) of infinitesimal cycle has radius squared: 
(-(_ u 2 c 2 _rf 2 _ 2 . udc) 2 )e 2 + 0(e 3 ) 

Image under cycle similarity of infinitesimal cycle has radius squared: 
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( (an 2 —l 2 +2ulk—u 2 k 2 ) 2 

Focus of the transormed cycle is from transformation of focus by: ( 

Orthogonality (leading term) to infinitesimal cycle is: 

(|m — ul + i u 2 k == 0) + O(e) 
s-Orthogonality of other cycle to infinitesimal: 

(u 2 kn — 2 uln + mn == 0 ) + O(e) 
s-Orthogonality of infinitesimal cycle to other: 

(0 == 0 ) + (0 == 0)e + (i( m - 2 ^+« 2fc - 2 ^» == 0 ))e 2 + 0(e 3 ) 

Det of Cayley-transformed infinitesimal cycle: ( - — Lp )e 2 + 0(e 3 ) 


Focus of the Cayley-transformed infinitesimal cycle displaced by: ( 

s-Orthogonality of Cayley transforms of infinitesimal cycle to other: 
(0 == 0 ) + (0 == 0 )e + (i( m ~ 2 ^+“ 2fc - 2 ^" == Q))e 2 + <D(e 3 ) 


) + ( 


)e + 0(e 


)e + 0(e 2 ) 


Appendix B. Example of the produced graphics 


An example of graphics generated by the program is given in Figure 1. This was produced by the part of program 
from the Section C.1.1. 






Figure 1. Orthogonality of the first kind in nine combinations. 


























Appendix C. Details of the Asymptote Drawing 
C.l. Drawing Orthogonality Conditions. 

C.1.1. First Orthogonality Condition. We define numeric values of all involved parameters first. 

30a (Drawing first orthogonality 30a) = (26b) 

numeric a;mm(-ll,4), xmax( 5), ymin(- 3), ymax = (si = 0?numeric(25, 4): 4); 

1st cycle-val = 1st (sign = numeric(si), signl = numeric(sif), 

k= numeric(2,3), 1= numeric(2,3), n = (si = l?numeric(-l):numeric(l,2)), m =numeric(-2)); 
cycle2D Cf = C.subs(cycle-val), Cg = C5.subs(cycle-val ), Cp =C2\ 

1st U, V: 
switch (si) { 

case -1: // points b, a, center, c, d 

U= numeric(ll,4), Cg.roots(half).op( 0), Cf.center().op(0).subs(cycle-val ), (l-^k).subs(cycle-val); 

V= Cf.roots(U.op(0 ), false). op(l), half, Cf.center().op(l).subs(cycle-val ), 

Cf.roots(l-h-k , false). op(0) . normal() .subs(cycle-val ); 

break; 
case 0: 

U = numeric(17,4), Cg.roots().op(0), Cf.center().op(0).subs(cycle-val ), ( Prk).subs(cycle-val ); 

V = Cf.roots(U.op(0), false).op(0), numeric(3,2), Cf.roots(l+k , false). op(0).subs(cycle-val), 

Cf.roots( (4- k, false). op(0) . normal(). subs( cycle-val ); 

break; 
case 1: 

U = numeric(12,4), Cg.roo£s(numeric(3,4)).op(0), Cf.center().op(0).subs(cycle-val), (l-trk).subs(cycle-val); 

V= Cf.roots(U.op(0), false).op(0), numeric(3,4), Cf.center().op(l).subs(cycle-val), 

Cf.roots( l-r- k , false). op(0) . normal () .subs( cycle-val ); 

break; 

} 

U. append(P.op(0).subs(cycle-val).subs(lst(u= U.op( 0), v= V.op(0))) .normal())' //Moebius transform of the first point 

V. append(P.op(l).subs(cycle-val).subs(lst(u = U.op(0 ), v= V.op(0))).normal()); 

asymptote -C endl <tC "erase ();" <C endl ; 

(Drawing orthogonal cycles 30b) 

asymptote -C "shipout (\"f irst-ort-" <C eph_names[si+ 1] <C eph_names[sil+ 1] <C -C endl ; 

Uses cycle2D 6g 11c 11c 34b 35b 39d 39d 39d 39d 39d 48b 48b 48b 48b 48b 56b 56b 57d, eph_names lOd, 1st lie, numeric lOd, 
points 65b, si lOe, sil lOe, u 62b, and v 62b. 

We start drawing from cycles. 

30b (Drawing orthogonal cycles 30b) = (30a 32a) 31a c> 

for (int j = 0; j< 2; j++) 
for (int i=0; i<(si= 1?4:5); *++) 

Cp.subs(lst(kl = (si= 0? numeric(3* *,2): numeric)*, 4)), nl = half, u= U.op(j), 
v = V.op(j))).subs(cycle-val).asy-draw(asymptote, xmin, xmax, ymin, ymax, 
lst(0.2, 0.2+j*(0.3+H-8.0), 0.2+(l-.?>(0.3+*4-8.0))); 

Cf.asy-draw(asymptote, xmin, xmax, ymin, ymax, lst(0.8, 0, 0), "+1"); 

Cg.asy-draw(asymptote, xmin, xmax, ymin, ymax, lst(0, 0, 0), "+0.3+dashed"); 
if (si = 0) 

C5.sub s(lst (sign =0, signl=0)).subs(cycle-val).asy-draw(asymptote, xmin, xmax, ymin, ymax, lst(0, 0, 0), 
"+dotted"); 


Uses 1st lie, numeric lOd, si lOe, u 62b, and v 62b. 


To finish we add some additional drawing explaining the picture. 

3la (Drawing orthogonal cycles 30b) += (30a 32a) < 30b 

asymptotes "pair[] z={(" -C ex_to<numeric>( U.op(0).evalf()).tO-double() S", " 
<C ea;_to<numeric>( V.op(Q).evalf(j).to-double() S ")*u"; 
for (int j = 1; j< 5; ]++) 

asymptotes ", (" <C ex-to<numeric>(U.op(j).evalf()).to-double() S ", " 

S ex-to<numeric>(V.op(j).evalf()).to-double() S ")*u" ; 

asymptote S " <C endl S " dot(z); " <C endl 

S (si = 0? " draw((z[2].x,0)—z [2], 0.3+dotted )S endl 
S (si = 0? " draw((z[3].x,0)—z [3], 0.3+dotted )S endl 
S" label(\"$a$\", z[l], NW )\" S endl 
S" label (\"$b$\" , z[0], SE )S endl 
S" label(\"$c$\", z [3], E) ;" < endl 

S" label" < "(\"$d$\", z [4] , " < (si =1?"NW) ; ":"NE) ; ") < endl] 

(Put units 3lc) 

(Draw axes 31b) 


Uses numeric lOd, si lOe, and u 62b. 

This chunk draws the standard coordinat axes. 

3lb (Draw axes 31b) = (31a 33-37) 

asymptote " draw_axes((" <C xmin.to-double () S " , " S ymin.to-double() 

•C "), ( " <C xmax.to-double() S ", " <C ymax.to-double() S <C endl, 

31c (Put units 31c) = (31a 36d) 

asymptote <C " label(\"$\\sigma=" S si S " , \\breve{\\sigma)-=" <C sil 

(0, " S ymin.to-double() S S); " <C endl S "draw( (1,-0.1) *u—(l,0.1)*u ) S endl 

S "draw((-0.1,1)*u—(0.1,l)*u)< endl 
S "label(\"$1$\", (1,0)*u, S) ; " < endl 
S "label (\"$1$\" , (0,1) *u, E) ; " < endl- 


Uses si lOe, sil lOe, and u 62b. 


C.1.2. Second Orthogonality Condition. We draw some Asymptote pictures to illustrate the second orthogonality 
relation. We define numeric values of all involved parameters first. 

32a (Drawing second orthogonality 32a) = (26b) 

numeric a:mm(-ll,4), xmax( 5), ymin(- 13,4), ymax = (si = 0?numeric(6): numeric(15,4)); 

1st cycle-val = 1st (sign = numeric(si), signl = numeric(sD), sign2 = numeric(l), //sign3 == jump_fnct(- 
si), //sign3 == (si > 0?numeric(-l):numeric(l)), 

k= numeric(2,3), 1= numeric(2,3), n= ( si= l?numeric(-4,3):/ia(/), m=(si= l?numeric(-9,3):numeric(- 

2 ))); 

cycle2D Cf = C.subs(cycle-val), Cg = C8.subs(cyclejval ), Cp =C6; 

1st U, V; 
switch (si) { 

case -1: // points b, a, center, c, d 

U = numerical,4), Cg.roots(half).op( 0), Cffocus().op(0).subs(cycle-val ), (l-^-k).subs(cycle-val); 

V= Cf.roots(U.op(0 ), false).op(l), half; Cf.focus().op(l).subs(cycle-val ), 

C1.roots(l-r-k , false).op(0). normal() .subs(cyde-vat ); 

break; 
case 0: 

U = numeric(4), Cf.roots().op( 0), Cf.focus().op(0).subs(cyde-val), (l+k).subs(cyde-val)- 1 
V= Cf.roots(U.op(0), false).op(0), numeric(3,2), Cf.focus().op(0).subs(cyde-val), 

Cl. roots( 1 -t- k, false). op(0). normal(). subs( cycle-val ); 

break; 
case 1: 

U = Cf.roots(numeric(l)).op(l ), Cg.roots( numeric(6, 4)).op(l), 

Cf.focusQ.op(0).subs(cyde-val ), (l+k).subs(cyde^val ); 

V = numeric(l), numeric(6, 4), Cf.focusQ.op(l).subs(cyde-vat), 

Cl. roots( Z-7- A:, false). op(0). normalQ .subs( cycle-val ); 

break; 

> 

U. append(Pl.op(0).subs(cycle-val).subs(\st(u= U.op( 0), v= V.op(0))). normalQ); //Moebius transform of U.opQf) 

V. append(Pl.op(l).subs(cycle-val).subs(lst(u = U.op( 0), v= V.op(0))).normalQ); 

asymptote endl <C "erase ();" //<< endl << ”size(250);” 

<C endl; 

(Drawing orthogonal cycles 30b) 

asymptote <C "shipout(\"sec-ort-" <C eph_names[si+ 1] <C eph_names[sil+ 1] -C -C endl; 

Uses cycle2D 6g 11c 11c 34b 35b 39d 39d 39d 39d 39d 48b 48b 48b 48b 48b 56b 56b 57d, eph_names lOd, jump_fnct 38a, 1st lie, 
numeric lOd, points 65b, si lOe, sil lOe, u 62b, and v 62b. 

C.2. Extra pictures from Asymptote. We draw few more pictures in Asymptote. 

32b (Extra pictures from Asymptote 32b) = (26c) 

numeric xmin(-5 ), xmax(5 ), ymin(- 13,4), ymax = numeric(6); 

(Three images of the same cycle 33a) 

(Centres and foci of parabolas 33b) 

(Zero-radius cycle implementations 34a) 

(Parabolic diameters 34b) 

(Distance as an extremum 35a) 

(Infinitesimal cycles draw 35b) 

(Cayley transform pictures 35c) 

(Three inversions 36e) 


Uses numeric lOd. 


C.2.1. Different implementations of the same cycle. A cycle represented by a four numbers (k,l,n,m looks different 
in three spaces with different metrics. 

33a (Three images of the same cycle 33a)= (32b) 

asymptote <C endl <C "erase ();" <C endl; 
cycle2D Clf C2f 
asymptote <C "pair[] z;"; 
for (int j = -1; j<2 ; j++) { 

Clf = cycle2D(l, lst(-2.5, 1), 3.75, diag-matrix{\st{-\, j))); 

C2f = cycle2D(l, lst(2.75, 3), 14.0625, diag-matrix{ lst(-l, j))); 

Clf.asy-draw(asymptote, xmin, xmax, ymin, ymax, lst(0, 1.0-0.4*0+1), 0.4*0+1)), "+.75", true, 7); 
C2f.asy-draw(asymptote, xmin, xmax, ymin, ymax, lst(0, 1.0-0.4*0+1), 0.4*0+1)), "+.75", true, 7); 
asymptote <C "z.push((" -C Clfcenter().op( 0) <C ", " -C Clf.center().op(l ) -C ")*u); z.push((" 

< C2f.center().op( 0) -C ", " -C C2f.center().op( 1) <C ")*u);" <C endl, 

} 

asymptote<€. "z.push((" <C Clf.roots().op( 0) <C " , 0)*u); z.push((" <C CIf roots f).op( 1) <C ", 0)*u);" <C endl 
-C " dot(z);" <C endl 

" for (int j = 0; j<2; ++j) {" 

•C " label(\"$c_e$\" , z[j], E) ; " <C endl 

<" label(\"$c_p$\", z[j+2] , SE);" < endl 

<C " label(\"$c_h$\", z[j+4], E);" <C endl 

<" label( (j==0?\"$r_0$\", z[j+6], (j==0? SW: SE));" < endl 
•C " draw(z[j]—z[j+4], . 3+dashed) ; " endl 

<C " >" <C endl, 

(Draw axes 31b) 

asymptote -C "shipout(\"same-cycle\")-C endl, 

Uses cycle 2a 2b 3e 3e 3e 3e 3e 46a 46b 47b 47b 48a 48a 48a 48a 48a, cycle2D 6g 11c 11c 34b 35b 39d 39d 39d 39d 39d 48b 48b 48b 48b 
48b 56b 56b 57d, 1st lie, and u 62b. 

C.2.2. Centres and foci of cycles. We draw two parabolas and their centres with three type of foci. 

33b (Centres and foci of parabolas 33b) = (32b) 

asymptote <C endl <C "erase (); " <C endl, 

Clf = cycle2D(l, lst(-1.5, 2), 3.75, par_matr); 

C2f = cycle2D(l, lst(2, 2), -3.5, par-matr ); 

Clf.asy-draw(asymptote, xmin, xmax, ymin, ymax, lst(0, 1.0-0.4, 0.4), "+.75", true, 7); 

C2f.asy-draw(asymptote, xmin, xmax, ymin, ymax, lst(0, 1.0-0.4, 0.4), "+.75", true, 7); 

asymptote -C "pair [] z= {(" <C Clf.center(-unit-matrix(2)).op(0) <C " , " <C Clfcenter(-unit-matrix(2)).op(l) 

< ")*u, (" <C C2f.center(-unit_matrix(2)).op(0 ) <C ", " <C C2f.center(-unit-matrix(2)).op(l ) <C ")*u, 
for (int j = -1; j<2 ; j++) { 
ex MS = diag-matrix( lst(-l, j )); 

1st FI = ex-to<\st>(Clf.focus(MS)), F2 = exJ,o< 1st >(C2f.focus(MS)); 
asymptote -C " (" Fl.op( 0) -C " , " <C Fl.op( 1) <C ")*u, (" 

<C F2.op{ 0) <C ", " -C F2.op( 1) -C ")*u" <C 0=1? endl; 

} 

asymptote -C " dot (z); " <C endl 

•C " draw(z[0]— z[l], dashed);" <C endl; 


asymptote -C "for (int 

j=l; j<3; 

++j) {" <C endl 

< 

" label(\"$c_e$\" 

, z[j-l] , 

N) ;" <C endl 

< 

" label(\"$f_e$\" 

, z[j+l] , 

E);" <C endl 

< 

" label(\"$f_p$\" 

1-1 

co 

+ 

•r-> 

N 

E);" -C endl 

< 

" label(\"$f_h$\" 

1—1 
LO 
+ 
•r-> 

N 

E);" -C endl 

< 

" draw(z[j+l] — z[j+5], dotted+0. 5);" <C endl 

< 

"}" <C endl; 




(Draw axes 31b) 

asymptote <C "shipout (\"parab-cent\");" <C endl; 

Uses cycle2D 6g 11c 11c 34b 35b 39d 39d 39d 39d 39d 48b 48b 48b 48b 48b 56b 56b 57d, ex 4a lib 39d 48a 48b 66a 66b 66c 67a, 1st lie, 
parjnatr 9e, and u 62b. 







C.2.3. Zero-radius cycles. Zero-radius cycles can look different in different EPH realisations, here is an illustration. 

34a (Zero-radius cycle implementations 34a)= (32b) 

asymptote <C endl <C "erase ();" <C endl 
-C "pairG z;" <C endl; 

{ 

numeric xmin(- 5), xmax{ 15), ymin{- 5), ymax( 5); 
for (int il=- 1; il< 2; il++) { 
for(int i2=- 1; i2<2; ii?++) { 

1st val(sign=il, signl=i2 , w=6*il+4, u=1.7); 

Zl.subs{val).asy-draw(asymptote, xmin, xmax, ymin, ymax , 1st(0.5+0.4*if, .5-0.3*i£, 0.5+0.3*iD),"", true, 7); 
asymptote <C "dot((" -C ex-to<num.eric>(Zl.focus(e).op(0).subs(val)).to_double() 

<C ", "<C ex-to<numeric>(Zl.focus(e).op(l).subs(val)).to-double() 

< ")*u, " < 0.4+0.4*il < "red+" 

<C .4-0.3*ii?<C "green+" 

<C 0.6+0.3*i£ <C "blue);" <C end/; 

} 

} 

(Draw axes 31b) 

> 

asymptote -C "shipout(\"zero-cycles\");"<C endl; 


Uses 1st lie, numeric lOd, u 62b, and v 62b. 

C.2.4. Diameters of cycles. The notion of diameter and related distance became strange in parabolic case. 

34b (Parabolic diameters 34b) = (32b) 

asymptote -C endl <C " erase ();" <C endl; 

CIO = cycle2D(l, lst((-4-l)-^2.0, 0.5), 4 ,par_matr); 

C10.asy-draw(asymptote, xmin, xmax, ymin, ymax , lst(0.1, 0, 0.6)); 

asymptote <C "pair [] z = {(" <C C10.roots().op(0) <C ", 0)*u, (" <C ClO.roolsQ.op(l) <C ", 0)*u> ; " <C endl; 
cycle2D(l, lst(5-^2.0, 0.5), 8 ,par_matr).asy-draw(asymptote, xmin, xmax, ymin, ymax, 
lst(0.1, 0.6, 0), true, 7); 

CIO =cycle2D(-l, lst(-5-^2.0, 0.5), 8-5.0*5-h2.0,par_mafr); 

C10.asy-draw(asymptote, xmin, xmax, ymin, ymax, lst(0.1, 0.6, 0), 

"+dashed ", true, 7); 

asymptote <C "z ,push( (" <C C10.roots().op( 1) <C ", 0) *u) ; z .push( (" <C C10.roots().op( 0) <C " , 0)*u); " <C endl; 
(Put labels on 22-23 34c) 

(Draw axes 31b) 

asymptote <C "shipout(\"parab-diam\") <C endl; 

Defines: 

cycle2D, used in chunks 6f, 7a, 12c, 14-18, 20a, 21g, 24a, 25b, 30a, 32, 33, 35-37, 39-41, 55-57, and 62-64. 

Uses 1st lie, par_matr 9e, and u 62b. 

Here is the common part of drawing points and labels on the figures 22-23. 

34c (Put labels on 22-23 34c) = (34b 35a) 

asymptote <C "z.push((z[2].x,0)); z,push((z[3],x,0));"<C endl 


< " 

dot(z);" -C endl 




< " 

draw(z [2] —z [3] , 

black+ 

.3);' 

<C endl 

< " 

draw(z [0] —z [1] , 

black+1.2) 

" -C endl 

< " 

draw(z [4] —z [5], 

black+1.2) 

" < endl 

< " 

label(\"$z_l$\" 

z [0] , 

NW) ; 

" < endl 

< " 

label(\"$z_2$\" 

z[l] , 

SE) ; 

" <C endl 

< " 

label(\"$zJ3$\" 

z [2] , 

SW) ; 

" <C endl 

<C " 

label(\"$z_4$\" 

z [3] , 

SE) ; 

" <C endl; 




C.2.5. Extremal property of the distance. To illustrate the variational definition of the distance [12, Defn.5.2] we draw 
several cycles which passes two given points. The cycles with the extremal value of diameter is highlighted in bold. 

35a (Distance as an extremum 35a) = (32b) 

asymptote <C endl <C "erase ();" <C endl ; 
for (int j=- 2; j < 3; j++) { 

ex-to<cycle2T)>(C.subject-to(lst(C.passing(lst(xmin+l, ymax- 5)), C.passing(lst(xmin+3, ymax- 6.5)), k= 1, 

1= xmin+2+0.5*j)).subs(sign = -1 )).asy-draw(asymptote, xmin, xmax, ymin, ymax, 
lst(0, 0A*abs(j), 1.0-0 A* abs(j)), (j = 0 ? "+1" : "+.3")); 
ex-to<cycle2D>(C.subject-to(lst(C.passing(lst(xmax-4:, ymax- 5)), C.passmg(lst(a;maa;-1, ymax- 2)), fc= 1, 

Z = xmax-2.5-0.2*(j+2))).subs(sign = O)).asy-draw(asymptote, xmin, xmax, ymin, ymax, 
lst(0.2*(j+2), 0, 1.0-0.2*(jT2)), (j= -2 ? " + 1" : " + .3"), true, 7); 

} 

asymptote <C "pair[] z ={ (" <C xmin+1 -C ", " ymax -5 <C ")*u, (" -C min+3 <C ", " 

<C ymax- 6.5 <C ")*u, (" <C xmax-1 <C ", " <C ymax -5 <C ")*u, (" <C £ma:r-l 

<C ", " <C ymax-2 <C ") *u>;" <C endl, 

(Put labels on 22-23 34c) 

asymptote " label (\"$d_e$\" , . 5z [0] + . 5z [1] , NE) ; " <C endl 
< " label(\"$d_p$\" , . 5z [4] +. 5z [5] , S) ; " < endl, 

(Draw axes 31b) 

asymptote <!C "shipout(\"dist-extr\")<IC endl, 

Uses cycle2D 6g 11c 11c 34b 35b 39d 39d 39d 39d 39d 48b 48b 48b 48b 48b 56b 56b 57d, 1st lie, and u 62b. 

C.2.6. Infinitesimal cycles. Here we draw a set of parabola with the same focus and the focal length tensing to zero. 

35b (Infinitesimal cycles draw 35b) = (32b) 

asymptote <C endl <C "eraseO;" <C endl, 

for (int j= 1; j < 5; j++) { 

cycle2D(lst(-2.5, 4.5), -unitjmatrix{ 2), I6.0*pow(2, -2*j)).asy-draw(asymptote, xmin, xmax, ymin, ymax, 

1st(0, 0.2 *abs(j), 1.0-0.2 *abs(j)), "+.3"); 

cycle2D(lst(l, 1.25), hyp_matr, 25*pow(1.8, -2*j)).asy-draw(asymptote, xmin, xmax, ymin, ymax-^-3, 
lst(0.2 *abs(j), 1.0-0.2 *abs(j), 0), "+.3", true, 5 +j); 
cycle2D(l, lst(2, pow(3,-j)), 2*2+2.0*pow(3,-j)-pow(3,-2*j), par_matr).asy-draw(asymptote, xmin, 
xmax, ymin, ymax, lst(l.0-0.17*), 0, 0.17*)), "+.3", true, 7); 

} 

asymptote " draw((2,l)*u—(2," -C ymax ")*u, blue+1);" endl, 

cycle2D(lst(l, 1.25), hypjmatr).asy-draw(asymptote, xmin, xmax, ymin, ymax-^3, lst(l, 0, 0), " + 1"); 
asymptote <C " dot((-2.5,4.5)*u);" <C endl 
<C " dot ((2, l)*u); " <C endl; 

(Draw axes 31b) 

asymptote <C "shipout(\"infinites\")<C endl, 

Defines: 

cycle2D, used in chunks 6f, 7a, 12c, 14-18, 20a, 21g, 24a, 25b, 30a, 32, 33, 35-37, 39-41, 55-57, and 62-64. 

Uses hypjnatr 9e, 1st lie, par_matr 9e, and u 62b. 

C.2.7. Pictures of the Cayley transform. We draw now pictures of Cayley transform, which shows that the unit cycle 
UC may be obtained as a reflection of the real line into the cycle Cl Of. 

35c (Cayley transform pictures 35c) = (32b) 36a i> 

xmin = -numeric(4,2); a:ma:r=numeric(4,2); T/mm=-numeric(7,2); ymax= numeric(3); 
cycle2D Cl Of, UC; 

ClOf = cycle2D(l, lst(0, sign2), sign, e); 

U C= real-line, cycle-similarity (Cl Of, es). normalizef); 


Uses cycle2D 6g 11c 11c 34b 35b 39d 39d 39d 39d 39d 48b 48b 48b 48b 48b 56b 56b 57d, 1st lie, and numeric lOd. 


Now we run cycles over signatures of point and cycle spaces and sign of sign2. 

36a (Cayley transform pictures 35c) += (32b) «35c 36b > 

for (s?=-l; si<2; si++) { 
for (sil=- 1; sil<2; sil++) 
if ((si =0 ) V (si = sil )) { 
asymptote -C endl -C " erase ();" -C endl; 
for (int si2=- 1; si2<2; si2=si2+ 2) { 

1st cycle-val = 1st (sign = s*, signl = sil , sign2^si2); 

Uses 1st lie, si lOe, and sil lOe. 

If point space is not parabolic, the unit cycle UC is the reflection of real line in Cl Of and we draw both of them. 
36b (Cayley transform pictures 35c) += (32b) «36a 36c > 

if (si ^ 0 ) { 

ex-to<cyc\e2D>(UC.subs(cycle-val, subs-options::algebraic \ subs-options::no-pattern )) 

. asy-draw(asymptote, xmin, xmax, ymin, ymax, lst(0, 0, 0.7), " + 1.5", true, 7); 

C10f.subs(cycle-val , subs-options:: algebraic \ subs-options::no-pattern).normalize() 

.asy-draw(asymptote, xmin , xmax , ymin, ymax, lst(0, 0.7, 0), (si2 =sil ? "+1" : "+Dotted "), true, 7); 

Uses cycle2D 6g 11c 11c 34b 35b 39d 39d 39d 39d 39d 48b 48b 48b 48b 48b 56b 56b 57d, 1st lie, si lOe, and sil lOe. 

In the parabolic space unit cycle obtained from the real line by cayley-parabQ procedure. 

36c (Cayley transform pictures 35c) += (32b) <36b 36dc> 

> else 

exJ,o<cyc\e2T)>(cayley-parab(reaLline,signl).subs(cycle-val, subs-options::algebraic \ subs_options\:no-pattern)) 
.asy-draw(asymptote, xmin, xmax, ymin, ymax, lst(0, 0, 0.7), " + 1.5", true, 7); 

> 

Uses cycle2D 6g 11c 11c 34b 35b 39d 39d 39d 39d 39d 48b 48b 48b 48b 48b 56b 56b 57d and 1st lie. 

The pictures are finished with standard stuff. 

36d (Cayley transform pictures 35c) += (32b) «36c 

(Put units 3lc) 

(Draw axes 31b) 

asymptote <C "shipout(\"cayley-"<C eph_names[si+1 ] <C eph-names[sil+1 ]<C"\");" endl, 

} 

} 


Uses eph_names lOd, si lOe, and sil lOe. 

C.2.8. Three types of inversions. We draw here pictures for three types of the inversions. First we make a rectangular 
grid. 

36e (Three inversions 36e)= (32b) 37a c> 

xmin=- 2; xmax= 2; ymin=- 2; ymax= 2; 

C^=cycle2D (1st (0,(1- abs(sign))+2) ,e, 1); 

Gl?=cycle2D(0,lst(/,«),m,,e); 

asymptote <C endl <C "eraseO; u=lcm;" <C endl; 
for(double i=- 4; i< 4; i+=A ) { 

C3.subs(\st(sign=-1, 1=0, n= 1, m=i)).asy-draw( 

asymptote, xmin, xmax, ymin, ymax, lst(0.5, .75, 0.5),"+0.25pt", true, 7); 

C3.subs(\st(sign=-\, 1=1, n= 0, m=i)).asy-draw( 

asymptote, xmin, xmax, ymin, ymax, lst(0.5, .5, 0.75),"+0.25pt", true, 7); 

} 

C2.subs(sign=-l).asy-draw(asymptote, xmin, xmax, ymin, ymax, lst(l,0,0),"+.75pt", true, 7); 

(Draw axes 31b) 

asymptote <C "shipout(\"pre-invers\"); " <C endl; 


Uses cycle2D 6g 11c 11c 34b 35b 39d 39d 39d 39d 39d 48b 48b 48b 48b 48b 56b 56b 57d, 1st lie, and u 62b. 


Now we define inversions of the grid lines in the unit cycle and draw them for three different metrics. 

37a (Three inversions 36e)+= (32b) «36e 37b > 

C4= C3. cyclesimilarity( C2 ); 
for(int si=- 1; si< 2; si++) { 
asymptote <C endl <C " erase ();" <C endl; 
for(double i=- 4; i< 4; i+=.4) { 

C4-subs(\st(sign=si, 1=0, n= 1, m=i)).asy-draw( 

asymptote, xmin, xmax, ymin, ymax, lst(0.5, .75, 0.5),"+0.25pt", true, 9); 

C4-subs(\st(sign=si, 1=1, n= 0, m=i)).asy-draw{ 

asymptote, xmin, xmax, ymin, ymax, lst(0.5, .5, 0.75),"+0.25pt", true, 9); 

> 

C2.subs(sign=si).asy-draw(asymptote, xmin, xmax, ymin, ymax, lst(l,0,0),"+.75pt", true, 7); 

Uses 1st lie and si lOe. 

We conclude by drawing the image of the cycle at infinity Zinf 

37b (Three inversions 36e)+= (32b) «37a 

ex-to<cycle2D>(Zinf.cyclesimilarity(C2)).subs(sign=si).asy-draw( 
asymptote, xmin, xmax, ymin, ymax, 1st(0,0,1), (se-1? "+3pt": "+.75pt")); 

(Draw axes 31b) 

asymptote <C "shipout(\"inversion-" <C eph-names[si+l\ <C <C endl, 

> 

Uses cycle2D 6g 11c 11c 34b 35b 39d 39d 39d 39d 39d 48b 48b 48b 48b 48b 56b 56b 57d, eph_names lOd, 1st lie, and si lOe. 

Appendix D. The Implementation the Classes cycle and cycle2D 
This is the main file providing implementation the Classes cycle and cycle2D. It is not well documented yet. 

D.l. Cycle and cycle2D classes header files. 

D.1.1. Cycle header file. This the header file describing the classes cycle and cycle2d. We start from the general 
inclusions and definitions and then defining those two classes. 

37c (cycle.h 37c) = 

^include <stdexcept> 

^include <ostream> 

^include <sstream> 
using namespace std ; 

^include <ginac/ginac ,h> 
using namespace GiNaC; 

#define CYCLELIB_MAJ0R_VERSI0N 1 
#define CYCLELIB_MIN0R_VERSI0N 1 

(Auxiliary functions headers 38a) 

(cycle class 38b) 

(cycle2D class 39b) 


Defines: 

CYCLELIB_MAJ0R_VERSI0N, never used. 
CYCLELIB_MIN0R_VERSI0N, never used. 


D.1.2. Some auxiliary functions. Here is the list of some auxiliary functions which are defined and used in the cycle. h. 

38a (Auxiliary functions headers 38a) = (37c) 

/* * Check of equality of two expression and report the string */ 

DEC LA RE_FUN CTI ON_lP(jump.fnct) 

const char *equality( const ex & E); 

inline const char *equality( const ex & El, const ex & E2) { return equality(El-E2 );} 
inline const char *equality( const ex & E, const ex & solnsl, const ex & solns2 ) 

{ ex e = E; return equality(e.subs(solnsl), e.subs(solns2 ));} 

/* * Return the string describing the case (elliptic, parabolic or hyperbolic) */ 

const char *eph_case( const numeric & sign); 

/* * Return even (real) part of a Clifford number */ 

inline ex scalar_part{ const ex & e) { return remove-dirac-ONE(normal(canonicalize-clifford(e 
+ clifford-bar(e))))~numeric(2);y 

/* * Return odd part of a Clifford number */ 

inline ex clifford-part( const ex & e) { return normal(canonicalize-difford(e - c(i(fonL&ar( e)))-^numeric(2);} 

/* * Produces a Clifford matrix form of element of SL2 */ 

matrix sl2-clifford(c onst ex & a, const ex & b, const ex & c, const ex & d, const ex & e, bool not_inverse= true); 
matrix sl2-clifford(const ex & M, const ex & e, bool notjinverse= true); 


Defines: 

eph.case, used in chunk 20f. 
equality, used in chunk 13. 

jump_fnct, used in chunks lOe, 15, 18, 32a, 50b, 56c, 66, and 67. 

Uses bool lie, ex 4a lib 39d 48a 48b 66a 66b 66c 67a, matrix 8c lOe Ilf 12a 25a, and numeric lOd. 


D.1.3. Members and methods in class cycle. The class cycle is derived from class basic in GiNaC according to the 
general guidelines given in the GiNaC tutorial, is defined through the general s 

38b (cycle class 38b) = (37c) 

/* * The class holding cycles kx~2-2<l,x>+m=0 */ 

class cycle : public basic 

{ 

GINAC-DECLARE_REGISTERED_CLASS{ cycle, basic) 

(cycle class constructors 2a) 

(service functions for class cycle 39a) 

(accessing the data of a cycle 3a) 

(specific methods of the class cycle 4b) 

(Linear operation as cycle methods 3d) 

protected: 

ex unit; //A Clifford unit to store the dimensionality and metric of the point space 
ex k; 
ex l: 
ex m; 

>; 

(Linear operation on cycles 3e) 


Uses cycle 2a 2b 3e 3e 3e 3e 3e 46a 46b 47b 47b 48a 48a 48a 48a 48a and ex 4a lib 39d 48a 48b 66a 66b 66c 67a. 


39a 


This is a set of the service functions which is required that a cycle is properly archived or printed to a stream. 

(service functions for class cycle 39a) = (38b) 

// internal constructors 

//cycle(const ex & kl, const ex & 11, const ex & ml, const ex & metr, const exvector & v, bool discard¬ 
able = false); 

//cycle(const ex & kl, const ex & 11, const ex & ml, const ex & metr, std::auto_ptr<exvector> vp); 

protected: 

void do-print( const print-dflt & c, unsigned level ) const; 
void do-print-dflt( const print_dflt & c, unsigned level ) const; 
void do_printJate i(const printJatex & c, unsigned level ) const; 

Uses bool lie, cycle 2a 2b 3e 3e 3e 3e 3e 46a 46b 47b 47b 48a 48a 48a 48a 48a, ex 4a lib 39d 48a 48b 66a 66b 66c 67a, and v 62b. 

D.1.4. The derived class cycle2D for two dimensional cycles. We derive a derived class cycle2D from cycle in order 
to add some more methods which only make sense in two dimensions. 

39b (cycle2D class 39b) = (37c) 

class cycle2D : public cycle 

{ 

GINAC.DECLARE_REGISTERED_CLASS{cycle2B, cycle) 

(constructors of the class cycle2D 6f) 

(methods specific for class cycle2D 7b) 

(duplicated methods for class cycle2D 39c) 

>; 

(duplicated linear operation on cycle2D 39d) 


Uses cycle 2a 2b 3e 3e 3e 3e 3e 46a 46b 47b 47b 48a 48a 48a 48a 48a and cycle2D 6g 11c 11c 34b 35b 39d 39d 39d 39d 39d 48b 48b 48b 
48b 48b 56b 56b 57d. 

The general framework developed in the cycle class have some duplicates for two dimensions. 

39c (duplicated methods for class cycle2D 39c) = (39b) 

inline cycle2D sw6s(const ex & e, unsigned options = 0) const { 
return ex-to<cyc\e2T»{inherited::subs)e, options))', } 

inline cycle2D normalize) const ex & k-new = numeric(l), const ex & e = 0) const { 
return ex-to<cycle2T)>{inherited::normalize{k_new, e)); } 
inline cycle2D normalize-det) const ex & e = 0) const { 

return ex-to<cyc\e2D>{inherited::normalize-det(e)); } 

inline cycle2D normal () const { return cycle2D (k.normal)), l.normal )), m.normalQ, unit.normalQ );} 


Uses cycle2D 6g 11c 11c 34b 35b 39d 39d 39d 39d 39d 48b 48b 48b 48b 48b 56b 56b 57d, ex 4a lib 39d 48a 48b 66a 66b 66c 67a, 
and numeric lOd. 

We also specialise for the derived class cycle2D all operations defined in § 2.3 

39d (duplicated linear operation on cycle2D 39d)= (39b) 

const cycle2D operator-t- (const cycle2D & Ih , const cycle2D & rh)', 
const cycle2D operator-(const cycle2D & Ih, const cycle2D & rh); 
const cycle2D operator*(const cycle2D & Ih, const ex & rh)', 
const cycle2D operator*(const ex & Ih, const cycle2D & rh)', 
const cycle2D operators (const cycle2D & Ih, const ex & rh); 
const ex operator*(const cycle2D & Ih, const cycle2D & rh); 


Defines: 

cycle2D, used in chunks 6f, 7a, 12c, 14-18, 20a, 21g, 24a, 25b, 30a, 32, 33, 35-37, 39-41, 55-57, and 62-64. 
ex, used in chunks 2-8, 11a, 13, 14, 17, 19-25, 33b, 38-40, 42, 44-47, 49-60, 63, 65a, 67, and 68. 


D.2. Implementation of the cycle class. We start from definitions of constructors in cycle class 

40a (cycle.cpp 40a) = 40b > 

^include <cycle.h> 

#define PRINT.CYCLE c.s « \ 

k. print(c, level)', \ 

c.s < ", \ 

l. print(c, level); \ 

c.s < ", \ 

m. print(c, level); \ 
c.s < ")"; 

GINA CL IMPL EMEN 7 R KG IS TER ED_ CL A .S'S'_ 0 P T( cycle. basic, 
print-func<print-dflt> (&cycle:: do-print). 
print-func<printJatex> (&cycle:: do-print.-latex)) 

GINACJMPLEMENT_REGISTERED_CLASS(cyc\e2D, cycle) 

//, print_func<print_dflt>(&cycle2D::do_print) 

cycle::cycle() : unitQ, k(), l(), mQ 

{ 

tinfo-key = &cycl e:\tinfo_static; 

> 


Defines: 

PRINT.CYCLE, used in chunk 47b. 

Uses cycle 2a 2b 3e 3e 3e 3e 3e 46a 46b 47b 47b 48a 48a 48a 48a 48a and cycle2D 6g 11c 11c 34b 35b 39d 39d 39d 39d 39d 48b 48b 48b 
48b 48b 56b 56b 57d. 

D.2.1. Main constructor of cycle from all parameters given. If all parameters of the cycle are given this constructor 
is used. 

40b (cycle.cpp 40a) += < 1 40a 40c [> 

cycle::cycle(const ex & kl, const ex & 11, const ex & ml, const ex & metr) // Main constructor 
: k[kl), m(ml) 

{ 

ex D, metric; 

Uses cycle 2a 2b 3e 3e 3e 3e 3e 46a 46b 47b 47b 48a 48a 48a 48a 48a and ex 4a lib 39d 48a 48b 66a 66b 66c 67a. 

The first portion of the code processes various form of presentation for l. 

40c (cycle.Cpp 40a) += <1 40b 40dl> 

if (is-a<indexed>(ll.simplify-indexed())) { 
l = ex-to<indexed>(ll.simplify_indexed()); 
if ( ex-to<indexed>(l).get-indices().size() = 1) { 

D = ex-to<idx> (ex-to<indexed> (1). get-indices()[0]). get_dim(); 
y else 

thro-w(std::invalid-argument(" cycle: :cycle() : the second parameter should be an indexed object" 
"with one varindex")); 

} else if (is_a<matrix>(W) A ( min(ex-to<m&trix>(ll).rows(), ex-to<ma.trix>(ll).cols()) =1)) { 

D = max(ex-to<ma.trix>(11).rows)), ex-to<m&trix>(ll).cols()); 
l = indexed(Zl, varidx((new symbol)—> setflag(status-flags::dynallocated), D)); 

> else if (ll.info(info_flags::list) A (ll.nops() > 0)) { 

D = ll.nopsQ; 

l = indexed(matrix(l, ll.nopsQ, ex-to<lst>(ll)), varidx((new symbol)—> setflag(status-flags::dynallocated), D)); 

Uses cycle 2a 2b 3e 3e 3e 3e 3e 46a 46b 47b 47b 48a 48a 48a 48a 48a, 1st lie, matrix 8c lOe Ilf 12a 25a, and varidx lOd. 

If 11 is zero we will try to get missing information from the metrix in the next chunk, otherwise throw an exception. 
40d (cycle.Cpp 40a) += <40c 41a l> 

} else if (not ll.is-zeroQ) 

th.vow(std::invalid-argument(" cycle: : cycleO : the second parameter should be an indexed object, " 
"matrix or list")); 


Uses cycle 2a 2b 3e 3e 3e 3e 3e 46a 46b 47b 47b 48a 48a 48a 48a 48a and matrix 8c lOe Ilf 12a 25a. 


Now we process the metric parameter, in case 11 did not provide information on the dimensionality we try to get it 
here. 

41a (cycle.Cpp 40a) += «40d 41b > 

if (is_a<clifFord>(mefr)) { 
if (D .is-zero ()) 

D = ex-to<idx>(metr.op(l)).get_dim()] 

unit = cliff orcLwmt(varidx(0, D ), eic_fo<clifford>( rnetr).get-metric ()); 

} else { 

if ( D.is-zeroQ ) { 
if ( is-d< indexed> ( metr) ) 

D = ex-to<idx>(metr.op(l)).get_dim()] 
else if (is_a<matrix>(mefr)) 

D = ea;_to<matrix>(metr).rows(); 
else { 

exvector indices = metr.get-free-indices(); 
if (indices.sizeQ = 2) 

D = ex_to<idx>(md*ces[0]).geLdzm(); 

> 

} 


Uses matrix 8c lOe Ilf 12a 25a and varidx lOd. 

For metric of unknown type we throw an exception. 

41b (cycle.Cpp 40a) += <d 41a 41c > 

if (D .is-zero()) 

thro’w(std::invalid-argument(" cycle: :cycle() : the metric should be either tensor, " 

"matrix, Clifford unit or indexed by to indices. " 

"Otherwise supply the through the second parameter.")); 

unit = clifford-unit( varidx(0, D ), metr)] 

> 

Uses cycle 2a 2b 3e 3e 3e 3e 3e 46a 46b 47b 47b 48a 48a 48a 48a 48a, matrix 8c lOe Ilf 12a 25a, and varidx lOd. 

Now we come back to the case 11 is zero and try to resolve it with new info on D. 

41c (cycle.cpp 40a) += < 1 41b 42a > 

if ( ll.is-zeroQ ) { 
if ( not D.is-zero ()) 

l = indexed(0, varidx((new symbol)—>setflag(status-flags::dynallocated), £>)); 
else 

throw(std::invalid-argument("cycle :: cycle () : the second argument is zero and the metric " 
"does not tell the dimensionality of space")); 

> 

(Set tinfo to dimension 4ld) 

} 


Uses cycle 2a 2b 3e 3e 3e 3e 3e 46a 46b 47b 47b 48a 48a 48a 48a 48a and varidx lOd. 

We set tinfo key for cycle according to its dimension 

4ld (Set tinfo to dimension 41d)= (41c) 

if ( 'is_«<numoric> ( I))') 
switch (ex-to<numeric>(D).to_int()) { 
case 2: 

tinfo-key = &cycle2D ::tinfostatic; 

break; 

default: 

tinfo-key = &cycl ewtinfostatic, 

break; 

> 

else 

tinfo-key = &cycl ewtinfostatic] 


Uses cycle 2a 2b 3e 3e 3e 3e 3e 46a 46b 47b 47b 48a 48a 48a 48a 48a, cycle2D 6g 11c 11c 34b 35b 39d 39d 39d 39d 39d 48b 48b 48b 48b 
48b 56b 56b 57d, and numeric lOd. 


D.2.2. Specific cycle constructors. Constructor for cycle with the given determinant rsquared, e.g. zero-radius cycle 
by default. 

42a (cycle.Cpp 40a) += «41c 42b > 

cycle::cycle(const 1st & l, const ex & metr, const ex & r_squared, const ex & e, const ex & sign) 

{ 

symbol m-temp; 

cycle C(numeric(l), 1, m_temp, metr); 

( *this ) = C.subject_to(\st{C.det(e, sign) = rsquared), 1st (m-temp)); 

} 


Uses cycle 2a 2b 3e 3e 3e 3e 3e 46a 46b 47b 47b 48a 48a 48a 48a 48a, ex 4a lib 39d 48a 48b 66a 66b 66c 67a, 1 st lie, and numeric lOd. 
This is the constructor of a cycle identical to the given one with replaced metric in the point space. 

42b (cycle.Cpp 40a) += «42a 42c > 

cycle::cycle(const cycle & C, const ex & metr) 

{ 

(*this) = metr.is-zero{)l C : cycle (C.get-kQ, C.getJ (), C.get-m{), metr)', 

} 


Uses cycle 2a 2b 3e 3e 3e 3e 3e 46a 46b 47b 47b 48a 48a 48a 48a 48a and ex 4a lib 39d 48a 48b 66a 66b 66c 67a. 

Constructor of a cycle from a matrix representations. First we check that matrix is in a proper form. 

42c (cycle.cpp 40a) += <42b 42d[> 

cycle::cycle(const matrix & M, const ex & metr, const ex & e, const ex & sign) 

{ 

if (not ( M.rows{) = 2 A M.colsf) = 2 A ( M.op(0)+M.op(3)).normal().is-zeroQ )) 

th.vow(std\:invalid-argument("cycle: :cycle() : the second argument should be square 2x2 matrix " 
"with M(1,1)=-M(2,2)")); 

(Create a Clifford unit 44a) 

varidx i0((ne w symbol)—> setflag(status-flags::dynallocated), D), 
il(( new symbol)— >setflag(status-flags::dynallocated), D, true); 

Uses cycle 2a 2b 3e 3e 3e 3e 3e 46a 46b 47b 47b 48a 48a 48a 48a 48a, ex 4a lib 39d 48a 48b 66a 66b 66c 67a, matrix 8c lOe Ilf 12a 25a, 
and varidx lOd. 

There are different options for sign , which should be checked. First we verify is it zero and use the default value in 
this case. 

42d (cycle.cpp 40a) += <42c 43a > 

if ( sign.is-zero ()) { 

try { 

(*this) = cycle(remove-dirac-ONE(M.op(2)), clifford-toJst(M.op(0), el), remove-dirac-ONE(M.op( 1)), metr)', 

} catch (std::exception Sxp) { 

(*this) = cycle(numeric(l), clifford-toJst(M.op(0)*clifford-inverse(M.op(2)), el), 
canonicalize-difford(M.op(l)*clifford-inverse(M.op(2))), metr); 

} 

> else { 

ex sign_m, conv; 
sigmm = sign.evalmQ; 

Uses catch 10b 25d 26b 26c 43b, cycle 2a 2b 3e 3e 3e 3e 3e 46a 46b 47b 47b 48a 48a 48a 48a 48a, ex 4a lib 39d 48a 48b 66a 66b 66c 67a, 
and numeric lOd. 


< 42d 43b > 


If sign is not zero we process different types which can supply it. 

(cycle.cpp 40a) += 

if ( is-a<tensor>(sign-m)) 

conv = indexed(ex_fo<tensor>(s*gn_m), iO, il): 
else if (is_a<clifford>(sign_m)) { 
if ( ex-to<idx>(sign-m.op(l)).get-dim () = D ) 
conv = ex-to<c\i&ord>{sign-m).get-metric(iO, il)] 
else 

thro'w(std::invalid-argument(" cycle: :cycle() : the sign should be a Clifford unit with " 

"the dimensionality matching to the second parameter")); 

} else if (is_a< indexed >(si<7n_m)) { 
exvector ind = ea;_to<indexed> (signjm).get-indices))', 

if ))ind.size)) = 2) A ( ex-to<\dx>)ind[d\). get-dim )) = D) A (eaLfo<idx>(md[l]).<?eLdim() = D)) 
conv = sign-m.subs(lst(ind[0\ = iO, ind[ 1] = il))', 

else 

t\iY<yw( y std\'.invalid-argument{"cycle :: cycle () : the sign should be an indexed object with two " 
"indices and their dimensionality matching to the second parameter")); 

} else if (is_a<matrix> (sign-m)) { 

if (( ex-to<ma.trix>(sign_m).cols() = D) A ( ex-to<m&trix>(sign-m).rows () = D )) 
conv = indexed(ea;_to<matrix>(s*(?n_m), iO, il)', 
else 

thro'w(std::invalid-argument("cycle:: cycle() : the sign should be a square matrix with the " 
"dimensionality matching to the second parameter")); 

} else 

throw(std:'.invalid_argument("cycle: :cycle() : the sign should be either tensor, indexed, matrix 
"or Clifford unit")); 

Uses cycle 2a 2b 3e 3e 3e 3e 3e 46a 46b 47b 47b 48a 48a 48a 48a 48a, 1st lie, and matrix 8c lOe Ilf 12a 25a. 

Then all blocks of the matrix are used to construct the cycle in main constructor. 

(cycle.cpp 40a) += <3 43a 44b > 

try { 

( *this ) = cycle(remove-dirac-ONE(M.op(2)), indexed(matrix(l, ex_fo<numeric> (D).toJnt(), 

clifford-toJst(M.op(0), el)), i0.toggle-variance())*conv, remove-dirac-ONE(M.op(l)), metr)', 

} catch (std::exception &p) { 

( *this) = cycle(numeric(l), indexed(matrix(l, ex-to<numeric>(D).toJnt(), clifford-toJst(M.op)O) 

* clifford_inverse(M.op(2)), el)), iO.toggle-variance)))* conv, 
canonicalize-clifford(M.op)l)*clifford-inverse{M.op{2))), metr)', 

} 

} 

} 


Defines: 

catch, used in chunks 42d and 50a. 

Uses cycle 2a 2b 3e 3e 3e 3e 3e 46a 46b 47b 47b 48a 48a 48a 48a 48a, matrix 8c lOe Ilf 12a 25a, and numeric lOd. 


We need the proper Clifford unit to decompose M(0,0) element into vector for l. 

44a (Create a Clifford unit 44a) = (42c) 

ex el, D\ 
if ( e.is-zero ()) { 
ex metrl ; 

if (is_a<matrix>(metr)) { 

D = ex-to<matrix>(metr).cols (); 
metrl = metr, 

} else if (is_a<clifford> (metr)) { 

D = ex-to<idx>(metr.op(l)).get_dim (); 
metrl = ex-to<c\i&ord>(metr).get-metric (); 

} else if (is_a<indexed>(metr)) { 

D = ex-to<idx.>(ex-to<indexed>(metr).getj,ndices()[0]).get-dim (); 
metrl = metr, 

> else 

throw(,std:: invalid-argument) "Could not determine the dimensionality of point space 
"from the supplied metric or Clifford unit")); 

el = c/ij[foreLumf(varidx((new symbol)— >setflag(status-flags::dynallocated), D), metrl)', 

} else { 
el = e; 

D = ex-to<idx>(e.op(l)).get-dim(); 

} 


Uses ex 4a lib 39d 48a 48b 66a 66b 66c 67a, matrix 8c lOe Ilf 12a 25a, and varidx lOd. 

D.2.3. Class cycle members access. Class cycle has four operands. 

44b (cycle.Cpp 40a) += «43b 45a > 

ex cycle: :op(size-t i ) const 
{ 

GINA C-A SSER T( i< nops ()); 

switch (i) { 
case 0: 
return k; 
case 1: 
return l; 
case 2: 
return m; 
case 3: 
return unit', 
default: 

throw (std'.'.invalidLargument^" cycle: : opO : requested operand out of the range (4)")); 

> 

} 


Uses cycle 2a 2b 3e 3e 3e 3e 3e 46a 46b 47b 47b 48a 48a 48a 48a 48a and ex 4a lib 39d 48a 48b 66a 66b 66c 67a. 


< 44b 45b > 


Operands may be set through this method. 

45a (cycle.cpp 40a) += 

ex & cyc\e::let-op{sizeJ i) 

{ 

GINA C-A SSER T( i< nops ()); 

ensureJfjmodifiablef); 

switch (i) { 
case 0: 
return k; 
case 1: 
return l; 
case 2: 
return m; 
case 3: 
return unit, 
default: 

thro'w(std::invalid-argument(" cycle: : let_op() : requested operand out of the range (4)")); 

> 

> 


Uses cycle 2a 2b 3e 3e 3e 3e 3e 46a 46b 47b 47b 48a 48a 48a 48a 48a and ex 4a lib 39d 48a 48b 66a 66b 66c 67a. 
Substitutions works as usual in GiNaC. 

45b (cycle.cpp 40a) += «45a 45c i> 

cycle cycle::su&s(const ex & e, unsigned options) const 

{ 

exmap m; 

if ( e.info(info-flags::list)) { 

1st l = ea;_to<lst>(e); 

for (1st: :const_iterator i = l.beginQ ; i ^ l.end()', ++i) 
m.insert(std::make-pair(i-^op(0 ), *—>op(l))); 

} else if (*s_a<relational>(e)) { 
m.insert(std::makejpair{e.op(0 ), e. op( 1))); 

> else 

t\iY(yw(std:'.invalid-argument{"cycle: :subs() : the parameter should be a relational or a 1st")); 

return ex-to<cyc\e>(inherited::subs(m, options)); 

} 


Uses cycle 2a 2b 3e 3e 3e 3e 3e 46a 46b 47b 47b 48a 48a 48a 48a 48a, ex 4a lib 39d 48a 48b 66a 66b 66c 67a, 1st lie, and relational lie. 

D.2.4. Service methods for the GiNaC infrastructure. Standard parts involving archiving, comparison and printing of 
the cycle class 

45c (cycle.Cpp 40a) += «45b 46a l> 

cycle::cycle(const archive-node &n, 1st &.symJst) : inherited(n, symJst) 

{ 

n./mcLex("k-param", k, symJst); 
n./mcLex("l-param", l, symJst); 
n./mcLex("m-param", m, symJst); 
n.find-ex{" unit", unit, symJst); 

} 


Uses cycle 2a 2b 3e 3e 3e 3e 3e 46a 46b 47b 47b 48a 48a 48a 48a 48a and 1st lie. 


Archiving routine. 

46a (cycle.cpp 40a) += «45c 46b > 

void cycle: :archive(archive-node &n) const 

{ 

inherited:: archive(n ); 
n.adc?_ea;("k-param", k); 
n.add-ex(" 1-param", t); 
n.add-ex(" m-param", m); 
n.add_ex{" unit", unit); 

} 

ex cycle: :unarchive(const archive-node &n, 1st SksymJst) 

{ 

return (new cycle(n, symJst))^setflag(status-flags::dynallocated); 

} 


Defines: 

cycle, used in chunks 2-9, 11-17, 20e, 24, 25, 33a, 38-47, 49-56, and 58. 

Uses ex 4a lib 39d 48a 48b 66a 66b 66c 67a and 1st lie. 

Comparison of cycles. 

46b (cycle.cpp 40a) += «46a 46c > 

int cycle: : compares ame-type(const basic Mother) const 

{ 

const cycle &o = static_cast<const cycle &>(ot/ier); 

if ((unit = o.unit) A (l*o.get-k() = o.getJ()*k) A (m*o.get-k() = o.get-m()*k)) 

return 0; 

else if ((unit < o.unit) 

V (l*o.get-k() < o.getJ()*k) V (m*o.get-k() < o.get-m()*k)) 

return -1; 
else 

return 1; 

> 


Defines: 

cycle, used in chunks 2-9, 11-17, 20e, 24, 25, 33a, 38-47, 49—56, and 58. 

Equality of cycles. 

46c (cycle.cpp 40a) += <46b 47a [> 

bool cycle: :is-equal(const basic & other) const 

{ 

if (not is-a<cyc\e>(other)) 

return false; 

const cycle o = ex-to<cyc\e>(other); 

if (-1 (unit.is-equal(o.unit) A (m*o.get-k()-o.get-m()*k).normal().is-zero())) 

return false; 

if ( is-a<numeric>(get-dim())) { 
int D = ex-to<numeric> (get-dimQ) .to-int(); 
for (int i= 0; i<D; ?++) 

if (-i (get-l(i)*o.get-k()-o.get-l(i)*k).normal().is-zero()) 

return false; 
return true; 

} else 

return (ko.get-k()).normal)).is-equal((o.getJ()*k),normal()); 

} 


Uses bool lie, cycle 2a 2b 3e 3e 3e 3e 3e 46a 46b 47b 47b 48a 48a 48a 48a 48a, and numeric lOd. 


< 1 46c 47b > 


A cycle is zero if and only if its all components are zero 

47a (cycle.cpp 40a) += 

bool cycle:: is-zero() const 

{ 

return ( k.is-zero( ) A l.is-ZeroQ A m.is-zero {)); 

> 


Uses bool lie and cycle 2a 2b 3e 3e 3e 3e 3e 46a 46b 47b 47b 48a 48a 48a 48a 48a. 

Printing of cycles. 

47b (cycle.Cpp 40a) += <47a 47c > 

void cycle: :do-print(const print-dflt & c, unsigned level ) const 

{ 

PRINT.CYCLE 

> 

void cycle: : do_printJatex(const printJatex & c, unsigned level) const 

{ 

PRINT.CYCLE 

} 


Defines: 

cycle, used in chunks 2-9, 11-17, 20e, 24, 25, 33a, 38-47, 49-56, and 58. 

Uses PRINT.CYCLE 40a. 

D.2.5. Linear operation on cycles. Here are linear operations on cycle defined as methods. 
47c (cycle.cpp 40a) += <47b 48a > 

cycle cycle: :add(const cycle & rh ) const 

{ 

return cycl e(get-k()+rh.get-k(), getJ()+rh.getJ(getJQ.op(l)), 
get-m()+rh.get-m (), getjmetric{))\ 

> 

cycle cycle: :su&(const cycle & rh) const 

{ 

return cycl e(get-k()-rh.get-k(), getJ()-rh.getJ(getJ(J.op(l)), 
get-m()-rh.get-m(), getjmetric{))\ 

} 

cycle cycle:: exmid(const ex & rh) const 

{ 

return cycl e(get-k()*rh, indexed((getJ() .op(0)*rh),evalm (), getJ().op(l)), 
get-m()*rh , get-metric())', 

> 

cycle cycle: :dw(const ex & rh) const 

{ 

return exmul(pow(rh, numeric(-l))); 

> 


Uses cycle 2a 2b 3e 3e 3e 3e 3e 46a 46b 47b 47b 48a 48a 48a 48a 48a, ex 4a lib 39d 48a 48b 66a 66b 66c 67a, and numeric lOd. 


The same linear structure is represented in operators overloading. 

48a (cycle.cpp 40a) += «47c 48b > 

const cycle operator+(const cycle & Ih, const cycle & rh) 

{ 

return lh.add(rh ); 

} 

const cycle operator-(const cycle & Ih, const cycle & rh) 

{ 

return lh.sub(rh ); 

> 

const cycle operator* (const cycle & Ih , const ex & rh) 

{ 

return Ih.exmul(rh)-, 

} 

const cycle operator* (const ex & Ih, const cycle & rh) 

{ 

return rh.exmul(lh); 

} 

const cycle operator-^(const cycle & Ih, const ex & rh) 

{ 

return Ih.div(rh)- 

} 

const ex operator* (const cycle & Ih, const cycle & rh) 

{ 

return Ih.mul(rh)-, 

} 

Defines: 

cycle, used in chunks 2-9, 11-17, 20e, 24, 25, 33a, 38-47, 49—56, and 58. 

ex, used in chunks 2-8, 11a, 13, 14, 17, 19-25, 33b, 38-40, 42, 44-47, 49-60, 63, 65a, 67, and 68. 
We make a specialisation of these operation for cycle2D class as well. 

48b (cycle.cpp 40a) += «48a 49a > 

const cycle2D operator-!-(const cycle2D & Ih, const cycle2D & rh) 

{ 

ret urn exJ,o< cycle2D > ( Ih. add{ rh )); 

} 

const cycle2D operator-(const cycle2D & Ih, const cycle2D & rh) 

{ 

return ex_to<cycle2D>( Ih.sub(rh))- 

} 

const cycle2D operator* (const cycle2D & Ih, const ex & rh) 

{ 

return ea;_fo<cycle2D> (Ih.exmul(rh )); 

} 

const cycle2D operator* (const ex & Ih, const cycle2D & rh) 

{ 

ret urn exJ,o< cycle2D > ( rh. exmul( Ih))-, 

} 

const cycle2D operator-^ (const cycle2D & Ih, const ex & rh) 

{ 

return ea;_to<cycle2D> (Ih.div(rh )); 

> 

const ex operator* (const cycle2D & Ih, const cycle2D & rh) 

{ 

ret urn exJ,o< cycle2D > ( Ih. mul{rh))\ 

} 


Defines: 

cycle2D, used in chunks 6f, 7a, 12c, 14-18, 20a, 21g, 24a, 25b, 30a, 32, 33, 35-37, 39-41, 55-57, and 62-64. 
ex, used in chunks 2-8, 11a, 13, 14, 17, 19-25, 33b, 38-40, 42, 44-47, 49-60, 63, 65a, 67, and 68. 


D.2.6. Specific methods for cycle. 


We oftenly need to normalise cycles to get rid of ambiguity in their definition. This is typically by prescribing a 
value to k. 

49a (cycle.Cpp 40a) += <48b 49b > 

cycle cycle: :normalize(const ex & k_new, const ex & e) const 

{ 

ex ratio = 0; 

if ( k_new.is-zero( )) // Make the determinant equal 1 
ratio = sqrt(det(e))\ 

else { // First non-zero coefficient among k, m, LO, LI, ... is set to k_new 
if (-> k.is-zero ()) 
ratio = k-h-k-new] 
else if (-i m.is-zero()) 
ratio = m-^-k-new, 
else { 

int D = ex-to<numeric>(get-dim()).to-intQ ; 
for (int i= 0; i<D ; *++) 
if (-il.subs(l.op( 1) = i).is-zero()) { 
ratio = l.subs(l.op( 1) = i)+k-new; 
i = D; 

> 

} 

} 

if (ratio.is-zero()) // No normalisation is possible 
return ( *this)\ 

if (ratio.is-zero()) 
return (*this)', 

return cycl e(k-^-ratio, indexed((l.op(0)-i-ratio).evalm().normal(), l.op( 1)), (m-f ratio).normal(), get-metric{))\ 

> 


Uses cycle 2a 2b 3e 3e 3e 3e 3e 46a 46b 47b 47b 48a 48a 48a 48a 48a, ex 4a lib 39d 48a 48b 66a 66b 66c 67a, and numeric lOd. 
The normalisation to unit determinant 

49b (cycle.cpp 40a) += <49a 50a > 

cycle cycle: :normalize-det(const ex & e) const 

{ 

ex d = det(e)', 

return (d.is-zero()7 *this: normalize(k-^-sqrt(d), e)); 

} 


Uses cycle 2a 2b 3e 3e 3e 3e 3e 46a 46b 47b 47b 48a 48a 48a 48a 48a and ex 4a lib 39d 48a 48b 66a 66b 66c 67a. 


This methods returns a centre of the cycle depending from the provided metric. 

50a (cycle.cpp 40a) += «49b 50b > 

ex cycle: : centeffconst ex & metr, bool return-matrix) const 
{ 

if ( is-a<rmmeric>(get-dim ())) { 
ex el, D = get-dim {); 
if ( metr.is-zero( )) 
el = unit ; 

else { 

if (?s_a<clifford>(mefr)) 

el = clifford-unit(v&ridx.(Q, D), ex-to<cliSord>(metr).get-metric()); 
else 
try { 

el = clifford-unit( varidx(0, D), metr)', 

} catch ( exception &p) { 

throw(std::invalid-argument(''cycle : : center () : supplied metric" 

" is not suitable for Clifford unit")); 

> 

> 


50b 


Uses bool lie, catch 10b 25d 26b 26c 43b, cycle 2a 2b 3e 3e 3e 3e 3e 46a 46b 47b 47b 48a 48a 48a 48a 48a, 
ex 4a lib 39d 48a 48b 66a 66b 66c 67a, numeric lOd, and varidx lOd. 

Finally, the centre is constructed for the cycle and given metric by the formula [10, Defn. 2.10]: 
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(cycle.Cpp 40a)+= < 50a 51a[> 

1st c; 

for(int ?=0; i<D; i++) 
if ( k.is-zero()) 
c.append(get-l(i))-, 

else 

//c.append(jump_fnct(-ex_to<clifford>(el).get_metric(varidx(i, D), varidx(i, D)))*get_l(i)/k); 
c.append(-ex-to<clif£ord>(el).get-metric(varidx.(i, D), varidx(j, D))*get-l(i)^rk)\ 
return ( return-matrix ? (ex)matrix(er_to<numeric>(I?).to_mf(), 1, c) : (ex)c); 

} else { 
return l+k; 

} 

} 


Uses ex 4a lib 39d 48a 48b 66a 66b 66c 67a, jump_fnct 38a, 1st lie, matrix 8c lOe Ilf 12a 25a, numeric lOd, and varidx lOd. 



D.2.7. Build cycle with given properties. We oftenly need cycles with prescribed properties, e.g. when converting of 
cycles to normalised form or matrix. This routine takes a system of linear equations with the cycle parameters and 
try to resolve it. The list of unknown parameters is either supplied or build automatically in a way suitable for most 
applications. 

51a (cycle.Cpp 40a) += «50b 51b > 

cycle cycle: :subject-to(const ex & condition, const ex & vars) const 
{ 

1st varsl ; 

if ( vars.info(info-flags::list ) A (vars.nopsQ ^ 0)) 
varsl = ex-to<\st>(vars)] 
else if (?s_o<symbol>(i>ars)) 
varsl = lst( vars) ; 

else if ((vars = 0) V ( vars.nops( ) = 0)) { 
if ( *s_a<symbol>(?7i)) 
varsl. append( m ); 
if ( is-a< numerio ( get-dim( ))) 

for (int i = 0; i < e:r_fo<numeric> (get-dim()).to-double()', i++) 
if (zs_a<symbol>(< 7 eiJ(i))) 
vars 1 .append(getJ(i))-, 
if (zs_a<symbol>(A;)) 
vars 1. append ( k ); 
if (varsl.nops() = 0) 

throw(std::invalid-argument("cycle: : subj ect_to () : could not construct the default list of " 
"parameters")); 

> else 

throw (std::invalid-argument(" cycle: : subject_to() : second parameter should be a list of symbols" 

" or a single symbol")); 

return subs(lsolve(condition.info(info-flags\:relation-equat)l 1st (condition) : condition, 
varsl), subs-options::algebraic \ subs-optionsv.no-pattern)', 

> 


Uses cycle 2a 2b 3e 3e 3e 3e 3e 46a 46b 47b 47b 48a 48a 48a 48a 48a, ex 4a lib 39d 48a 48b 66a 66b 66c 67a, 1st lie, and numeric lOd. 

D.2.8. Conversion of the cycle to the matrix form. This method is inverse to the constructor of the cycle from its 
matrix, see (2.2) and [10, § 3.1]. First, we process the supplied e to the standard form of the Clifford unit. 

51b (cycle.cpp 40a) += < 51a 52a > 

matrix cycle::to-matrix(const ex & e, const ex & sign) const 
{ 

ex one, es, conv, D = get-dim(); 

varidx id((new symbol )-^>setflag(status-flags::dynallocated), D), 
il(( new symbol)— > setflag(status-flags::dynallocated), D, true); 
if ( e.is-zero ()) { 
one = dirac-ONE (); 

es = cliff'ord-unit(il.toggle-variance(), get-metric ()); 

} else if (is.a< clifford>(e)) { 

one = dirac-ONE(ex-to<c\i&ord>(e).get-representationJabel()); 
es = e.subs(e.op( 1) = il.toggle-variance ()); 

} else if (zs_a<tensor>(e) V *s_a<indexed>(e) V is-a< matrix>(e)) { 
one = dirac-ONE()- 

es = clifford-unit(il.toggle-variance(), e); 

} else 

throw(std::invalid-argument(" cycle: :to_matrix() : expect a Clifford number, matrix, tensor or " 
"indexed as the first parameter")); 

Uses cycle 2a 2b 3e 3e 3e 3e 3e 46a 46b 47b 47b 48a 48a 48a 48a 48a, ex 4a lib 39d 48a 48b 66a 66b 66c 67a, matrix 8c lOe Ilf 12a 25a, 
and varidx lOd. 


<1 51b 52b > 


Then we work out the sign, which should be used. 

(cycle.cpp 40a) += 

ex signjrn = sign.evalmQ; 

if (is_a<tensor > {signjrn)) 

conv = indexed(ex_to<tensor>(sign_m), iO, il ); 
else if (*s_a<clifford>(sign_m)) { 
if (ex-to<idx>(sigri-m.op(l)).get-dim() = D) 
conv = exJ,o<c\if£ord>(signjm).get_metric{iO, il)', 
else 

t\iv(yw(std::invalid-argument(" cycle : :to_matrix() : the sign should be a Clifford unit with the " 
"dimensionality matching to the second parameter")); 

} else if (is_a<indexed>(si^n_m)) { 
exvector ind = ex-to<indexed>(sign-m).get-indices()', 

if (( ind.size() = 2) A (ex-to<idx>(ind[0)).get_dim() = D) A (ex-to<idx>(ind[l]).get-dim() = D)) 
conv = sign _m. sub s(lst(md[0] = iO , ind[ 1] = il))', 
else 

thro'w(std::invalid-argunient("cycle: :to_matrix() : the sign should be an indexed object with two 
"indices and their dimensionality matching to the second parameter")); 

} else if (is_a<matrix>(sign_m)) { 

if ((ex_to<matrix> (sign_m).cols(J = D) A ( ex-to<ma.trix>(sign-m).rows() = D)) 
conv = indexed(ea;_to<matrix> (sign-m), iO, il ); 

else 

t\iv(yw{std::invalid-argument("cycle: :to_matrix() : the sign should be a square matrix with the " 
"dimensionality matching to the second parameter")); 

} else 

thro'w(std::invalid-argunient(" cycle: :to_matrix() : the sign should be either tensor, indexed, " 
"matrix or Clifford unit")); 

Uses cycle 2a 2b 3e 3e 3e 3e 3e 46a 46b 47b 47b 48a 48a 48a 48a 48a, ex 4a lib 39d 48a 48b 66a 66b 66c 67a, 1st lie, 
and matrix 8c lOe Ilf 12a 25a. 

When all components are ready the matrix is build in few lines. 

(cycle.cpp 40a) += < 52a 53a o 

ex aOO = expand-dummysum(l.subs(ex-to<indexed>(l).get_indices()[0\ = i0.toggle-variance()) * conv* es); 

return matrix(2, 2, 1st [aOO, m * one, k * one, -aOO)); 

} 


Uses ex 4a lib 39d 48a 48b 66a 66b 66c 67a, 1st lie, and matrix 8c lOe Ilf 12a 25a. 


D.2.9. Calculation of a value of cycle at a point. This is used in the construction of a relational cycle '/.passing 
describing incidence of a point to cycle. 

53a (cycle.Cpp 40a) += «52b 53b > 

ex cycle::ua/(const ex & y ) const 

{ 

ex yO , D = get-dim(); 

varidx iO, il; 
if (is_a<indexed>(y)) { 

iO = ex-to< varidx>(e:r_io< indexed >(y).getJndices()[0]); 
if ((ex-to<indexed>(y).getJndices().size() = 1) A ( i0.get_dim () = D)) { 
yO = ex-to< indexed >{y); 

il = varidx((new symbol)—> setflag(status-flags::dynallocated), D); 

} else 

thx<yw (std::invalid-argument{" cycle: :val() : the second parameter should be a indexed object with " 
"one varindex")); 

} else if ( y.info(info-flags::list ) A (y.nopsQ = D)) { 

10 = varidx((new symbol)— >setflag(status-flags::dynallocated), D ); 

11 = varidx((new symho\)^setflag{status-flags::dynallocated) 1 D); 
yO = indexed(matrix(l, y.nops {), ex-to<lst>(y )), i0)\ 

} else if (is_a<matrix>( 2 /) A ( mm{exJ,o< matrix>( y).rrrws (). ex-to<m.&trix>(y).cols()) =1) 

A (D = max(ex-to<ma.trix>(y).rows(), ex-to<m.atrix>(y).cols()))) { 

10 = varidx((new symbol )—>setflag(status-flags::dynallocated) 1 D ); 

11 = varidx((new symbol )^>setflag(status-flags::dynallocated), D); 
yO = indexed (y, iO ); 

} else 

thi'cyw(std::invalid-argument{"cycle: :val() : the second parameter should be a indexed object, " 
"matrix or list")); 

return expand-dummy_sum(-k*y0*y0.subs(i0 = il)*get_metric(i0.toggle-variance (), il.toggle-variancef)) 

-2* I*y0.subs(i0 = ex-to<\raridx>(ex-to<indexed>(l).get-indices()[0]).toggle-variance ()) +m); 

} 


Uses cycle 2a 2b 3e 3e 3e 3e 3e 46a 46b 47b 47b 48a 48a 48a 48a 48a, ex 4a lib 39d 48a 48b 66a 66b 66c 67a, 1st lie, 
matrix 8c lOe Ilf 12a 25a, and varidx lOd. 

D.2.10. Matrix methods for cycle. The method det() may be defined in several ways. An alternative to the present 
definition is pseudodeterminant [2, (4.9)] 

ex cycle::det(const ex & e = 0, const ex & sign = (new tensdelta)^setflag(status-flags::dynallocated )) const 
{ex M = normalize^).to-matrix(e, sign); 

return remove-dirac-ONE(M.op(0)*clifford-star(M.op(3))-M.op(l)*cliffordstar(M.op(2))) ; > 

However due to the structure of matrix this coincides with the usual determinant of the matrix. 

53b (cycle.cpp 40a) += «53a 54a > 

ex cycle: :def(const ex & e, const ex & sign , const ex & kjnorm) const 
{ 

return remove-dirac-ONE({k_norm.is-zero() r !*this:normalize(k_norm)).to-matrix(e, sign).determinant^)); 

y 


Uses cycle 2a 2b 3e 3e 3e 3e 3e 46a 46b 47b 47b 48a 48a 48a 48a 48a and ex 4a lib 39d 48a 48b 66a 66b 66c 67a. 


Multiplication of cycles in the matrix representations and their similarity with respect to elements of SL 2 (WL) and 
other cycles. 

54a (cycle.Cpp 40a) += «53b 54b > 

matrix cycle::mid(const ex & C, const ex & e, const ex & sign, const ex & signl) const 
{ 

if (is_a<cycle>((7)) 

return ea;_to<matrix>( canonicalize-clifford)to-matrix) e, sign).mul)ex-to< cycle>( C).to-matrix)e, 
signl. is-zero))? sign: signl))))', 

else if (zs_o<matrix>((7) A (e:r_fo<matrix>( C).rowsQ = 2) A (ea;_to<matrix>( C).cols)) = 2)) 
return ex-to< matrix>( canonicalize-clifford)to-matrix)e, sign).mul(ex-to< matrix>( C) ))); 
else 

th.row(std::invalid-argument("cycle: :mul() : cannot multiply a cycle by anything but a cycle " 

"or 2x2 matrix")); 

> 


Uses cycle 2a 2b 3e 3e 3e 3e 3e 46a 46b 47b 47b 48a 48a 48a 48a 48a, ex 4a lib 39d 48a 48b 66a 66b 66c 67a, and matrix 8c lOe Ilf 12a 25a. 

D.2.11. Actions of cycle as matrix, cycle in the matrix form can act on other objects, or matrices can acts on cycle. 
Any 2 x 2-matrix acts on a cycle by the similarity: M : C i—> MCM~ X . 

54b (cycle.cpp 40a) += «54a 54c > 

cycle cycle::matrixsimilarity)const ex & M, const ex & e, const ex & sign, bool not-inverse) const 
{ 

if ( not (is_a<matrix>(M) A ex-to<m&trix> )M) ,rows))=2 A ex-to<m&trix> )M). cols))=2)) 

thro'w)std::invalid-argument)" cycle: :matrix_similarity () : the first parameter sgould be a 2x2 matrix")) 

return matrixsimilarity(M.op{ 0), M.op( 1), M.op(2), M.op(3), e, sign, not_inverse); 

> 


Uses bool lie, cycle 2a 2b 3e 3e 3e 3e 3e 46a 46b 47b 47b 48a 48a 48a 48a 48a, ex 4a lib 39d 48a 48b 66a 66b 66c 67a, 
and matrix 8c lOe Ilf 12a 25a. 

The same method works if the matrix is provided by its four elements. 

54c (cycle.Cpp 40a) += <54b 54dl> 

cycle cyc\e::m,atrixsimilarity(const ex & a, const ex & b, const ex & c, const ex & d, const ex & e, 
const ex & sign, bool not_inverse) const 

{ 

return cycle(ex-to<m.a.trix>(canonicalize-difford(ma.trix(2,2,not-inverse7lst(a, b, c, d) : 1st( cliffordstar(d), - 
cliffordstar(b), - cliffordstar(c), cliffordstar(a))) 

.mul(mul(m&trix(2,2,not-inverse7\st(clifford-star(d), - cliffordstar(b), - cliff’ordstar(c), clif- 
ford_star(a)):lst(a, b, c, d)), e, sign)) 

.evalmQ).normal))), get-metric )), e, sign)', 

} 


Uses bool lie, cycle 2a 2b 3e 3e 3e 3e 3e 46a 46b 47b 47b 48a 48a 48a 48a 48a, ex 4a lib 39d 48a 48b 66a 66b 66c 67a, 1st lie, 
and matrix 8c lOe Ilf 12a 25a. 

For elements of SX 2 (R) we have a specific method which make the proper “cliffordization” of the matrix first. 

54d (cycle.cpp 40a) += «54c 55a > 

cycle cycle: :sl2similarity )const ex & a, const ex & b, const ex & c, const ex & d, const ex & e, const ex & sign, 
bool not-inverse) const 

{ 

return cycle)ex-to<matrix>)canonicalize-clifford) 

sl2-clifford)a, b, c, d, e.is-zero))7unite, not-inverse) 

.mul)mul)sl2-difford)a, b, c, d, e.is-zero))7unite, -> not-inverse ), e, sign)) 

.evalm)).subs)c*b = )d*a- 1), subs-options::algebraic \ subs-options::no-pattern)).normal))), get-metric)), e, sign)', 

} 


Uses bool lie, cycle 2a 2b 3e 3e 3e 3e 3e 46a 46b 47b 47b 48a 48a 48a 48a 48a, ex 4a lib 39d 48a 48b 66a 66b 66c 67a, 
and matrix 8c lOe Ilf 12a 25a. 


55a 


(cycle.cpp 40a) += < 1 54d 55b > 

cycle cycle: :sl2similarity (const ex & M, const ex & e, const ex & sign, bool notjinverse) const 

{ 

if ( is_a<matrix>(M) V M.info)info_flags::list )) 

return sl2similarity)M.op) 0), M.op) 1), M.op) 2), M.op) 3), e, sign, not_inverse); 

else 

throw(std::OT^al*d_ar<7wne7it("sl2_cliffordO : expect a list or matrix as the first parameter")); 

} 


Uses bool lie, cycle 2a 2b 3e 3e 3e 3e 3e 46a 46b 47b 47b 48a 48a 48a 48a 48a, ex 4a lib 39d 48a 48b 66a 66b 66c 67a, 
and matrix 8c lOe Ilf 12a 25a. 

cycle acts on other cycle by the similarity: C : C\ i—> CCiC, see [10, (4.8)]. 

55b (cycle.cpp 40a) += «55a 55c > 

cycle cycle:: cycle-similarity (const cycle & C, const ex & e, const ex & sign , const ex & signl) const 
{ 

return cyc\e)ex-to<ma.trix>)canonicalize-clifford)C.mul)mul)C, e, sign, signl.is-zeroQ?sign:signl), e, 
sign l.is-zero))? sign: signl))), get-metric)), e, sign); 

} 


Uses cycle 2a 2b 3e 3e 3e 3e 3e 46a 46b 47b 47b 48a 48a 48a 48a 48a, ex 4a lib 39d 48a 48b 66a 66b 66c 67a, and matrix 8c lOe Ilf 12a 25a. 

55c (cycle.cpp 40a) += <55b 55di> 

relational cycle::iss-orthogonal)const cycle & C, const ex & e, const ex & sign) const 
{ 

return )C. cycle-similarity)* this, e, sign). get-l)get-dim))-l).normal)) = 0); 

> 

Uses cycle 2a 2b 3e 3e 3e 3e 3e 46a 46b 47b 47b 48a 48a 48a 48a 48a, ex 4a lib 39d 48a 48b 66a 66b 66c 67a, and relational lie. 

D.3. Implementation of the cycle2D class. The derived class cycle2D for two dimensional cycles. Here con¬ 
structors, archiving, and comparison come first. 

55d (cycle.cpp 40a) += «55c 55ec> 

cycle2D::cycle2D() : inherited () 

{ 

tinfoJzey = &cycle2D wtinfostatic; 

> 

Uses cycle2D 6g 11c 11c 34b 35b 39d 39d 39d 39d 39d 48b 48b 48b 48b 48b 56b 56b 57d. 

55e (cycle.cpp 40a) += <55d 55fo 

cycle2D::cycle2D(const ex & kl, const ex & 11, const ex & ml, const ex & metr) 

: inherited)kl, 11, ml, metr) 

{ 

if )get-dim{) ^ 2) 

thro'w)std::invalid-argument)" cycle2D: : cycle2D() : class cycle2D is defined in two dimensions")); 
tinfo-key = &cycle2D ::tinfo_static; 

} 

Uses cycle2D 6g 11c 11c 34b 35b 39d 39d 39d 39d 39d 48b 48b 48b 48b 48b 56b 56b 57d and ex 4a lib 39d 48a 48b 66a 66b 66c 67a. 

55f (cycle.cpp 40a) += <55e 56a > 

cycle2D::cycle2D(const 1st & l, const ex & rsquared, const ex & metr, const ex & e, const ex & sign) 

: inherited(l, rsquared, metr, e, sign) 

{ 

if ( get-dim)) ^ 2) 

thro'w)std::invalid-argument)" cycle2D: : cycle2D() : class cycle2D is defined in two dimensions")); 
tinfo-key = &cycle2D wtinfostatic; 

> 

Uses cycle2D 6g 11c 11c 34b 35b 39d 39d 39d 39d 39d 48b 48b 48b 48b 48b 56b 56b 57d, ex 4a lib 39d 48a 48b 66a 66b 66c 67a, 
and 1st lie. 


< 55f 56b > 


56a (cycle.cpp 40a) += 

cycle2D::cycle2D(const cycle & C , const ex & metr) 

{ 

(*this ) = cycle2D (C.get-k(), C.getJ(), C.get-m (), ( metr.is-zero()? C.get-metricQ: metr)); 

} 

Uses cycle 2a 2b 3e 3e 3e 3e 3e 46a 46b 47b 47b 48a 48a 48a 48a 48a, cycle2D 6g 11c 11c 34b 35b 39d 39d 39d 39d 39d 48b 48b 48b 48b 
48b 56b 56b 57d, and ex 4a lib 39d 48a 48b 66a 66b 66c 67a. 

56b (cycle.cpp 40a) += «56a 56c > 

void cyc\e2T)::archive(archive-node &n) const 

{ 

inherited:: archive(n ); 

} 

cycle2D::cycle2D(const archive-node &n, 1st &csymJst) : inherited(n, symJst) {; } 

ex cycle2D::unarchive(const archive-node &n, 1st Sxsym-lst) 

{ 

return (new cycle2D(n, symJst))^setflag(status-flagsr.dynallocated); 

} 

int cycle2D::comparesame-type(const basic & other) const 

{ 

return inherited:: comparesame-type(other); 

> 


Defines: 

cycle2D, used in chunks 6f, 7a, 12c, 14-18, 20a, 21g, 24a, 25b, 30a, 32, 33, 35-37, 39-41, 55-57, and 62-64. 

Uses ex 4a lib 39d 48a 48b 66a 66b 66c 67a and 1st lie. 

D.3.1. The member functions of the derived class cycle2D. The standard definition of the focus for a parabola is 

/1 m l 2 n\ 

2n 2 nk + 2k) 

We calculate focus of a cycle based on its determinant in the corresponding metric. 

56c (cycle.cpp 40a) += <56b 56di> 

ex cycle2D::/ocus(const ex & e, bool return-matrix) const 

{ 

1st f=\st{jump-fnct{-get-metric(ysLV\Ax.(Q, 2), varidx(0, 2)))*getJ(0)-hfc, 

(-det(e, (new tensdelta)—>setflag(status-flags::dynallocated), k)-^-(2*getJ(l)*k)).normal^)); 
return ( return-matrix ? (ex)matrix(2, 1 , f) : (ex)/); 

} 

Uses bool lie, cycle2D 6g 11c 11c 34b 35b 39d 39d 39d 39d 39d 48b 48b 48b 48b 48b 56b 56b 57d, ex 4a lib 39d 48a 48b 66a 66b 66c 67a, 
jump_fnct 38a, 1st lie, matrix 8c lOe Ilf 12a 25a, and varidx lOd. 

56d (cycle.cpp 40a) += «56c 57a > 

1st cycle2D::roots(const ex & y, bool first) const 

{ 

ex D = get-dim(); 

1st ksign = lst(-fct=geLmetnc(varidx(0, D ), varidx(0, D)), -fc*geLmetnc(varidx(l, D), varidx(l, I?))); 
int iO = (first? 0:1), il = (first? 1:0); 

ex c = ksign.op(il)*pow(y , 2) - numeric(2)*getJ(il)*y+m; 
if ( ksign.op(iO).is-zero()) 

return (getJ(iO).is-zero() ? lst() : lst(c-^<7etJ(z0)-^numeric(2))); 
else { 

ex disc = sqrt(pow(get-l(i0) , 2) - ksign.op(i0)*c); 

return 1st ((get-l(i0)-disc)-^ksign.op(i0), (getJ(i0)+disc)-i-ksign.op(i0)); 

} 

> 

Uses bool lie, cycle2D 6g 11c 11c 34b 35b 39d 39d 39d 39d 39d 48b 48b 48b 48b 48b 56b 56b 57d, ex 4a lib 39d 48a 48b 66a 66b 66c 67a, 
1st lie, numeric lOd, and varidx lOd. 


«56d 57b > 


57a (cycle.Cpp 40a) += 

1st cycle2T)::line-intersect(const ex & a, const ex & b) const 

{ 

ex D = get-dim()', 

ex pm = -fc*geLmetnc(varidx(l, D), varidx(l, D))\ 
return cycle2D(k*(numeric(l)+pm*pow(a,2)).normal(), lst((getJ(0)+getJ(l)*a-pm*a*b).normal (), 0), {m- 
2*getJ(l)*b-\-pm*pow( b, 2)). normally). rootsQ ; 

} ' 

Uses cycle2D 6g 11c 11c 34b 35b 39d 39d 39d 39d 39d 48b 48b 48b 48b 48b 56b 56b 57d, ex 4a lib 39d 48a 48b 66a 66b 66c 67a, 1st lie, 
numeric lOd, and varidx lOd. 


D.3.2. Drawing cycle2D. Some auxiliary functions used for drawing 

57b (cycle.cpp 40a) += «57a 57c > 

inline ex max( const ex ha, const ex hb) {return ex-to<numeric>((a-b).evalj{)).is-positive()'!a:b',y 
inline ex m*n( const ex ha, const ex hb) {return ex-to<mimeric>((arb).evalfO).is-positive()?b:a;y 


Uses ex 4a lib 39d 48a 48b 66a 66b 66c 67a and numeric lOd. 

The most complicated member function in the class cycle2D 

57c (cycle.Cpp 40a) += «57b 57d[> 

^define PAIR(X, Y) ex_to<numeric> ( (X) . evalf () ) . to_double () << << \ 

ex-to< numerio (( Y). evalf {)). to.doubleQ 

^define DRAW_ARC(X, S) u = ex_to<numeric> ( (X) . evalf ()). to_double 0 ; \ 

v= ex_to<numeric>(roots(X, -<not-swapped).op(zero-or-one).evalfO).to-double()', \ 
du = dir*(-k-d* signv*v+lv)', \ 
dv = dir*(k_d*signu*u-lu ); \ 

if ( notswapped) \ 

ost<t: £< ")*u{" < \ 

else \ 

ost <C S' <C v <C "," <C u <C ") *u{" <C (sign = 0? dv : -dv) <C " , " <C (sign = 0? du : -du) <C 
Defines: 

DRAW_ARC, used in chunk 65b. 

PAIR, used in chunks 60-62 and 64b. 

Uses du 62b, dv 62b, k_d 62b, numeric lOd, u 62b, v 62b, and zero_or_one 62b. 

The main drawing routine for cycle2D. 

57d (cycle.cpp 40a) += «57c 58a o 

void cycle2D ::metapost_draw(ostream & ost , const ex & xmin , const ex & xmax , const ex & ymin , const ex & ymax , 
const 1st & color, const char * more-options, bool with_header, 

int points-per_arc, bool asymptote, char * picture, bool only-path, bool is-continuation) const 
{ 

exD= get-dim()', 

ostringstream drawstart, draw-options', 
char* already-drawn', 

if ( is-continuation ) already-drawn ="““ (" ; 

else already-drawn = " ("; // Was any arc already drawn? 

drawstart <C "draw" <C (asymptote ? "(" : " ") -C pictured (int(*picture)=t)l "" : -C 

ios-base::fmtflags keep-flags = ost.flagsQ; // Keep stream’s flags to be restored on the exit 
draw-options.flags(keep-flags)', // Synchronise flags between the streams 
draw-options.precision(ost.precision()); // Synchronise flags between the streams 

Defines: 

cycle2D, used in chunks 6f, 7a, 12c, 14-18, 20a, 21g, 24a, 25b, 30a, 32, 33, 35-37, 39-41, 55-57, and 62-64. 

Uses bool lie, ex 4a lib 39d 48a 48b 66a 66b 66c 67a, and 1st lie. 


Each drawing command is concluded by options containing color, etc. They are formatted differently for Asymptote 
and MetaPost. 

58a (cycle.cpp 40a) += «57d 58b > 

ost < fixed; 
draw-options <C fixed; 
if (color.nopsQ = 3) 
if ( asymptote ) 
draw-options <C ",rgb(" 

<C ex-to<numeric>(color.op(0)).to-double() <C 
<C ex-to<nuineric>(color.op(l)).to-double() -C" , " 

<C ex-to<numeric>(color.op(2)).to-double() <C 
else 

draw-options <C showpos -C " withcolor " 

<C ex-to<numeric>(eolor.op(0)).to-double() <C "*red" 

<C ex-to<numeric>(color.op(l)).to-double() <C"*green" 


<C ex-to<numeric>(color.op(2)).to-double() <C "*blue 
draw-options <C more-options <C (asymptote ? ;" : ";<C endl; 


Uses numeric lOd. 

A drawing command can be also preceded by a human-readable comment describing the cycle to be drawn. 

58b (cycle.cpp 40a) += «58a 58c > 

if ( with-header) 

ost <C (asymptote ? "// Asymptote" : Metapost") <C " data in [" -C xmin <C 

•C xmax -C "] x [" -C ymin -C " , " 

<C ymax <C "] for " -C (ex)passm^(lst(symbol("u"), symbol("v"))); 

if (k .is-zero() A l.subs(l.op(\) = 0 ).is-zero() A l.subs(l.op( 1) = 1 ).is-zero() A m.is-zeroQ) { 

ost <K " zero cycle, (whole plane) " <C endl; 
ost.flags(keep-flags); 

return; 

> 


Uses cycle 2a 2b 3e 3e 3e 3e 3e 46a 46b 47b 47b 48a 48a 48a 48a 48a, ex 4a lib 39d 48a 48b 66a 66b 66c 67a, 1st lie, u 62b, and v 62b. 
There are several parameters which control the output. Their values depend from either we draw cycle in the original 
coordinates or swap the u and v 

58c (cycle.Cpp 40a) += «58b 59a > 

ex xc = center().op( 0), yc = center)).op( 1); // the center of cycle 

double signO = ea;_to<numeric>(-gef_mefr*c(varidx(0, D), varidx(0, D)).evalf()).to-double (), 
signl = ea;_to<numeric>(-^eLmetnc(varidx(l, D), varidx(l, D)).evalfifi).to-double(), 
sign = signO * signl; 

numeric determinant = ea;_to<nuineric> (det().evalf()); 

bool notswapped = (sign>0 V signl= 0 V ((sign <0) A ~^determinant.is-positive())); 

double signu = (notswapped! signO: signl), signv = (notswapped?signl:signO); 

int iu = (notswapped? 0:1), iv = (notswapped? 1:0); 

ex umin = (notswapped ? xmin : ymin), umax = (notswapped ? xmax : ymax), 
vmin = (notswapped 1 ymin: xmin), vmax = (notswapped ? ymax : xmax), 
uc = (notswapped 1 xc: yc), vc = (notswapped ? yc : xc); 

1st bsoots = roots(vmin, notswapped), tsoots = roots(vmax, notswapped); 

Uses bool lie, cycle 2a 2b 3e 3e 3e 3e 3e 46a 46b 47b 47b 48a 48a 48a 48a 48a, ex 4a lib 39d 48a 48b 66a 66b 66c 67a, 1st lie, 
numeric lOd, and varidx lOd. 


Here is the outline of the rest of the method. It effectively splits into several cases depending from the space metric 
and degeneracy of cycle2D. 

59a (cycle.Cpp 40a) += <58c 66ac> 

(Draw a straight line 59b) 

(Find intersection points with the boundary 60b) 
if ( sign > 0) { // elliptic metric 
(Draw a circle 60c) 

> else { // parabolic or hyperbolic metric 
(Draw a parabola or hyperbola 62b) 

} 

ost -C endl; 
ost.flags(keep-flags ); 

return; 

} 

If line is detected we identify its visible portion. 

59b (Draw a straight line 59b) = (59a) 59c c> 

if ( b-roots.nops() ^ 2) { // a linear object 
if ( with_header ) 

ost -C " (straight line)" <C endl; 
ex ul, u2, vl, v2; 

if ( b-roots.nops() = 1){ // a ”non-horisontal” line 
ul = max{min)b-roots.op(0), umax), umin)\ 
u2 — min(max(t-roots.op(0) : umin ), umax); 

} else { // a ’’horisontal” line 
ul = umin; 
u2 = umax; 

} 


Uses ex 4a lib 39d 48a 48b 66a 66b 66c 67a. 

Vertical lines case. 

59c (Draw a straight line 59b) += (59a) <59b 59dc> 

if (getJ(iv).is-zeroQ) { // a vertical line 

if ( ex-to< numeric>(( bjroots. op( 0)- umin).evalf()).is-positive() 

A ex-to<numeric>((umax-b-roots.op(0)).evalf()).is-positive ()) { 
vl = vmin; 
v2 = vmax; 

} else { // out of scope 
ost. flags) keep-flags); 

return; 

> 


Uses numeric lOd. 

Look for the visible portion of generic line. 

59d (Draw a straight line 59b) += (59a) <59c 60a o 

> else { 

vl = roots)ul, -inotswapped).op)0); 
v2 = roots(u2, ^notswapped).op(0); 
if (( max{vl , v2) > vmax) V ( min)vl , v2) < vmin)) { 
ost.flags(keep-flags); 
return; //out of scope 
> 

> 


Actual drawing of the line. 
60a (Draw a straight line 59b) += 


(59a) < 59d 



<C ")*u—(" -C PAIR{notswapped ? u2: v2, notswapped ? v2: u2 ) <C ")*u" 

<C ( only-pathl "" : draw-options.strQ); 
already_drawn= " ~ “ ("; 
if ( with_header) 
ost <tC endl; 
ost.flags(keep-flags ); 
return; 

} 

Uses PAIR 57c and u 62b. 

Make initially this intervals (left[i], right[i]) irrelevant for drawing by default, if necessary, it will be redefined letter 


on. 

60b (Find intersection points with the boundary 60b) = 
ex left[ 2] = { max(min(uc, umax ), umin), 
max(min(uc, umax), umin )}, 


(59a) 



if ( ea;_<o<numeric>( b_roots. op(0). evalf)) ). is-real()) { 
if ( ex-to<numeric>((b_roots.op(0)-b-roots.op(l)).evalf()).is-positive ()) 
b_roots = lst(6_roo£s.op(l), b_roots.op{ 0)); // rearrange to have minimum value first 
left{0] = min(max(b_roots.op(0), umin), umax); 
right[ 0] = max(min(b_roots.op( 1), umax), umin); 

> ' 

if ( ex-to< numeric> (t~roots. op(0). evalfQ).is-real()) { 
if ( e:r_fo<numeric>((Lroo£s. op(0)-t-roots. op(l)). evalf))). is-positive()) 
t_roots = lst(Lroo£s.op(l), t_roots.op(0)); // rearrange to have minimum value first 
left[l] = min(max(t_roots.op(0), umin), umax); 
right) 1] = max(min(t-roots.op(l), umax), umin); 

} ' 

Uses ex 4a lib 39d 48a 48b 66a 66b 66c 67a, 1st lie, and numeric lOd. 

We start from the most involved case of a circle with a positive radius. Two this end we calculate coordinates w[2][4] 
and v[2][4] of endpoints for up to four arcs making the circle. The x-components of intersection points with vertical 
boundaries are rearranged appropriately. 

60c (Draw a circle 60c) = (59a) 60do 

if ( determinant.is-positiveQ) { 
exr= sqrt(det()), w[2][4], f[2][4]; 
if ( with_header) 

ost <!C " /circle of radius " C r < endl; 

Uses ex 4a lib 39d 48a 48b 66a 66b 66c 67a, u 62b, and v 62b. 

Depending from the y-position of the centre we draw different arcs. 

60d (Draw a circle 60c)+= (59a) «60c 61a c> 


if ( ex-to< numeric >((vc-umax).evalf()).is-positive ()) { 
w[0][2] = left) 1]; m[0][3] = right) 1]; 
u[l][2\ = left) 0]; m[1][3] = right) 0]; 
m[ 0][0] = it[l][0] = uc; 
w[0][l] = m[1][1] = uc; 


Uses numeric lOd and u 62b. 


(59a) <60d 61b > 


61a (Draw a circle 60c) += 

} else if (ex-to<numeric>((vc-vmiri) .evalf()) .is.positiveQ) { 
m[0][0] = left[ 1]; w[0][l] = right[ 1]; u[0][2] = right\ 0]; m[0][3] = left[0\- 
if ( ex.to< numeric>( ( uc-r-umin).evalfO).is.positivel )) 
u[l] [0] = u[l][3] = uc-r, 

else 

m[ 1][0] = u[l][3] = umin; 

if ( ex.to< numerio (( umax-uc- r) . evalfl )) . is.positivel )) 
u[l][l] = «[1][2] = uc+r, 
else 

u[l][l] = u[l][ 2 ] = umax\ 

Uses numeric lOd and u 62b. 

61b (Draw a circle 60c)+= (59a) «61a 61c > 

> else { 

u[ 0 ][ 0 ] = left[ 1 ]; «[ 0 ][ 1 ] = right[ 1 ]; 
m[ 1 ][ 0 ] = left{ 0 ]; u[l][l] = right[ 0 ]; 
w[0][2] = m[1][2] = mc; 
m[ 0][3] = u[l][3] = uc; 

> 


Uses u 62b. 

We calculate now the y-components of the endpoints corresponding to x-components found before.. 

6lc (Draw a circle 60c) += (59a) <6lb 6ld > 

1st y.roots', 

for (int j= 0; j< 2; j++) 
for (int i= 0; i< 4; z++) 

if ((u[j\[i]-uc).evalf().is-zero()) // Touch the horizontal boundary? 
v[j\[i] = (i= 0 V i =1? vc+r : vc-r ); 

else if ((u[j][i]-uc-r).evalf().is.zero{) V (u[j\[i]-uc+r).evalf().is-zero()) // Touch the vertical boundary? 
v[j\[t\ = vc\ 

else { 

y.roots = roots(u[j][i ], false); 

if ( ex_fo<numeric>( yjroots.op(0) . evalfQ).is.realQ) { // does circle intersect the boundary? 
if (i<2) 

v[j\[i] = min{max{y-roots.op(0), yjroots.op{ 1 )), vmax ); 

else 

x[)][j] = max(min{y.roots.opl 0 ), yjroots.op{ 1 )), vmin ); 

} else 
v[j\ [i] = vc ; 

> 


Uses 1st lie, numeric lOd, u 62b, and v 62b. 

Now we drawing up to four arcs which make the visible part of the circle. Each arc is defined through its two endpoints 
and tangent vector in them. 

6ld (Draw a circle 60c) += (59a) «6lc 62a > 

for (int v= 0 ; i< 4; *++) {// actual drawing of four arcs 
int s = (i= 0 V i =2? -1:1); 

if ((«[ 0 ][j] 7 ^ «[!][*]) v (v[ 0 ] [z] 7 ^ x[l][*])) //do not draw empty arc 

osf<C " " <C [only.pathl already-drawn: draw.start, stri)) <C PAIR{u[0][i], r[0][*]) <C ")*u{" -C PAIR(s*(v[0] [*]- 

vc), s*(itc-u[ 0 ][*])) 

-C (asymptote ? "}: : : "}. . . {") 

// << (asymptote ? ”}.. tension atleast 1.5 : ”>...{”) // temporary fix 

<C PAIR(s*(v[l][i]-vc), s*(mc-m[1][z])) -C "X" <C PAIR(u[ l][i], r[l][i]) <C ")*u" 

•C ( only.path ? "" : draw-options. str())', 
already.drawn =" ~ ~ ( 11 ; 

} 


Uses PAIR 57c, u 62b, and v 62b. 


Finally, for zero-radius circles we draw a point and do not draw anything for circles with an imaginary radius. 

62a (Draw a circle 60c)+= (59a) «6ld 

> else if ( det().is-zero ()) { 
if ( with-header ) 

ost -C " /circle of zero-radius" <C endl; 

ost <C (asymptote ? ( only-path ? already-drawn : "dot (") : "draw " ) <C picture <C (int(*p*cfure)= 0 ? "" : ",") 

<C (only-path ? "" : "(") 

<C PAIR(uc, vc) < ")*u" < ( only-path ? : draw-options. str())-, 

already-drawn= 

} else 

if ( with-header) 

ost -C " /circle of imaginary radius—not drawing" -C endl] 

Uses PAIR 57c and u 62b. 

First we look if the parabola or hyperbola are degenerates into two lines, then treat two types of cycles separately. 
62b (Draw a parabola or hyperbola 62b) = (59a) 

double u , v, du, dv, k-d = ex-to<numeric>(k.evalf()).to-double(), 
lu = ex-to<nuineric>{getJ(iu).evalf()).to-double(), 

Iv = ex-to<num.eric>(get-l(iv).evalf(J).to-double()] 

bool change-branch = ( sign ^ 0 ); // either to do a swap of branches 

int zero-or-one = ( sign e 0 V k-d*signv > 0 ? 0 : 1); // for parabola and positive k take first 

if (sign = 0 ) { 

(Treating a parabola 62c) 

} else { 

(Treating a hyperbola 64c) 

> 


Defines: 

du, used in chunk 57c. 

dv, used in chunk 57c. 

k_d, used in chunks 57c, 63b, and 65a. 

u, used in chunks 10, 11, 15d, 16d, 19—25, 30—36, 57c, 58b, 60-62, and 64b. 

v, used in chunks 10, 11, 15d, 16d, 19—23, 30, 32a, 34a, 39a, 57c, 58b, 60, and 61. 
zero_or_one, used in chunks 57c and 65. 

Uses bool lie and numeric lOd. 

For parabolas degenerated into two parallel lines we draw them by the recursive call of this function 
cycle2D:: metapost-drawQ. 

62c (Treating a parabola 62c) = (62b) 63a i> 

if (signO = 0 A getJ(0).is-zero()) { 
if ( with-header) 

ost <!C " /parabola degenerated into two horizontal lines" <C endl ; 
cycle2D(0, lst(0,1), 2*b-roots.op(0 ), get-metricQ).metapost-draw(ost , xmin 7 xmax , ymin , ymax , color , more-options , 
false, 0 , asymptote , picture , only-path , is-continuation ); 

cycle2D(0, lst(0,1), 2*b-roots.op(l), get-metricQ).metapost-draw(ost, xmin , xmax , ymin , ymax , color , more-options , 
false, 0, asymptote , picture , only-path 7 true); 
if ( with-header) 
ost <C endl ; 
ost.flags(keep-flags)] 
return; 


Uses cycle2D 6g 11c 11c 34b 35b 39d 39d 39d 39d 39d 48b 48b 48b 48b 48b 56b 56b 57d and 1st lie. 


Two vertical lines are drawn here 


63a (Treating a parabola 62c) += (62b) «62c 63b > 

} else if (signl = 0 A getJ(l).is-Zero()) { 
if ( with_header ) 

ost -C " /parabola degenerated into two vertical lines" <C endl\ 
cycle2D(0, lst(l, 0), 2*bjroots.op(0), get-metricQ) .metapost_draw(ost, xmin, xmax, ymin, ymax , color, more-options , 
false, 0, asymptote, picture , only-path, is-continuation ); 

cycle2D(0,1st(1, 0), 2*6_roots.op(l), get-metricQ) .metapost_draw(ost , xmin, xmax , ymin, ymax, color, more-options, 
false, 0, asymptote, picture, only_path, true); 
if ( with_header ) 
ost -C end/; 
ost.flags(keep-flags); 
return; 

} 

Uses cycle2D 6g 11c 11c 34b 35b 39d 39d 39d 39d 39d 48b 48b 48b 48b 48b 56b 56b 57d and 1st lie. 

If a proper parabola is detected we rearrange intervals appropriately in order to draw pieces properly. 

63b (Treating a parabola 62c) += (62b) «63a 63c > 

if ( withJieader) 
ost < " /parabola" <C endl; 

if ( ex.to<immeric>((right[0]-left[0]).evalJ{)).is.positive() A ex-to<numeric>((right[l\-left[l)).evalf()).is-positive()) 
if ( k-d*signu > 0) { //rearrange intervals 
ex e = left[ 1]; left[ 1] = right[0 ]; right[0 ] = left[0 ]; (e/t[0] =e; 

> else { 

ex e = left[ 1]; Ze/t[l] = ng/it[l]; right[l] = ng , /i<[0]; ng/itfO] =e; 

> 


Uses ex 4a lib 39d 48a 48b 66a 66b 66c 67a, k_d 62b, and numeric lOd. 

Parabolas can be exactly represented by a cubic Bezier arc if the second and third control points correspondingly are: 
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63c (Treating a parabola 62c) += 
for (int i =0; i < 2; z++) { 

if ( ex-to<rmm.eric>((right[i]-left[i]).evalf{)).is-positive()) { // a proper branch of a parabola 


ex cp[8]; 

if ( not_swapped ) { 
cp[0] = left[i]; 

cp[ 1] = vaZ(lst(cp[0],0))-r2-^getJ(l); 
cp[ 6] = n<?/it[z]; 

cp[7] = vaZ(lst(cp[6],0))-r2-^getJ(l); 

cp[2] = numeric(2,3)*cp[0]+numeric(l,3)*cp[6]; 

cp[3] = (numeric(l,6)*cp[0]*cp[0]*fc + numeric(l,3)*cp[0]*cp[6]*A; 

- numeric(2,3)*cp[0]*<7etJ(0)- numeric(l,3)*<7etJ(0)*cp[6]+?rH-2)-^<7etJ(l); 
cp[4] = numeric(l,3)*cp[0]+numeric(2,3)*cp[6]; 
cp[ 5] = (numeric(l,3)*cp[0]*fc*cp[6]-numeric(l,3)*cp[0]*getJ(0) 
-numeric(2,3)*<?etJ(0)*cp[6]+numeric(l,6)*fc*cp[6]*cp[6]+77H-2)-h<7efJ(l); 


Uses ex 4a lib 39d 48a 48b 66a 66b 66c 67a, 1st lie, and numeric lOd. 


The similar formulae for swapped drawing. 

64a (Treating a parabola 62c) += (62b) «63c 64b > 

} else { 

cp[ 1 ] = left[i\; 

cp[0] = raZ(lst(0,cp[l]))-^2-^getJ(0); 
cp[ 7] = right[i]; 

cp[ 6 ] = val(lst(0,cp[7]))+2+getJ(0)] 
cp[3] = numeric(2,3)*cp[l]+numeric(l,3)*cp[7]; 
cp[ 2] = (numeric(l, 6 )*cp[l]*cp[l]*fc + numeric(l,3)*cp[l]*cp[7]*fc 
- numeric( 2 , 3 )*cp[l]*< 7 etJ(l)- numeric(l, 3 )*< 7 etJ(l)*cp[ 7 ]+?rH- 2 )- 4 -getJ( 0 ); 
cp[ 5] = numeric(l,3)*cp[l]+numeric(2,3)*cp[7]; 
cp[ 4] = (numeric(l,3)*cp[l]*fc*cp[7]-numeric(l,3)*cp[l]*getJ(l) 
-numeric( 2 , 3 )*<?etJ(l)*cp[ 7 ]+numeric(l, 6 )*fc*cp[ 7 ]*cp[ 7 ]+? 7 H- 2 )-y< 7 efJ( 0 ); 

} 


Uses 1st lie and numeric lOd. 

The actual drawing of the parabola arcs. 

64b (Treating a parabola 62c) += (62b) «64a 

ost -C ( only-path ? already-drawn : drawstart.strQ) -C PAIR(cp[ 0], cp[l]) -C ")*u . . controls ("; 
if ( asymptote ) 

ost <C PAIR(cp[2], cp[3]) < ")*u and (" <C PAIR(cp[ 4], cp[5]) <C ")*u . . ("; 
else 

ost <C "(" -C PAIR(cp[2], cp[3]) -C ")*u) and ((" <C PAIR(cp[4\, cp[ 5]) <C ")*u) . . ("; 
ost <C PAIR(cp[6 ], cp[7]) < ")*u" <C (only-path ? "" : draw-options.str{)); 
already-drawn= " ~ “ ("; 

} 

} 


Uses PAIR 57c and u 62b. 

If a hyperbola degenerates into a light cone we draw it as two separate lines. 

64c (Treating a hyperbola 64c)= (62b) 65a c> 

if ( determinant.is-zeroQ ) { 
if ( with-header ) 

ost -C " / a light cone at (" <C xc <tC "," <C yc <C")" -C endl ; 
cycle2D(0, lst(l, 1), 2 *(uc+vc), get-metric()).metapost-draw(ost, xmin, xmax, ymin, ymax, color, more-options, 
false, 0, asymptote, picture, only-path, is-continuation)] 
cycle2D(0, lst(l, -1), 2*(itc-uc), get-metricQ) .metapost-draw(ost, xmin, xmax, ymin, ymax, color, more-options, 
false, 0, asymptote, picture, only-path, true); 
if ( with-header ) 
ost < endl, 
ost.flags(keep-flags)] 
return; 


Uses cycle2D 6g 11c 11c 34b 35b 39d 39d 39d 39d 39d 48b 48b 48b 48b 48b 56b 56b 57d and 1st lie. 


Otherwise we rearrange the interwals for hyperbola branches. 

65a (Treating a hyperbola 64c) += (62b) «64c 65b > 

} else { 
if ( with-header ) 

ost " /hyperbola" <C endl; 

if ( ex-to< numeric > (( vmin- vc ). evalf ()). is-positive ()) { 
ex e = left[ 1]; left[ 1] = right[ 0]; right[0] = left[ 0]; left[ 0] =e; 
change-branch = false; 
zero-or-one = ( k-d*signv > 0 ? 1 : 0); 

> 

if ( ex-to< numeric >((vc-vmax ) .evalfQ).is-positive ()) { 
ex c= left[ 1]; Ze//[1] = right[ 1]; right[ 1] = ng/if[0]; n^/it[0] =e; 
change-branch = false; 
zero-or-one = [k-d*signv > 0 ? 0 : 1); 

> 

> 


Uses ex 4a lib 39d 48a 48b 66a 66b 66c 67a, k_d 62b, numeric lOd, and zero_or_one 62b. 

Two arcs of the hyperbola are drown now 

65b (Treating a hyperbola 64c) += (62b) «65a 

int points = ( points-per-arc = 0? 5 : points-per-arc)', 
for (int i =0; i < 2; «++) { 

double dir = ex-to<numeric>(csgn(signv*(2*zero-or-one-1))).to-double(); //direction of the tangent vectors 
// double dir = ((sign == 0? lv : signv*(2*zero_or_one-l))<0?-l:l); direction of the tangent vectors (sec¬ 
ond alternative) 

if ( ex-to<rmm.eric>((right[i]-left[i]).evalf()).is-positive()) { // a proper branch of the hyperbola 
DRAW-ARC(left[i], ( only-path ? already-drawn : draw-start.str()))', 
for (int j= 1; j<points ; j++) { 

DRAW-ARC(left[i]*(1.0-j+(points-1.0))+right[i]*j+(points-1.0), (asymptote ? " : : (" : "...(") ); 

} 

ost <C (only-path ? "" : draw-options.strQ); 
already-drawn= " ~ ~ ("; 

> 

if ( change-branch ) 

zero-or-one = 1 - zero-or-one\ // make a swap for the next branch of hyperbola 

} 


Defines: 

points, used in chunks lOd, 11a, 24b, 30a, and 32a. 
Uses DRAW_ARC 57c, numeric lOd, and zero_or_one 62b. 


D.4. Auxiliary functions implementation. The auxiliary functions defined as well. 


D.4.1. Heaviside function. We define Heaviside function: x( x ) = 1 for x > 0 and xi x ) = 0 for x < 0. 

66a (cycle.cpp 40a) += «59a 66b > 

////////// 

// Jump function 

////////// 

static ex jump-fnct-evalf(cor\st ex & arg) 

{ 

if ( is-exactly-a<numeric>(arg )) { 

if {(ex-to<num.eric>(arg).is-real() A ex-to<numeric>(arg).is-positive()) V ex-to<numeric>(arg).is-zero()) 

return numeric(l); 
else 

return numeric(-l); 

> 

return jump-fnct(arg) .holdf); 

} 


Defines: 

ex, used in chunks 2-8, 11a, 13, 14, 17, 19-25, 33b, 38-40, 42, 44-47, 49-60, 63, 65a, 67, and 68. 

Uses jump_fnct 38a and numeric lOd. 

66b (cycle.cpp 40a) += «66a 66c > 

static ex jump-fnct-eval( const ex & arg) 

{ 

if (■ is-exactly-a<mimev\c>(arg )) { 

if ((ex-to< numeric >(arg).is.real() A ex-to<numeric>(arg).is-positive( )) V ex-to<numeric>(arg).is-zero()) 

return numeric(l); 
else 

return numeric(-l); 

} else if ( is-exactly-a<mul>(arg ) A 

is-exactly_a< numerio( arg.op{arg.nopsQ- 1))) { 
numeric oc = ex-to<numeric>(arg.op(arg.nops()-l )); 
if ( oc.is-realQ) { 
if (oc > 0) 

// jump_fnct(42*x) -> jump_fnct(x) 
return jump-fnct(arg-h-oc).hold (); 

else 

// jump_fnct(-42*x) -> -jump_fnct(x) 
return -jump_fnct(arg 4- oc).hold(); 

} 

} 

return jump-fnct(arg).hold(); 

} 


Defines: 

ex, used in chunks 2-8, 11a, 13, 14, 17, 19-25, 33b, 38-40, 42, 44-47, 49-60, 63, 65a, 67, and 68. 
Uses jump_fnct 38a and numeric lOd. 

66c (cycle.cpp 40a) += <66b 67a > 

static ex jump-fnct-Conjugate (const ex & arg) 

{ 

return jump-fnct(arg)', 

} 


Defines: 

ex, used in chunks 2-8, 11a, 13, 14, 17, 19-25, 33b, 38-40, 42, 44-47, 49-60, 63, 65a, 67, and 68. 
Uses jump_fnct 38a. 


<66c 67b > 


67a (cycle.cpp 40a) += 

static ex jump-fnct-power) const ex & arg , const ex & exp) 

{ 

if ( js_a<numeric>(ea:p) A e*_to<numeric>(ea;p). is-integer{)) { 
if ( ex-to< numeric > ( exp) .is-even()) 

return numeric(l); 
else 

return jump-fnct{arg)\ 

} 

if ( is_a<numeric>(ea;p) A ex-to<numeric>(-exp).is-positive()) 
return exJ,o<has\c>)pow)jump-fnct)arg ), -exp)), hold)); 
return exJ.o<has\c>)pow)jump-fnct)arg) 1 exp)).hold)); 


Defines: 

ex, used in chunks 2-8, 11a, 13, 14, 17, 19-25, 33b, 38-40, 42, 44-47, 49-60, 63, 65a, 67, and 68. 
Uses jump_fnct 38a and numeric lOd. 

67b (cycle.Cpp 40a) += «67a 67c o 

static void jump-fnct-print-dflt_text)const ex & x, const print-context & c) 

{ 

c.s C "H("; x.print(c ); c.s 

} 

static void jump-fnct-printJatex( const ex & a;, const print-context & c) 

{ 

c.s <C "\\chi("; x.print(c ); c.s <C 

> 


Defines: 

jump_fnct_print_df It .text, used in chunk 67c. 
jump_fnct_print_latex, used in chunk 67c. 

Uses ex 4a lib 39d 48a 48b 66a 66b 66c 67a. 

All above methods are used to register the function now. 

67c (cycle.cpp 40a) += <67b 67do 

REGISTER_FUNCTION(jump_fnct , evaLfunc{jumpJncEeval). 
evalf-func{jump-fnct-evalf). 
latex-name( " \\ chi"). 

//text_name(” H”). 

print-func<print-dflt> (jump-fnct-print-dflt_text). 
print-func<printJatex> ( jump-fnct-printdatex). 

/ / derivative _func(2*delta). 
power_func(jump-fnct-power). 
conjugate-func(jump_fnct-Conjugate )); 


Uses jump_fnct 38a, jump_fnct_print_dflt.text 67b, and jump_fnct_print_latex 67b. 

This function prints if its parameter is zero in a prominent way. 

67d (cycle.cpp 40a) += «67c 68a o 

const char * equality) const ex & E) 

{ 

if ( normal(E). is-zero {)) 

return "-equal-"; 
else 

return "DIFFERENT!!!"; 

> 


Defines: 

equality, used in chunk 13. 

Uses ex 4a lib 39d 48a 48b 66a 66b 66c 67a. 


<67d 68b > 


This function decodes metric sign into human-readable form. 

68a (cycle.cpp 40a) += 

const char *eph_case( const numeric & sign) 

{ 

if (numeric(s« 5 ri-(-l)).is_zero()) 

return "Elliptic case (sign = -1)"; 
if (numeric(sign). is-zeroQ) 
return "Parabolic case (sign = 0)"; 
if (numeric(sign-l).is-zero()) 
return "Hyperbolic case (sign = 1)"; 
return "Unknown case!!!!"; 

} 


Defines: 

eph.case, used in chunk 20f. 

Uses numeric lOd. 

Elements of SL 2 (S>) are transformed into appropriate “cliffordian” matrix. 

68b (cycle.cpp 40a) += «68a 68c > 

matrix sl2-difford( const ex & a, const ex & b, const ex & c, const ex & d, const ex & e, bool not-inverse ) 
{ 

if (is_a<clifford>(e)) { 

ex eO = e.subs(e.op( 1) = 0); 

ex one = dirac-ONE(ex-to<c\if¥ord> (e). get_representationJ,abel{))\ 
if ( not_inverse ) 
return matrix(2, 2, 

1st (a * one, b * eO, 
c * pow{e0, 3), d * one)); 

else 

return matrix(2, 2, 
lst(d * one, -b * eO, 

-c * pow(e0, 3), a * one)); 

} else 

throxv(std::invalid-arguTnent("sl2_clif ford() : expect a Clifford numeber as a parameter")); 

} 


Uses bool lie, ex 4a lib 39d 48a 48b 66a 66b 66c 67a, 1st lie, and matrix 8c lOe Ilf 12a 25a. 

68c (cycle.Cpp 40a) += «68b 

matrix sl2-clifford(const ex & M, const ex & e, bool not-inverse) 

{ 

if (zs_a<matrix>(M) V M.info(info_flags::list)) 

return sl2-clifford{M.op( 0), M.op{ 1), M.op{ 2), M.op{ 3), e, notJnverse); 

else 

throw(std::invalid-argument("sl2_clif ford() : expect a list or matrix as the first parameter")); 

> 


Uses bool lie, ex 4a lib 39d 48a 48b 66a 66b 66c 67a, and matrix 8c lOe Ilf 12a 25a. 
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Appendix E. Index of Identifiers 

bool: 3c, 4c, 5c, 5d, 5e, 7b, 7c, 7e, 7f, 8a, lie, 12d, 21f, 38a, 39a, 46c, 47a, 50a, 54b, 54c, 54d, 55a, 56c, 56d, 57d, 58c, 62b, 
68b, 68c 

catch: 10b, 25d, 26b, 26c, 42d, 43b, 50a 
check_conformality: 19c, 19e, 20c . 22c, 23b 

cycle: 2a, 2b, 2c, 2d, 3c, 3d, 3e, 3e, 3e, 3e, 3e, 4a, 4b, 4f, 5c, 5d, 5e, 5f, 6b, 6c, 6d, 7a, 7b, 8d, 9e, 11c, 12e, 12f, 13b, 13c, 13d, 

14a, 14b, 15a, 16b, 16c, 17a, 17d, 17e, 17f, 20e, 24a, 24b, 24c, 24d, 24e, 24f, 25c, 25d, 25e, 33a, 38b, 39a, 39b, 40a, 40b, 40c, 

40d, 41b, 41c, 41d, 42a, 42b, 42c, 42d, 43a, 43b, 44b, 45a, 45b, 45c, 46a, 46b, 46c, 47a, 47b, 47b, 47c, 48a, 48a, 48a, 48a, 

48a . 49a, 49b, 50a, 51a, 51b, 52a, 53a, 53b, 54a, 54b, 54c, 54d, 55a, 55b, 55c, 56a, 58b, 58c 

cycle2D: 6f, 6g, 7a, 11c , 11c , 12c, 14b, 15e, 16b, 16c, 16d, 17f, 18g, 20a, 21g, 24a, 25b, 30a, 32a, 33a, 33b, 34b . 35a, 35b , 35c, 
36b, 36c, 36e, 37b, 39b, 39c, 39d, 39d, 39d, 39d, 39d, 40a, 41d, 48b, 48b, 48b, 48b, 48b, 55d, 55e, 55f, 56a, 56b, 56b, 56c, 
56d, 57a, 57d, 62c, 63a, 64c 
CYCLELIB_MAJ0R_VERSI0N: 37c 
CYCLELIB_MIN0R_VERSI0N: 37c 
debug: lie, 14d, 15d, 17g, 18c, 18g, 20a, 21f 
DRAW.ARC: 57c, 65b 
du: 57c, 62b 
dv: 57c, 62b 
eph.case: 20f, 38a , 68a 
eph_names: lOd . 30a, 32a, 36d, 37b 
equality: 13b, 13d, 38a . 67d 

ex: 2a, 2b, 2c, 2d, 3a, 3b, 3c, 3d, 3e, 4a, 4b, 4c, 4d, 4e, 4f, 4g, 5a, 5b, 5c, 5d, 5e, 5f, 6a, 6b, 6c, 6d, 6f, 6g, 7a, 7b, 7c, 7d, 7e, 
7f, 8a, 11a, lib, 13e, 13f, 14a, 14b, 17b, 17c, 17d, 17e, 17f, 19a, 19b, 20c, 21d, 22a, 22g, 23d, 23g, 24e, 25b, 33b, 38a, 38b, 
39a, 39c, 39d, 40b, 42a, 42b, 42c, 42d, 44a, 44b, 45a, 45b, 46a, 47c, 48a, 48b, 49a, 49b, 50a, 50b, 51a, 51b, 52a, 52b, 53a, 
53b, 54a, 54b, 54c, 54d, 55a, 55b, 55c, 55e, 55f, 56a, 56b, 56c, 56d, 57a, 57b, 57d, 58b, 58c, 59b, 60b, 60c, 63b, 63c, 65a, 66a, 
66b . 66c . 67a . 67b, 67d, 68b, 68c 
f ocal_length_check: 22d, 22e, 22f, 22g 
hypjnatr: 9e, 35b 

infinitesimal.calculations: 23e, 23f, 23g 

jump.fnct: lOe, 15e, 15f, 15g, 18c, 18e, 18f, 32a, 38a . 50b, 56c, 66a, 66b, 66c, 67a, 67c 

jump_fnct_print_df lt.text: 67b , 67c 
jump_fnct_print_latex: 67b . 67c 
k_d: 57c, 62b . 63b, 65a 

1st: 2b, 6g, 7b, 7c, 7d, 7e, 7f, 8a, 8b, 9e, lOe, lib, 11c, lid, lie, Ilf, 12c, 12d, 14b, 14c, 14e, 14g, 15b, 15d, 15e, 15f, 15g, 
16b, 16c, 16d, 17f, 17g, 18a, 18b, 18c, 18d, 18e, 18f, 18g, 19a, 19c, 19d, 19e, 20a, 20b, 21a, 21d, 22a, 22b, 22g, 23f, 24a, 24b, 
24d, 25a, 25c, 25d, 26a, 30a, 30b, 32a, 33a, 33b, 34a, 34b, 35a, 35b, 35c, 36a, 36b, 36c, 36e, 37a, 37b, 40c, 42a, 43a, 45b, 
45c, 46a, 50b, 51a, 52a, 52b, 53a, 54c, 55f, 56b, 56c, 56d, 57a, 57d, 58b, 58c, 60b, 61c, 62c, 63a, 63c, 64a, 64c, 68b 
main: 10a 

math.string: 9f, 12e, 13e, 13f, 14a, 14b, 14d, 15c, 15d, 17a, 17b, 17c, 17d, 17e, 17f, 17g, 18c, 19d, 20b, 21f, 21g, 22c, 22g, 24a, 
24b, 24c, 24d, 24e, 24f, 25c, 25d, 25e 

matrix: 2d, 4g, 5b, 8c, lOe, Ilf, 12a, 16d, 18b, 25a, 25d, 38a, 40c, 40d, 41a, 41b, 42c, 43a, 43b, 44a, 50b, 51b, 52a, 52b, 53a, 
54a, 54b, 54c, 54d, 55a, 55b, 56c, 68b, 68c 

numeric: 4b, 5a, 6b, lOd , 11c, 18g, 20a, 21a, 30a, 30b, 31a, 32a, 32b, 34a, 35c, 38a, 39c, 41d, 42a, 42d, 43b, 46c, 47c, 49a, 
50a, 50b, 51a, 56d, 57a, 57b, 57c, 58a, 58c, 59c, 60b, 60d, 61a, 61c, 62b, 63b, 63c, 64a, 65a, 65b, 66a, 66b, 67a, 68a 
PAIR: 57c , 60a, 61d, 62a, 64b 
par_matr: 9e, 33b, 34b, 35b 
points: lOd, 11a, 24b, 30a, 32a, 65b 
PRINT.CYCLE: 40a, 47b 
print.perpendicular: 19c, 19e, 21g, 22c 
realsymbol: lOd . lOe 


si: lOe . 16b, 20c, 20d, 20f, 20g, 21a, 21c, 21f, 26a, 30a, 30b, 31a, 31c, 32a, 36a, 36b, 36d, 37a, 37b 
sil: 10e . 20f, 21a, 21c, 26a, 30a, 31c, 32a, 36a, 36b, 36d 

u: lOd, 11a, 11c, Ilf, 15d, 16d, 19b, 19d, 19e, 20b, 21d, 21g, 22b, 22c, 22d, 22e, 22f, 22g, 23a, 23b, 23c, 24a, 24b, 25d, 30a, 
30b, 31a, 31c, 32a, 33a, 33b, 34a, 34b, 35a, 35b, 36e, 57c, 58b, 60a, 60c, 60d, 61a, 61b, 61c, 61d, 62a, 62b , 64b 
v: lOd, 11a, 11c, Ilf, 15d, 16d, 19b, 20a, 20b, 21d, 21g, 22b, 22c, 22d, 22e, 22f, 23a, 23b, 23c, 23d, 23g, 30a, 30b, 32a, 34a, 
39a, 57c, 58b, 60c, 61c, 61d, 62b 

varidx: lOd , 40c, 41a, 41b, 41c, 42c, 44a, 50a, 50b, 51b, 53a, 56c, 56d, 57a, 58c 

wspaces: 9f, 12e, 13e, 13f, 14a, 14b, 16b, 16c, 17b, 17c, 17d, 17e, 17f, 18g, 19d, 20a, 20b, 20d, 20e, 20f, 21f, 21g, 22c, 23a, 
23b, 23c, 24e, 24f, 25e 
zero_or_one: 57c, 62b . 65a, 65b 


School of Mathematics, University of Leeds, Leeds LS2 9JT, UK 
E-mail address: kisilv@maths.leeds.ac.uk 
URL: http://maths.leeds.ac.uk/”kisilv/ 



